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I. 


INTRODUCTION 


This  manual  describes  the  fourth  generation  of  SHABERTH, 

• "A  Computer  Program  for  the  Steady  State  and  Transient  Analyses 
of  Shaft  Bearing  Systems  Including  Ball  Cylindrical  and  Tapered 
Roller  Bearings."  The  original  program,  AE72Y003,  was  developed 
by  Kellstrom  (1)  under  U.S.  Army  Contract  DAAD05-73-C-0011, 
sponsored  by  the  Ballistics  Research  Laboratory  at  Aberdeen 

« Proving  Grounds.  The  present  version  was  developed  under  U.S. 

Army  Contract  DAAD05-75-C-0747 , also  sponsored  by  the  Ballistics 
Research  Laboratories. 

m 

The  master  program  consists  of  the  following  major  sub- 
programs : 

1)  Bearing  Analysis.  These  subprograms  are  largely 
based  upon  tne  methods  of  Harris,  (2,3). 

2)  Three  Demensional  Shaft  Deflection  Analysis  de- 
veloped by  Noriander  and  Friedrichson. 

* 

3)  Bearing  Dimensional  Change  Analysis  based  on  the 
methods  of  Timoshenko,  (4) , and  adapted  to  the 
shaf t-bearing-housing  system  by  Crecelius,  (5). 

4)  Generalized  Steady  State  and  Transient  Temperature 
Mapping  and  Heat  Dissipation  Analyses  based  on  the 
methods  of  Harris,  (6),  Fernlund,  (7)  and  Andreason, 
(8). 

Although  the  primary  function  of  all  four  generations 
of  the  program  is  to  predict  general  bearing  performance  char- 
acteristics, and  the  bulk  of  the  coding  reflects  this  emphasis, 
the  steady  state  and  transient  heat  dissipation  and  temperature 
mapping  sub-program  may  be  used  on  a stand  alone  basis  to  model 
the  thermal  behavior  of  any  system  which  can  be  represented  by 
discrete  temperature  nodes. 

The  differences  between  the  successive  generations  of 
the  program  reflect  the  development  and  installation  of  im- 
proved bearing  lubrication  and  friction  models,  improved  analy- 
sis of  the  bearing  cage  and  improvements  in  the  program  struc- 

• ture  which  increased  the  program  versatility  and  solution  pro- 
cedures. In  the  fourth  generation,  program  capabilities  were 
expanded  to  include  tapered  and  cylindrical  roller  bearings 

• with  flanges. 

The  first  generation  of  the  program  used  the  Newtonian 
lubricated  friction  models  developed  by  Harris  (2,3). 
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The  second  generation  of  the  program  which  carries  the 
designation  AT74Y001  was  created  by  Crecelius,  Liu  and  Chiu 
under  Air  Force  Contract  No.  F33615- 72-01467  and  Navy 
MIPR  No.  M52376- 3- 000007  and  is  documented  by  McCool,  et  al 
(9)*.  In  this  effort,  with  Program  AE72Y003as  the  basis, 
the  ball  bearing  subprogram  was  modified  to  include  new  models 
as  follows: 

1)  An  EHD  film  thickness  model  that  accounts  for  , 

i)  thermal  heating  in  the  contact  inlet  using  a 

regression  fit  to  results  obtained  by  Cheng  (10) 

and  ii)  lubricant  film  starvation  using  theoretical 

results  derived  by  Chiu  (11)  . *- 

2)  A new  semi- empirical  model  for  fluid  traction  in  an 
EHD  contact  (9) , is  combined  with  an  asperity  load 
sharing  model  developed  by  Tallian  (12)  to  yield  a 
model  for  traction  in  concentrated  contacts  which 
reflects  various  conditions  of  lubrication.  Dry, 
partial  and  full  EHD  regimes  are  considered. 

3)  A model  for  the  hydrodynamic  rolling  and  shear  forces 
in  the  inlet  zone  of  lubricated  contacts  accounting 
for  the  degree  of  lubricant  film  starvation,  (9). 

4)  Normal  and  friction  forces  between  a ball  and  a cage 
pocket  are  modelled  in  a way  that  accounts  for  the 
transition  between  the  hydrodynamic'  and  elasto- 
hydrodynamic  regimes  of  lubrication  (9) . 

5)  A model  for  the  effect  on  fatigue  lirfe  of  the  ratio 
of  the  EHD  plateau  film  thickness  to  the  composite 
surface  roughness,  (9). 

Additionally,  models  for  temperature  viscosity  and  pressure 
viscosity  variation  as  functions  of  temperature  given  by  Walther 
(13)  and  Fresco  (14)  respectively,  were  adopted. 

Program  AT74Y001  is  capable  of  analyzing  only  a single 
axially  loaded  ball  bearing.  The  program  cannot  be  used  to 
analyze  a multi-bearing  system.  All  other  capabilities  are 
present  however. 

• 

“Due  to  the  similarities  between  major  segments  of  SHABERTH  AT77Y001  5 
AT74Y001 , many  sections  of  (9)  have  been  included  in  this  f 

text,  without  modification. 
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The  basis  for  the  third  version  AT75Y004,  was  AT74Y001. 

The  latent  capability  for  the  analysis  of  up  to  five  ball  and 
roller  bearings  subjected  to  general,  (5  degrees  of  freedom) 
loading,  was  utilized.  The  models  added  to  AT74Y001  were  used 
in  the  calculation  of  both  the  ball  and  roller  bearing  friction 
forces  and  frictional  heat  generation  rates.  The  third  version 
also, includes  a model  for  the  hydrodynamic  rolling  and  slip 
forces  in  the  inlet  zone  of  lubricated  line  contacts,  based  on 
the  work  of  Chiu,  Ref.  (15).  Additionally,  a cage  model  deve- 
loped under  NASA  Contract  No,  NAS3-19739  Ref.  (15)  was  added 
which  allows  the  cage  to  move  with  up  to  three  degrees  of 
freedom  versus  the  one  degree  of  freedom  permitted  in  Program 
AT74Y001.  This  cage  model  is  used  in  the  analysis  of  both 
ball  and  roller  bearings. 

This  fourth  generation  program,  SHABERTH,  AT77Y001,  differs 
from  the  preceding  versions  in  that  it  can  accomodate  tapered 
roller  bearings  as  well  as  ball  and  cylindrical  roller  bearings. 
Additionally,  the  capability  of  the  cylindrical  roller  bearing 
analysis  has  been  expanded  to  permit  load  to  be  carried  by 
inner  and  outer  ring  flanges. 

SHABERTH  is  intended  to  be  as  general  as  possible  with  the 
following  limits  on  system  size. 

Number  of  bearings  supporting  the  shaft  - five  (5) 
maximum 

Number  of  rolling  elements  per  bearing  - thirty  (30) 
maximum 

Number  of  temperature  nodes  used  to  describe  the  system- 
one  hundred  (100)  maximum 

The  program  structure  is  modular  and  has  been  designed  to 
permit  substitution  of  new  mathematical  models  and  refinements 
to  the  existing  models  as  the  needs  and  opportunities  develop. 
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2.  PROBLEM  FORMULATION  AND  SOLUTION 

The  purpose  of  the  program  is  to  provide  a tool  with  which 
the  shaft  bearing  system  performance  characteristics  can  be 
determined  as  functions  of  system  temperatures.  These  system 
temperatures  may  be  a function  of  steady  state  operation  or  a 
function  of  time  variant  conditions  brought  on  by  a change 
in  the  system  steady  state  condition.  Such  a change  would  be 
the  termination  of  lubricant  supply  to  the  bearings  and  other 
lubricated  mechanical  elements. 

The  program  is  structured  with  four  nested,  calculation 
schemes  as  follows: 

1.  Thermal',  steady  state  or  transient  temperature  calcu- 
lations which  predict  system  temperatures  at  a given 
operating  state. 

2.  Bearing  dimensional  equilibrium  which  uses  the  bearing 
temperatures  predicted  by  the  temperature  mapping 
subprograms  and  the  rolling  element  raceway  load  distri- 
bution, predicted  by  the  bearing  subprograms,  to  calcu- 
late bearing  diametral  clearance  at  a given  operating 
state. 

3.  Shaft-bearing  system  load  equilibrium  which  calculates 
bearing  inner  ring  positions  relative  to  the  respective 
outer  rings  such  that  the  external  loading  applied  to 
the  shaft  is  equalized  by  the  rolling  element  loads 
which  develop  at  each  bearing  inner  ring  at  a given 
state. 

4 . Bearing  rolling  element  and  cage  load  equilibrium 
which  calculates  the  rolling  element  and  cage  equili- 
brium  positions  and  rotational  speeds  based  upon  the 
relative  inner-outer  ring  positions,  inertia  effects 
and  friction  conditions,  which  if  lubricated,  are 
temperature  dependent. 

The  above  program  structure  allows  complete  mathematical 
simulation  of  the  real  physical  system.  The  program 
has  been  coded  to  allow  various  levels  of  program  execution 
which  prove  useful  and  economical  in  bearing  design  studies. 

These  levels  of  execution  are  explained  fully  in  Sections 

3,  4,  and  5. 
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The  structure  of  the  program  and  the  nesting  of  the  solu- 
tion loops  noted  above  can  be  seen  clearly  in  the  Program  Flow 
Chart  which  is  discussed  in  Appendix  A. 

The  sections  below  present  the  systems  of  field  equations 
which  are  solved  in  each  of  the  nested  calculation  schemes. 

A more  detailed  discussion  is  contained  in  (1,9  and  15). 

2 . 1 Temperature  Calculations 

Subsequent  to  each  calculation  of  bearing  generated  heat 
rates,  either  the  steady  state  or  transient  temperature  mapping 
solution  scheme  may  be  executed.  This  set  of  sequential  calcu- 
lations is  terminated  as  follows: 

1.  For  the  steady  state  case;  when  each  system  temperature 
is  within  EPI  °Centrigrade  of  its  previously  pre- 
dicted value,  EPI  is  specified  by  the  user.  If  it  is 
zero  or  left  blank,  a default  value  of  1°  Centigrade  is 
used.  This  criteria  implies  that  the  steady  state 
equilibrium  conditions  has  been  reached. 

2.  The  transient  calculation  terminates  when  the  user 
specified  time  up  is  reached  or  when  one  of  the  system 
temperatures  exceeds  600°C. 

2.1.1  Steady  State  Temperature  Map 

The  mechanical  structure  to  be  analyzed  is  thought  of  as 
divided  into  a number  of  elements  or  nodes,  each  represented  by 
a temperature.  The  net  heat  flow  to  node  i from  the  surrounding 
nodes  j,  plus  the  heat  generated  at  node  i,  must  numerically 
equal  zero.  This  is  true  for  each  node  i,  i going  from  1 to  n, 
n being  the  number  of  unknown  temperatures. 

After  each  calculation  of  bearing  generated  heat,  which 
results  from  a solution  of  the  shaft-bearing  system  portion 
of  the  program,  a set  of  system  temperatures  is  determined 
which  satisfy  the  system  of  equations: 

q^  = qQ^  + qg^  =,  0 for  all  temperature  nodes  i (2.1) 

where  qQ^  is  the  heat  flow  from  all  neighboring  nodes  to 
node  i 

q . is  the  heat  generated  at  node  i.  These  values 
g may  be  input  or  calculated  by  the  shaft  bearing 
program  as  bearing  frictional  heat 
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This  scheme  is  solved  with  a modified  Newton- Raphson 
method  which  successfully  terminates  when  either  of  two 
condition  are  met: 


A*-i  < EP2  for  all  nodes  i (2.2) 


where  : At  represents  the  New  ton- Raphson  correction  to  the 
temperature  t at  a given  iteration  such  that, 
tN+l  = tN  + At  and  N + 1,  and  N,  refer  to  the 
next  and  current  iteration  respectively. 


EP2  is  a user  specified  constant.  If  EP2  is  left  blank 
or  set  to  zero  (0)  a default  value  of  0.001  is  used. 


A second  convergence  criterion  dependent  upon  EP2  is  also 
used.  In  the  system  of  equations,  q0^  + qgi  = 0 for  all  nodes 
i,  absolute  convergence  would  be  obtained  if  the  right  hand 
side  (EQ)  in  fact  reduced  to  zero  (0).  Usually  a small  residue 
remains  at  each  node,  such  that  (qoi  + Qgi)  = (EQ)i* 

The  second  convergence  criterion  is  satisfied  if: 


n 


l 

l 


(2.3) 


where  n = number  of  equations  in  bearing  solution 


2.1.2  Transient  Temperatures 

In  the  transient  case  the  net  heat  q^  transferred  to  a 
node  i heats  the  element.  It  is  thus  necessary  for  heat  balance 
at  node  i that  the  following  equations  are  satisfied. 


dt  i 

pi  Cpi  Vi  ITT  = qi 


(2.4) 


< 
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where  p - density 

Cp  = specific  heat 

V=  volume  of  the  element 
t=  temperature 
T=  time 

The  temperatures,  t0^,  at  the  time  of  initiation  T = Ts  are 
assumed  to  be  known,  that  is 

ti (Ts)  = t q i i = 1 , 2,  »»•>  n (2.5) 

The  problem  of  calculating  the  transient  temperature 
distribution  in  a bearing  arrangement  thus  becomes  a problem 
of  solving  a system  of  non-linear  differential  equations  of 
the  first  order  with  certain  initial  values  given.  The 
equations  are  non-linear  since  they  contain  terms  of  radiation 
and  free  convection,  which  are  non-linear  with  temperature  as 
will  be  shown  later.  The  simplest  and  most  economical  way  of 
solving  these  equations  is  to  calculate  the  rate  of  temperature 
increase  at  the  time  T = T^  from  equation  2.4  and  then  calculate 
the  temperatures  at  time  T^  + AT  from 

‘k*1  = iT  ' *k  * pwv  dT  (2-6) 


If  the  time  step  AT  used  as  program  input  is  chosen 
too  large,  the  temperatures  will  oscillate,  and  if  it  is  chosen 
too  small  the  calculation  will  be  costly.  It  is  therefore 
desirable  to  choose  the  largest  possible  time  step  that  does 
not  give  an  oscillating  solution.  The  program  optionally  cal- 
culates such  a time  step.  The  step  is  obtained  from  the 
condition,  {16} 

dti.k+l  >0  i ■ 1,  2 n (2.7) 

dti,k 

If  this  derivative  were  negative,  the  implication  would 
be  that  the  local  temperature  at  node  i has  a negative  effect 
on  its  future  value.  This  would  be  tantamount  to  asserting 
that  the  hotter  a region  is  now,  the  colder  it  will  be  after 
an  equal  time  interval.  An  oscillating  solution  would  result. 
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Differentiating  equation  (2.6)  for  node  i,  one  obtains 


dt 

dt 


i,k+l 
i ,k 


1 + 


ATi- 


dq 
* dt 


i 

i >k 


i = 1 , 2 , . . .n 


The  derivative  dq/dt^  is  calculated  numerically 
dqi,k  = qi ,k+l  - qitk 


(2.8) 


(2.9) 


For  each  node,  the  value  of  AT  is  calculated  which  gives 
a value  of  zero  to  the  right  hand  side  of  Eqn.  (2.8)  and  thus 
Eqn.  (2.7)  as  well. 

A value  of  T rounded  off  to  one  significant  digit  smaller 
than  the  AT^  given  by  Eqn.  (2.8)  is  used. 

2.1.3.  Calculation  of  Heat  Transfer  Rate 


The  transfer  of  heat  within  a medium  or  between  two  media 
can  occur  by  conduction,  convection,  radiation  and  fluid  flow. 

All  these  types  of  heat  transfer  occur  in  a bearing  appli- 
cation as  the  following  examples  show. 

1.  Heat  is  transferred  by  conduction  between  inner  ring 
and  shaft  and  between  outer  ring  and  housing. 

2.  Heat  is  transferred  by  convection  between  the  surface 
of  the  housing  and  the  surrounding  air. 

3.  Heat  is  transferred  by  radiation  between  the  shaft 
and  the  housing. 

4.  When  the  bearing  is  lubricated  and  cooled  by  circulating 
oil,  heat  is  transferred  by  fluid  flow. 

Therefore,  in  calculating  the  net  flow  to  a node  all  the 
above  mentioned  modes  of  heat  transfer  will  be  considered. 

2 . 1 . 3 . 1 Generated  Heat  « 

There  may  be  a heat  source  at  node  i giving  rise  to  a heat 
flow  to  be  added  to  the  heat  flowing  from  the  neighboring  nodes. 
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In  the  case  that  the  heat  source  is  a bearing,  it  may 
either  be  considered  to  produce  a known  amount  of  heat,  in  which 
case  constant  numbers  are  entered  as  input  to  the  program,  or 
the  shaft-bearing  program  may  be  used  to  calculate  the  bearing 
generated  heat  as  a function  of  bearing  temperatures. 

2 . 1 . 3 . 2 Conduction 

The  heat  flow  q • . which  is  transferred  by  conduction 
from  node  i to  node  J,  ’is  proportional  to  the  difference  in 
temperature  (t^  - tj)  and  the  cross  sectional  area  A and 
is  inversely  proportional  to  the  distance  i between  the  two 
points,  thus 


^ CV  tp  (2.10) 

where  X-  the  thermal  conductivity  of  the  medium. 

2. 1.3.3  Free  Convection 


Between  a solid  medium  such  as  a metallic  body  and  a 
liquid  or  gas,  heat  transfer  is  by  free  or  forced  convection. 
Heat  transfer  by  free  convection  is  caused  by  the  setting  in 
motion  of  the  liquid  or  gas  as  a result  of  a change  in  density 
arising  from  a temperature  differential  in  the  medium.  With 
free  convection  between  a solid  medium  and  air,  the  heat 
energy  q . . transferred  between  nodes  i and  i can  be  calculated 
from  the  elation,  (2.11) 

qvi,j  = “v*  Cti-tj^  d*  SIGNCt.-tj)  (2.11) 


where 


A 

d 


the  film  coefficient  of  heat  transfer  by  free 
convection 

the  surface  area-  of  contact  between  the  media 
is  an  exponent,  usually  = 1.25,  but  any  value 
can  be  specified  as  input  to  the  program 


SIGN 


f1  if  tj 

v;1  if  V 'j 


» 
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The  last  factor  is  included  to  give  the  expression  qvi^.j 
a correct  sign. 

The  value  of  Oy  can  be  calculated  for  various  cases,  see 
Jacob  and  Hawkins,  (16). 

2. 1.3. 4 Forced  Convection 

Heat  transfer  by  forced  convection  takes  place  when  liquid 
or  gas  moves  around  a solid  body,  for  example,  when  the  liquid 
is  forced  to  flow  by  means  of  a pump  or  when  the  solid  body  is 
moved  through  the  liquid  or  gas.  The  heat  flow  q . . transferred 

by  forced  convection  can  be  obtained  from  the  following  equation. 


Iwi.j  =aw  A(li  - V 


(2.12) 


where  a is  the  film  coefficient  of  heat  transfer  during 
forced  convection.  This  value  is  dependent  on 
the  actual  shape,  the  surface  condition  of  the 
body,  the  difference  in  speed,  as  well  as  the 
properties  of  the  liquid  or  gas. 

In  most  cases,  it  is  possible  to  calculate  the  coefficient 
of  forced  convection  from  a general  relationship  of  the  form, 


Nu  = aRebPrc 


(2.13) 


where  a,  b,  and  c are  constants  obtained  from  handbooks 

such  as  (17).  R and  P are  dimensionless  numbers 
defined  by  e 


Nu  = Nusselt  number  = aw  L/A 

L = characteristic  length 
A = conductivity  of  the  fluid 
R = Reynold's  number  = ULp/n 

U = characteristic  speed 
p = density  of  the  fluid 
ri  = dynamic  viscosity  of  the  fluid 
Pr  = Prandtl  number  = nCp/A 

Cp  = specific  heat 


< 
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The  program  can  use  a value  of  the  coefficient  of  convec- 
tion, or  let  it  vary  with  actual  temperatures,  the  variation 
being  determined  by  how  the  viscosity  varies.  Input  can  be 
given  in  one  of  the  following  ways,  for  each  coefficient. 

Constant  viscosity 

1.  Values  of  the  parameters  of  equation  (2.13)  are  given 
as  input  and  a constant  value  of  is  calculated 

by  the  program. 

Temperature  dependent  viscosity 

2.  The  coefficient  for  turbulent  flow  and  heating  of 
petroleum  oils  is  given  by 

aw  - kg  . n(t)k10  (2.14) 

where  kg  and  k-^g  are  given  as  input  together  with  viscosity 
at  two  different  temperatures. 

3.  Values  of  the  parameters  of  equation  (2.13)  are  given 
as  input  together  with  viscosity  at  two  different 
temperatures . 

2 . 1 . 3 . 5 Radiation 

If  two  flat  parallel,  similar  surfaces  are  placed  close 
together  and  have  the  same  surface  area  A,  the  heat  energy 
transferred  by  radiaiton  between  nodes  i and  j representing 
those  bodies,  will  be, 

qRi,j  = eaA  (ti  + 273) 4 - (t j + 273) 4 (2.15) 

where  e is  the  surface  emissivity.  The  value  of  the 
coefficient  e is  an  input  variable  and  varies 
between  1 for  a completely  black  surface  and  0 
for  an  absolutely  clean  surface.  In  addition 
a is  Stefan-Bolt zmann ' s radiation  constant 
which  has  the  value  5.76  x 10“8  watts/m2- (°K) 4 
and  t^  and  tj  are  the  temperatures  (°C)  at  points 
i and  j . J 

Heat  transfer  by  radiation  under  other  conditions  can  also 
be  calculated,  (16) . The  following  equation,  for  instance 
applies  between  two  concentric  cylindrical  surfaces. 


11 


(tj  +273) 


(2.16) 


n . eA[cV273)4  ' 

R1,j  1 * (1-eT  (A~ /AeT 

where  is  the  area  of  the  inner  cylindrical  surface 
A is  the  area  of  the  outer  cylindrical  surface 

v 

2. 1.3. 6 Fluid  Flow 

Between  nodes  established  in  fluids,  heat  is  transferred 
by  transport  of  the  fluid  itself  and  the  heat  it  contains. 


€k  ♦ 

r*t- 

r i 

^ / / , 7 / / < 


i 1 

! ff)  s 1 

‘ « ‘-y- 

, 1 1 

L, y ■,  ! 

i— — ^ ® J ; 

‘ j 

7 — 7 — 7 — -/'// — 7 — 7" 


Figure  2.1  Convective  Heat  Transfer 


Figure  2.1  shows  nodes  i and  j at  the  midpoints  of  consecu- 
tive segments  established  in  a stream  of  flowing  fluid. 

The  heat  flow  through  the  boundary  between  nodes 

i and  j can  be  calculated  as  the  sum  of  the  heat  flow  q.^  through 

the  middle  of  the  element  i,  and  half  the  heat  flow  qQ^  trans- 
ferred to  node  i by  other  means,  such  as  convection. 

The  heat  carried  by  mass  flow  is, 

qfi  = pi  Cpi  Vi  ti  = Kiti  (2.17) 

where  V-^  = the  volume  flow  rate  through  node  i 
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The  heat  input  to  node  i is  the  sum  of  the  heat  generated 
at  node  i (if  any)  and  the  sum  over  all  other  nodes  of  the 
heat  transferred  to  node  i by  conduction,  radiation,  free  and 
forced  convection. 


(2.18) 


The  heat  flow  between  the  nodes  of  Fig.  2.2  is  then 


qui,j  " qfi  + qoi/2 

If  the  flow  is  dividing  between  node  i and  j, 
then  the  heat  flow  is  calculated  from 


(2.19) 
Figure  2.2 


q . . = K . . 

HU1,J  1J 


(qfi  * qoi/2) 


(2.20) 


where  K- . = the  proportion  of  the  flow  at  i going  to 
3 node  j,  0 <K. . <1.  K. - is  specified  at 

input.  3 3 


2. 1.3. 7 Total  Heat  Transferred 


The  net  heat  flow  rate  to  node  i can  be  expressed  as, 


qi  ~ qG,i  + j ^qci, j + qui , j + qvi,j  + qwi,j^ 


(2.21) 
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The  summation  should  include  all  nodes  j , both  with 
unknown  temperatures  as  well  as  boundary  nodes,  at  which  the 
temperature  is  known  so  long  as  they  have  a direct  heat  exchange 
with  note  i. 

This  expression  is  a non-linear  function  of  temperatures 
because  of  the  terms  qw  and  q^.  Therefore  the  equations  to  be 
solved  for  a steady  state  solution  are  non-linear. 


2.1.4  Conduction  Through  a Bearing 

As  described  in  Section  2. 1.3. 2 the  conduction  between 
two  nodes  is  governed  by  the  thermal  conductivity  parameter  X 
of  the  medium  through  which  conduction  takes  place.  The  value 
of  X is  specified  at  input. 

An  exception  is  when  one  of  the  nodes  represents  a . 
bearing  ring  and  the  other  a set  of  rolling  elements.  In  this 
case  the  conduction  is  separately  calculated  using  the  princi- 
ples described  below.  Note  that  separate  calculations  are 
performed  for  the  rolling  element  raceway  contacts  and  the  roll- 
ing element-flange  contacts.  The  methods  for  both  calculations 
are  identical  and  are  performed  within  the  program. 

2 . 1 . 4 . 1 . Thermal  Resistance 

It  is  assumed  that  the  rolling  speeds  of  the  rolling 
elements  are  so  high  that  the  bulk  temperature  of  the  rolling 
elements  are  the  same  at  both  the  inner  and  outer  races,  except 
in  a volume  close  to  the  surface.  The  resistance  to  heat  flow 
can  then  be  calculated  as  the  sum  of  the  resistance  across 
the  surface  and  the  resistance  of  the  material  close  to  the 
surface . 


The  resistance  Q is  defined  implicitly  by 

A t = fl*  q (2.22) 

where 

At  is  temperature  difference 
q is  heat  flow 

The  resistance  due  to  conduction  through  the  EHD  film  is 
calculated  as 


fli  = (h/  X )•  A 


(2.23) 


where  h is  taken  to  be  the  calculated  plateau  film  thickness 
A is  the  Hertzian  contact  area  at  the  specific  rolling 
element-ring  ocntact  under  consideration. 

X is  the  conductivity  of  the  oil. 

The  geometry  is  shown  in  Figure  2.3(a).  Asperity  conduction 
is  not  considered. 

So  far,  a constant  temperature  difference  between  the 
surfaces  has  been  assumed.  But  during  the  time  period  of 
contact,  the  difference  will  decrease  because  of  the  finite 
thermal  diffusivity  of  the  material  near  the  surface.  Fig.  2.3 
(b)  . 


Two  pnints  at  a distance  from  the  surface  this  phenomenon 
will  have  the  same  effect  as  an  additional  resistance  ^ acting 
in  series  with  . 

This  resistance  was  estimated  in  {18}  as, 

l r m/j-A 

ft2"  *Are,i  ^2F7v'  (2.24) 


where  £-re  = contact  length,  or  inthe  case  of  an 

elliptical  contact  area,  0.8  times  the 
major  axis 

A = heat  conductivity 

ip  = thermal  diffusivity  = A/(p.Cp) 

p = density 


Cp  = specific  heat 
b = half  the  contact  width 
V = rolling  speed 


The  resultant  resistance  is 


a 


res 


+ ^2 


TEMPERATURE  tj 


AVERAGE  FILM 
THICKNESS  - h 


direction  of  rolling 


(b)  Temperature  Distribution  at  Rolling, 
Concentrated  Contact  Surfaces 


FIGURE  2.3 

CONTACT  GEOMETRY  AND  TEMPERATURES 

t 
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There  is  one  such  resistance  at  each  rolling  element. 
They  all  act  in  parallel.  The  resultant  resistance ^tres , 
is  thus  obtained  from 


= l 1 (2.26) 

^ re's  i=l  “res , i 


2 . 2 Bearing  Dimensional  Change  Analysis 

The  program  calculates  the  changes  in  bearing  diametral 
clearances  according  to^  the  analysis  described  in  {5},  and 
expressed  in  generalized  equation  form  as, 

ADCL  = f { (Fits)m,  ti,  m,  CQr)m},  m = 1,2  for  inner  (2*27 

and  outer 

rings  respectively 

i = 1,2, 3,4, 5 for 
shaft,  inner  ring, 
outer  ring,  housing 
and  rolling  element 
respectively 
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where:  ADCL  is  the  change  in  bearing  diametral  cearance 
Fits  are  the  cold  mounted  shaft  and  housing  fits, 
t^  are  the  component  temperatures 

refers  to  the  ring  rotational  speeds 
Q™  refers  to  the  radial  component  of  the  minimum 
rolling  element-race  normal  force 


A bearing  clearance  change  criterion  is  satisfied  when 
the  change  in  bearing  diametral  clearance  remains  within 
a narrow,  user  specified  range,  for  two  successive  iterations 
as  follows: 

|(ADCL)n  - (ADCL)n_1  < EPSFIT  for  all  bearings  (2.28) 


D 

where:  N denotes  the  most  recent  iteration  and 
N-l  denotes  the  previous  iteration, 

D denotes  the  ball  or  roller  diameter  and 
EPSFIT  is  a user  specified  value,  = .0001D 

It  should  be  noted  that  although  ring  rotational  speeds, 
and  initial,  i.e.  cold,  shaft  and  housing  fits  are  considered 
in  the  clearance  change  analysis,  these  two  factors  are  fixed 
at  input  and  remain  constant  through  the  entire  solution. 

Although  component  temperatures  may  change  as  a consequence 
of  the  thermal  solution,  temperatures  remain  constant  through 
a complete  set  of  clearance  change  iterations.  As  a result,  only 
the  change  in  bearing  load  distribution  affects  the  change  in 
bearing  clearance  within  a set  of  clearance  change  iterations. 

2 . 3 Bearing  Inner  Ring  Equilibrium 

The  bearing  inner  ring  equilibrium  solution  is  obtained 
by  solving  the  .system: 

-+•  •* 

CFM^) i - (FMg)i  = 0 for  all  bearings,  i (2.29) 


where:  FM^,  denotes  a vector  of  bearing  loads  and 

moments  resulting  from  rolling  element/race 
forces  and  moments. 


Forces 


(2.30) 


FMbi  = 


bxi 


byi 

bzi 


M. 


byi 


M. 


bzi 


Moments 


If  the  bearing  solution  considers  friction,  FM.  is  com 
prised  of  the  ball  race  friction  forces  as  well  as  the 
normal  forces. 

If  the  beaming  solution  is,  at  the  user's  option, 
frictionless,  FM,  is  comprised  only  of  rolling  element/race 
normal  contact  forces. 

FM  • denotes  a similar  vector  of  loads,  exerted  on  the 
inner  ring  by  the  shaft. 


Forces 


(2.31) 

Moments 


FM  • 
si 


sxi 


syi 


szi 


M. 


syi 


M 


szi 


The  variables  in  this  system  of  equations  are  ths+bearing 
inner  ring  deflections  A,  and  the  shaft  displacements  A at 
all  bearing  locations.  The  bearing  loads  may  be  expressed  as 
a function  of  the  inner  ring  deflections. 

FMb‘ = FMb  (Ab)  (2.32) 
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The  deflection  Ab  of  a bearing  is  described  by  two  radial 
deflections  $ y and  5z,two  angular  deflections  €>  y and£>z  and 
axial  deflection  4 x*  The  axial  deflection  is  assumed  to  oe 
the  same  for  all  bearings . 

The  solution 
-*• 

6 (A)  • • 

J i.1  < 


If  for  some  i or 
(2.37) . 


NBRG  denotes 
EPS1  or  EPS2 

tions  are  fricitonless  or  include  friction,  respectively.  If 
the  bearing  deflections  are  extremely  small,  computer- generated 
numerical  inaccuracies  may  prevent  convergence  according  to  the 
above  criteria  although  a perfectly  good  solution  has  been  ob- 
tained. To  overcome  this  problem,  the  iteration  is  terminated 
if  all  angular  deflections  are  less  than  2 x 10-6  radians 
and  all  linear  deflections  are  less  than  5 x 10-1 2 * * * * * 8  inches.  Any 
one  of  the  above  criteria  imply  that  inner  ring  equilibrium  is 
satisfied. 

2 . 4 Bearing  Quasi-Dynamic  Solution 

The  bearing  quasi-dynamic  solution  is  obtained  through  a 
two  step  process: 

1)  Elastic  Solution  - considering  rolling  element  centri- 
fugal force,  plus  the  gyroscopic  moment  for  a tapered 
roller. 

2)  Elastic  and  Quasi-dynamic  Solution*  • 

*Quas i- dynamic  equilibrium  is  used  to  connote  that  the  true 

dynamic  equilibrium  terms  containing  first  derivatives  of  the 

ball  rotational  speed  vectors  and  the  second  derivatives  of 

rolling  element  position  vectors  with  respect  to  time  are 

replaced  by  numerical  expressions  which  are  position  rather 

than  time  dependent. 


(2.37) 


scheme  is  ended  when 
EPS1  (frictionless) 

EPS2  (friction) 

i = 1,...  (Number  of  bearings) 

j =1,5  for  the  3 linear  and  two 
angular  deflections  at 
each  bearing 

j,  (A)  * 0,  Eq  (2.38)  is  used  in  place  of 


EPS 1 (frictionless) 


(2.38) 


(0.001  x NBRG)  ESP2  (friction) 
the  number  of  bearings  in  the  system. 

is  used  depending  on  whether  the  bearing  solu- 
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The  equations  which  define  rolling  element  quasi-dynamic 
force  equilibrim  take  the  form  v 


zr  rm~ 

SU  «. + 

m 

->■ 

where :Q_  is  the 
m of  the 


m = 1-3  refers  to  the  outer 
inner  and  cage  rolling  element 
contacts  respectively 
m = 1-4  for  tapered  rollers  with 
-i  one  roller  end-flange  contact. 

+ F |+F=0  m =1-7  for  cylindrical  rollers 

m*J  where  four  roller  and  flange  (2.39) 

contacts  are  possible. 

vector  normal  load  per  unit  length 
contact.  See  Ref.  (1).  * 


f is  the  vector  of  friction  force  per  unit  length 
of  the  contact.  See  Ref.  (9). 


F is  the  vector  of  inertia  and  drag  forces.  See  Ref.(l) 
s is  a coordinate  along  the  contact 

perpendicular  to  the  direction  of  rolling  (usually 
the  major  axis) 

a is  half  the  contact  length.  See  Ref.  (1). 

Fm  is  the  vector  sum  of  the  hydrodynamic  forces 

acting  on  the  rolling  element  at  the  m-th  contact. 

For  ball-raceway  contact  see  Ref.  (9).  For  the 
roller- raceway  contact,  see  Ref.  (16). 


Z 

m 


Rolling  element  moment  equilibrium  is  defined  by: 


-a 


m 


• 

«C 

+ f 

->■ 

f 

F 

a 

m 

, m. 

m, 

-V 

r 

m 

m 


+ MT  = 


(2.40) 


and  s are  defined  above,  Mj  is  a vector 
of  inertia  moments.  For  the  definition  of  Mj 

refer  to  Ref . (1)  . 

is  a vector  from  the  rolling  element 
center  to  the  point  of  contact. 


In  the  frictionless  elastic  solution  F and  f = 0 

m m 

Additionally,  the  only  rolling  element  inertia  term  considered 
in  the  frictionless  solution  is  centrifugal  force,  plus  the  gyro- 
scopic moment  for  tapered  rollers.  As  a consequence  only  the 
axial  and  radial  force  equilibrium  equations  are  solved  for  each 
ball.  For  each  roller  the  radial  and  axial  force  equilibrium 
and  the  tilting  moment  about  the  z axis  of  Fig.  2.4  is  solved 
A dummy  equation  for  axial  force  equilibrium  is  included  in  the 
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FIGURE  Z.4 

BEARING  INERTIAL  (XYZ)  AND  ROLLING  ELEMENT  (xyz) , 
COORDINATE  SYSTEMS 
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solution  matrix  which  keeps  the  roller  centered  with  respect 
to  the  outer  race  if  the  cylindrical  roller  bearing  has  no  flanges 
* 

The  friction  solution  determines  ball  quasi-dynamic  equili- 
brium for  six  degrees  of  freedom.  A roller  is  permitted  four 
degrees  of  freedom.  The  rolling-element  variables  in  this 
solution  are  xl,  yl,  Wx,  Wy,  Wz,  and  Wo. 

where  x^  is  the  rolling  element  axial  position  relative 
to  the  outer  race. 

y^  is  the  rolling  element  radial  position  relative 
to  the  outer  race, 

Wx,  Wy,  Wz  are  orthogonal  rolling  element  rotational 
speeds  relative  to  the  cage  speed,  about  the  x,y,  and 
z axes  and  W0  is  the  rolling  element  orbital  speed. 

For  the  roller,  Wz  is  a dummy  variable. 

The  variables  Xi  and  y^  are  the  ball  variables  in  the  friction- 
less solution.  The  variables  in  the  roller  frictionless  solution 
are  x-^,  y^  and0z  = tan  - 1 (Wy/Wx) . 

The  details  of  the  cage  analysis  are  contained  in  SKF 
Report  No.  AL76P003.  Either  one  or  three  cage  equilibrium 
equations  are  considered  depending  upon  cage  type  and  the  degree 
of  rolling  element  orbital  speed  variation. 

The  cage  equations  and  cage-rolling  element  interactions 
are  not  considered  when  the  friction  forces  are  omitted  from 
the  rolling  element  equilibrium  equations. 

The  number  of  degrees  of  freedom  given  to  the  cage  is  one, 
if  the  cage  will  tend  to  rotate  concentrically  with  respect  to 
the  ring  on  which  it  is  riding.  This  condition  is  determined 
as  a function  of  the  rolling  element  orbital  speed  variation 
and  prevails  with  most  roller  bearings  and  with  ball  bearings 
subjected  only  to  axial  loading.  In  both  cases,  orbital  speed 
variation  is  often  inconsequential.  Also,  a single  degree  of 
freedom  is  allowed  when  the  cage  is  rolling  element  riding. 

The  single  degree  of  freedom  corresponds  to  a small  angular 
rotation  about  the  bearing  axis,  measured  with  respect  to  rolling 
element  1.  The  angular  displacement  is  converted  to  a linear 
dimension  by  a multiplication  by  the  bearing  pitch  diameter 
and  is  noted  in  Fig.  2.4  as  •j-c.  When  a single  degree  of  freedom 
is  permitted,  the  sum  of  moments  acting  on  the  cage  about  the 
bearing  x axis  is  required  to  be  zero.  This  moment  equation 
considers  the  cage-rolling  element  normal  and  friction  forces 
as  well  as  the  torque  generated  at  the  cage-ring  surface. 

If  there  is  significant  rolling  element  orbital  speed 
variation,  the  cage  is  permitted  to  .move  to  an  eccentric  posi- 
tion with  respect  to  the  land  on  which  it  is  piloted.  Two 
additional  degrees  of  freedom  are  required  to  describe  the 
eccentric  position.  These  additional  degrees  of  freedom  are 
the  cage  center  of  mass  radial  displacement,  e,  and  the  angu- 
lar displacement  of  the  center  of  mass,  with  respect  to  the 
bearing  Y axis,  Oq , See  Figures  (2-5  and  2-6).  These  radial 
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Figure  2-  5 : Inner  Ring-Cage  Land  Contact  Geometry 
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Figure  2-6  : Outer  Ring-Cage  Land  Contact  Geometry 
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friction  forces  as  well  as  the  pressure  buildup  between  the 
cage  and  its  piloting  surface  are  considered  in  the  equilibrium 
equations.  . The  effect  of  the  cage  mass  is  neglected. 

The  ball  bearing  friction  solution  is  thus  obtained  by 
soving  6Z+(1  or  3)  equations  where  Z is  the  number  of  rolling 
elements.  The  ball  bearing  frictionless  solution  is  obtained 
by  solving  1,  (Z/2)  (Z/2+1)  or  Z sets  of  2 equations,  depending 
upon  the  number  of  rolling  elements  in  the  bearing  and  the 
degree  of  load  symmetry  which  prevails.  The  various  symmetry 
conditions  are  explained  below. 

The  roller  bearing  friction  solution  contains  4Z+(1  or  3) 
equations  and  the  frictionless  solution  contains  Z/2,  Z/2+1 
or  Z sets  of  three  equations  again  depending  upon  the  number 
of  rolling  elements  and  whether  or  not  load  symmetry  exists. 

The  various  load  symmetry  conditions  are  as  follows.  Axial 
symmetry  is  utilized  if  the  load  is  axial  only,  then  only  one 
set  of  two  or  three  equations  is  solved  for  the  frictionless 
case  and  six  rolling  elements  and  one  cage  equilibrium  equation 
is  solved  when  friction  is  included.  All  rolling  elements  are 
assumed  to  behave  identically. 

Radial  load  symmetry  is  utilized  if  the  non-axial  shaft 
loading  is'  comprised  of  only  radial  components  parallel  to  the 
Y axis  and  moment  components  parallel  to  the  Z axis  and  the 
position  of  the  first  rolling  element  is  utilized.  When  this 
symmetry  exists,  only  half  the  rolling  elements  need  be  con- 
sidered if  the  number  of  rolling  elements  is  even  and  one 
half  plus  one  need  be  considered  if  the  number  is  odd.  Because 
of  inertia,  terms,  radial  load  symmetry  can  only  be  utilized  in 
the  frictionless  calculations. 

If  load  symmetry  is  not  present,  then  Z sets  of  two  (ball 
bearing)  or  Z sets  of  three  (roller  bearing)  equations  must  be 
solved  to  obtain  the  frictionless  solution., 

As  with  the  steady  state  temperature  mapping  scheme,  the 
Newton- Raphson  scheme  in  subprogram  S0LV13  is  used  to  solve 
the  sets  of  equations  for  each  bearing.  The  iteration  scheme 
terminates  when  either: 
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EPS1  frictionless 
EPS2  friction 

EPS1  frictionless 
EPS2  friction 
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Experience  has  shown  that  the  second  criterion  is  usually 
responsible  for  terminating  the  solution.  However,  when  rolling 
element  loads  are  extremely  large,  on  the  order  of  105  Newtons, 
it  becomes  difficult  to  reduce  the  equation  residues  to  less 
than  10  Newtons.  In  those  instances,  the  first  criteria 
usually  terminates  the  iteration  scheme. 
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3.  PROGRAM  INPUT 


3.1  Types  of  Input  Data 

A complete  set  of  input  data  is  comprised  of  four  distinct 
categories  of  data.  Within  these  categories,  cards  which  con- 
vey specific  kinds  of  information  are  referred  to  as  card 
types.  Depending  on  the  complexity  of  the  problem,  the  input 
data  set  may  contain  none,  one  or  several  cards  of  a given  type. 

The  categories  are  listed  below. 

I.  Title  Cards 

A title  card  plus  a second  card  which  provides  the 
program  control  information  for  the  shaft-bearing 
solution. 

II .  Bearing  Data  Cards 

A set  of  up  to  sixteen  (16)  card  types.  Each  set  describes 
one  bearing  in  the  assembly.  All  bearings  must  be  so 
described.  The  card  sets  must  be  input  sequentially 
in  order  of  increasing  distance  from  a selected  end  of 
the  shaft. 

III..  Thermal  Data  Cards 

A set  of  up  to  nine  (9)  card  types  to  describe  the 
thermal  model  of  the  assembly. 

IV.  Shaft  Data  Cards 

A set  of  three  (3)  card  types  to  describe  the  shaft 
geometry,  bearing  locations  on  the  shaft  and  shaft 
loading. 

If  the  program  is  being  used  to  predict  .the  performance 
of  a bearing  assembly,  cards  from  all  four  sets  must  be  included 
in  the  runstream.  If  the  program  is  being  used  to  thermally 
model  a mechanical  system  wherein  no  bearing  heat  generation 
rates  are  required,  and  therefore  no  bearing  calculations 
need  be  performed,  the  cards  from  sets  II  and  IV  are  omitted. 

The  review  of  required  input  information  which  follows  is 
broken  into  the  four  sets  of  data  categories  given  above,  with 
special  emphasis  on  program  control  data. 

The  input  data  instructions  are  given  in  Appendix  C,  and 
are  for  the  most  part,  self  explanatory.  They  are  laid  out 
in  the  format  of  an  eighty  column  data  card.  A description 
of  the  variables  is  given  in  the  input  instruction  forms. 
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The  units  used  for  input  data  are  as  follows: 

Linear  Dimensions  - (mm) 

Angles  - (degrees) 

Surface  Roughness  (microns) 

Bearing  Angular  Mounting  Errors  - (radians) 

Rotational  Speeds  - (RPM) 

Force  - (Newtons)  (N) 

Moments  - (N-mm) 

Pressure,  Elastic  Modulus  - (N/mm^) 

Density  - (gm/cm^) 

Kinematic  Viscosity  - (cs) 

Temperature  - (degrees  centigrade)  (°C) 

Coefficient  of  Thermal  Expansion  - (°C“1) 

Thermal  Conductivity  - (Watts/m/°C) 

3. 2 Data  Set  L - Title  Cards 

3.2.1  Title  Card  1 . 

This  card  should  contain  the  computer  run  title  and  any 
information  which  might  prove  useful  for  future  identification. 
The  full  eighty  (80)  columns  are  available  for  this  purpose. 

The  title  will  appear  at  the  top  of  each  page  of  Program  output. 

3.2.2  Title  Card  2 


This  card  provides  the  control  information  for  the  shaft 
bearing  solution. 

Item  1;  Shaft  Speed  in  rpm,  GOV  (1) . All  bearings  have 
the  same  shaft,  inner  ring,  and  speed. 

Item  2:  Number  of  Bearings  on  the  shaft  (NBRG) , a minimum 
of  zero  is  permitted  if  no  bearing  solution  is  being  sought.  A 
maximum  of  five  is  permitted.  Note  that  a bearing  is  defined  as 
a single  row  of  rolling  elements.  Thus  a double  row  bearing  is 
treated  as  two  separate,  single  row  bearings. 

Item  3: • Print  Flag  (NPRINT) , > NPRINT  equal  to  zero  is  normal 
and  will  result  in  no  intermediate  or  debug  output.  With  a value 
of  one,  a low  level  intermediate  print  is  obtained  at  the  end  of 
each  shaft  bearing  iteration.  The  values  of  the  variables,  the 
inner  ring  displacements  (DEL) , and  the  equation  residues  are 
printed. 

At  the  end  of  each  bearing  iteration,  wherein  the  rolling 
element  and  cage  equilibrium  equations  are  solved,  an  error 
parameter  is  printed  which  has  the  value: 
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Error  Parameter  = &X^/  ^n-1 

**N  is  the  change  in  the  variable  X specified  at 
iteration  N. 

is  the  value  of  the  variable  specified  at 
the  previous  iteration. 

The  Error  Parameter  is  calculated  for  each  of  the  bearing 
variables,  but  only  the  largest  one  is  printed. 

Additionally,  at  the  end  of  each  Clearance  Change  itera- 
tion, the  clearance  change  error  parameter  is  printed.  This 
error  is  defined: 


DCL^  — DCL^_ 

Error  Parameter  - ^oITIiig  Element  Diameter 


where  DCLjsj  and  DCLj^.j  denote  the  clearance  changes  calculated 
at  the  current  and  previous  iterations  respectively. 

If  NPRINT  is  set  at  2 all  of  the  above  information  is 
printed.  Additionally  the  variable  values  and  residue  values 
are  printed  for  each  iteration  of  the  rolling  element  and  cage 
equilibrium  solution. 

Item  4:  ITFIT  controls  the  number  of  iterations  allowed  to 
satisfy  the  bearing  clearance  change  iteration  scheme.  If 
ITFIT  is  set  to  zero  (0),  or  left  blank,  the  clearance  change 
portion  of  the  program  is  not  executed.  If  a positive  integer 
is  input,  the  clearance  change  scheme  is  utilized  with  a maximum 
iteration  limit  of  five  (5).  If  a negative  integer  is  input, 
the  scheme  is  used  Vrith  a maximum  iteration  limit  equal  to 
the  absolute  value  of  the  negative  integer. 

Item  5:  ITMAIN  limits  the  number  of  iterations  attempted 
during  the  solution  of  the  shaft  and  bearing  inner  ring  equili- 
brium problems-  for  establishing  the  equilibrium  of  bearing 
reactions  and  applied  shaft  loads.  If  ITMAIN  is  left  blank, 
set  to  zero,  or  to  a positive  integer,  then  (15)  iterations 
are  permitted.  If  ITMAIN  is  set  to  a negative  integer  the 
number  of  iterations  is  limited  to  the  absolute  value  of  that 
integer. 
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Item  6 : G0V(2)  or  EPSFIT  is  the  convergence  criterion 
for  the  diametral  clearance  change  portion  of  the  analysis. 

As  mentioned  under  item  3 above,  this  error  parameter  is 
defined  by  Eq.  2.28. 

The  iteration  scheme  is  terminated  when  the  error  parameter 
is  less  than  the  input  value  of  EPSFIT,  If  EPSFIT  is  left 
blank  or  is  set  to  zero  (0),  the  program  default  value  of 
0.0001  is  used. 

Items  7 § 8:  Main  loop  accuracy  for  frictionless  elastic 
(EPSl^  and  friction  solution  (EPS2).  These  accuracy  values 
control  the  accuracy  of  the  shaft  bearing  deflection  solution 
as  well  as  the  quasi-dynamic  solution  of  the  component 
dynamics.  If  EPS1  and  EPS 2 are  left  blank  or  set  to 

zero  CO) , default  values  of  0.001  and  0.0001  respectively 
are  used. 

Item  9:  IMT,  if  set  to  1,  tR.e  material  properties  for 
both  bearing  rings  and  the  rolling  elements  are  to  be  input 
on  card  types  B 11  through  B 19.  If  IMT  is  zero  or  blank,  the 
rings  and  rolling  elements  are  assumed  to  be  52100  bearing  steel. 
Card  types  B 11  through  B 14  are  required  only  if  the  change 
in  bearing  diametral  clearance  is  to  be  calculated. 

Item  10:  NPASS  controls  the  level  of  the  bearing  solution: 

0 Elastic  Contact  Forces  are  calculated.  No'  lubrication 
or  friction  effects  are  considered. 

1 Elastic  Contact  Forces  are  calculated.  Lubrication 
and  friction  effects  are  considered  using  raceway 
control  (ball  bearing)  or  epicyclic  (roller  bearing) 
assumptions  to  estimate  rolling  element  and  cage  speeds. 

2 Inner  Equilibrium  is  satisfied  considering  only  the 
Elastic  Contact  Forces.  Using  the  inner  ring  positions 
thus  obtained,  rolling  element  and  cage  equilibrium 
are  determined  considering  friction. 

3 Complete  Solution.  The  inner  ring,  rolling  element 
and  cage  equilibrium  is  determined  considering  all 
elastic  and  friction  forces. 

3. 3 Data  Set  II  - Bearing  Data 

Most  of  the  input  instructions  are  self-explanatory. 

Where  certain  items  are  deemed  to  require  more  explanation 
than  given  in  the  input  data  format  instructions  they  are 
treated  on  an  individual  basis  by  card  type  and  item  number. 
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Most  of  the  bearing  input  data  is  read  into  a two  dimen- 
sional array  named  "BD",  which  has  the  dimensions  (1830,  5). 

For  each  of  the  five  bearings  permitted  on  a shaft,  a total 
of  1830  pieces  of  data  may  be  stored.  Denoting  BD(I,J),.  I 
represents  a specific  piece  of  bearing  data,  J represents 
the  bearing  number.  The  bearing  input  data  of  Data  Set  II 
occupies  the  first  106  locations  of  the  1830  allotted.  On 
the  input  data  format  sheets  the  designation  BD(I)  where 
1=1... 106,  denotes  the  location  within  the  BD  array  where 
each  piece  of  input  data  is  stored. 

3.3.1  Card  Type  1 - Bearing  Type  and  Material  Designations 

Item  1 : Bearing  type,  columns  1-10  must  be  specified, 
left  justified,  i.e.,  "B",  "C"  or  "T"  in  column  1.  This  format 
must  be  followed  since  the  Program  recognition  of  bearing  type, 
(ball,  cylindrical  or  tapered  roller  bearing),  is  derived  from 
reading  the  "B" , "C",  or  "T"  in  the  first  column  of  this  card. 

Items  283:  Columns  11^30  and  31-50,  "Steel  designations", 
inner  and  outer  rings  respectively.  The  alphameric-literal 
description  of  the  steel  types  such  as  "M-50"  or  "AISI  52100" 
is  input. 

Items  485:  Columns  51-60  and  61-70,  the  numbers  input 
for  items  4 and- 5 are  used  to  account  for  improved  materials 
and  multiply  the  raceway  fatigue  lives  as  determined  by 
Lundberg-Palmgren  methods.  Typical  life  factor  values  for  modern 
steels  are  in  the  neighborhood  of  2.0  to  3.0.  If  the  ASME 
Publication  Life  Adjustment  Factors  for  Ball  and  Roller  Bearings, 
is  referenced  by  the  user , the  ^Material  Factor  D and  the  ~~ 
Material  Process  Factor  E should  be  used  as  multipliers 
as  inputs  for  items  4 and  5.  Additionally  if  the  user  is 
accustomed'  to  using  a lubricant  life  multiplier  he  must  also 
multiply  the  material  factor  by  the  lubricant  life  multiplier. 

The  program  considers  EHD  film  thickness  and  RMS  surface  rough- 
ness but  generates  a life  multiplier  having  a maximum  value  of 
1 and  a miminum  of  0.479,  i.e.  Lube-Life  Factor  Programmed  only 
serves  to  reduce  predicted  Fatigue  Life; 

Item  6:  Columns. 71-78,  "Orientation  angle  of  the  first 
rolling  element".  (£1)  (degrees).  Refer  to  Fig.  (2.4).  The 
quasi-dynamic  rolling  element  bearing  problem  has  an  infinite 
number  of  solutions  which  fall  within  a narrow  envelope  having 
a periodic . shape.  The  solution  obtained  is  a function  of 
the  rolling  element  positions  relative  to  the  bearing  system 
coordinate  axes.  4 1 s places  a rolling  element  on  the  Y 
axis  and  is  the  choice  customarily  made,  ^.can  be  desig- 
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nated  as  any  value  0<4><  where  Z is  the  number  of  foiling 
elements.  For  each  cTiTferent  value  assigned  to  $1  a differ- 
ent, although  similar,  bearing  solution  will  be  obtained.  To 
take  advantage  of  bearing  symmetry  and  the  computer  time  savings 
which  result  must  be  specified  as  zero  or  left  blank. 

Item:  7 Column  80,  a signal,  termed  the  crown  drop  flag, 

which  specifies  for  a cylindrical  or  tapered  roller  bearing, 
whether  the  roller-race  crown  drops  will  be  calculated,  or  read 
directly.  If  item  7 is  blank  or  zero,  the  crown  drops  are  cal- 
culated based  on  the  roller-race  crown  radius,  and  effective  flat 
length  input  information.  If  the  crown  drop  flag  is  other  than 
zero  or  blank  the  non-uniform  separation  of  the  roller  and  race- 
way must  be  specified  at  the  center  of  each  slice  into  which  the 
roller-raceway  effective  contact  length  is  divided.  The' slice 
widths  are  identical.  The  number  of  slices  is  input  as  item  7 
card  type  B4.  The  non-uniform  roller-raceway  separation  is  in- 
put on  card  types  B5  and  B6. 

3.3.2  Card  Type  B2  - Bearing  Geometry  and  Outer  Ring  Speed 
3. 3. 2.1.  Ball  Bearing  Geometry 

Items  1,  2 and  6 are  self-explanatory.  Item  5 pertains 
only  to  a tapered  roller  bearing,  as  discussed  later.  Items  3 
and  4 require  explanation  however. 

Through  the  proper  specification  of  the  diametral  clearance 
and  contact . angle,  the  program  can  properly  handle  deep  groove, 
split  inner,  and  angular  contact  ball  bearings. 

The  deep  groove  ball  bearing  requires  the  specification 
of  zero  contact  angle  and  either  the  operating  diametral  clear- 
ance Pd  or  the  off  the  shelf  diametral  clearance,  if  the  dimen- 
sional change  analysis  is  utilized. 

The  angular  contact  bearing  is  fully  described  through 
specification  of  the  contact  angle  measured. under 
axial  load.  However,  this  method  of  input  does  not  accurately 
define  the  system  if  there  is  more  than  one  angular  contact 
supporting  the  shaft  and  at  least  one  of  those  bearings  has 
its  grooves  offset  in  the  direction  opposite  to  the  other 
bearings  and  if  the  shaft  is  capable  of  axial  and/or  radial 
play.  In  other  words,  if  what  are  known  as  angular  contact 
ball  bearings  are  mounted  such  that  some  diametral  shaft 
play  is  permitted,  an  auxilliary  angle . as  well  as  the  diametral 
play  must  be  specified  at  input.  The  angle  input  is  not  the 
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manufacturer's  designated  contact  angle, ©< , but  an  auxilliary 
angle, o<  0,  the  calculation  for  which  shall  be  demonstrated. 

Refer  to  Figure  3.1.  The  manufacturer's  contact  angle  is 
calculated  as  follows: 

* -iTaA-Pil 

o CoS  L ZA  J (3.1) 

^ - rD  + v\  *-D  (3.2) 


where:  r0  and  r^  are  the  outer  and  inner  raceway  groove 

radii  respectively 

D is  the  ball  diameter 

Under  a gauge  axial  load,oCis  obtained  at  both  inner  and 
outer  raceways  for  each  ball.  Under  this  condition,  the  outer 
and  inner  raceways  are  axially  offset  an  amount^. 


(3.3) 


When  angular  contact  ball  bearings  are  mounted  with  some 
diametral  play,  the  grooves  are  offset  an  amount  Sa0  such  that 
S®(0<S«(.  The  diametral  play  which  obtains  at  this  condition 
is  Sd.  • This  diametral  play  is  usually  known  by  the  engineer 
or  designer  and  is  usually  required  to  allow  some  forgiveness 
when  thermal  gradients  are  encountered.  Assuming  that  the  user 
has  the  values  for  «<  , r0,  r^,  D and  S*(0  then: 


where:  Pd  and  A may  be  calculated  from  Eqs.  (3.1)  and  (3.2). 


If  Sa0  i-s  unknown,  the  following  equation  may  be  solved 
for  a0 • 
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, FIC.  3.1  ANGULAR  CONTACT  BALL  BEARING  GEOMETRY 


FIG.  3.1  ANGULAR  CONTACT  BALL  BEARING  GEOMETRY 


Geometry  Showing  the  as  Mounted 
Groove  Alignment  and  the  Input 
Diametral  Play  (Sd) . 


Geometry  Showing  Groove 
Alignment  with  Rings 
Subjected  to  a Gage  Radial 

Load,  i.e.,  with  the  Radial 
Play  Removed.  Also  shown, the 
Auxilliary  Input  Contact 
Angle  o<0 


(3.5) 


z a 

•z.  A 


-1 


In  order  that  the  Program  properly  handle  split  inner 
ring  ball  bearings  an  auxilliary  angle  and  diametral  play  must 
be  input.  Referring  to  Figure  3.2,  the  auxilliary  angle 0 and 
diametral  play  Sd  must  be  determined  and  input.  Typically 
the  values  of  D,  r0,  r^,^(s  and  Sd’  are  known  and  Pd  may  be 
calculated  as  follows: 

YA-  + (-2.rC--D)(  l-  Costs')  (3.6) 


The  unloaded  half  of  the  inner  ring  must  be  removed  from  con- 
sideration and  the  ball  moved  such  that  its  center  lies  on  the 
line  connecting  the  origins  of  r^  and  r0  and  positioned  such 
that  the  auxilliary  clearance  Sd/4  exists  at  both  the  inner 
and  outer  raceways.  The  auxilliary  angle  is 


-t  ^ /i) 


(3.7) 


The  angle  associated  with  each  ball  bearing  must  be 
specified  with  the  correct  sign.  A positive  contact  angle 
, allows  the  bearing  to  accept  a positively  directed  axial  load 
* transmitted  by  the  shaft. 
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WITH  SHIM  IN  PLACE 


(a) 


3. 3. 2. 2 Tapered  Roller  Geometry 

Items  2 and  6 are  self-explanatory. 

Item  1;  Bearing  Pitch  Diameter 

The  tapered  roller  bearing  pitch  diameter  may  be  calculated 
through  specification  of  the  roller  measured  large  end  diameter, 
the  roller  large  end  corner  radius,  one  half  the  cup  included 
angle,  the  roller  included  angle  and  the  roller  total  length. 

This  calculation  is  performed  within  the  program.  The  user  there 
fore  need  not  specify  the  tapered  roller  bearing  pitch  diameter. 

Item  3:  Bearing  Axial  Play 

For  a tapered  roller  bearing,  the  bearing  axial  play  rather 
than  diametral  clearance  must  be  specified.  (See  Fig.  3.3)  Note 
that  this  end  play  pertains  to  only  the  bearing  in  question.  For 
two  identical  tapered  roller  bearings  on  a shaft,  one  half  the 
total  shaft  axial  play  should  be  specified  for  each  bearing.  For 
two  tapered  roller  bearings  of  dissimilar  size,  the  total  shaft 
axial  play  should  be  apportioned  according  to  the  bearing  size 
such  that  the  sum  of  the  axial  plays  specified  for  the  two 
bearings  equals  the  total  for  the  shaft. 

Item  4:  Bearing  Contact  Angle 


For  the  tapered  roller  bearing,  one -half  the  included  CUP 
angle  (a  ) is  input  as  the  contact  angle.  (See  Fig,  3.3)  inis 
angle  mu$t  be  specified  with  the  correct  sign.  A positive  . 
angle  allows  the  bearing  to  accept  a positively  directed  axial 
load  transmitted  by  the  shaft  and  vice  versa  for  a negative  angle. 


Item  5:  Tapered  Roller  Bearing  Flange  Angle 

The  flange  angle  is  shown  by  a f in  Fig.  3.3  The  flange 
angle  must  always  be  positive. 


3. 3. 2. 3 Cylindrical  Roller  Bearing  Geometry 

Items  1,  2,  3 and  6 are  self  explanatory.  Both  items  4 
and  5 should  be  left  blank. 
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FIG.  3.3  TAPERED  ROLLER  BEARING  GEOMETRY 


-p* 

o 


Note:  Circled  nomenclature  represent 
input  variable. 


3.3.3  Card  Type  B3  - Rolling  Element  Geometry 


3. 3. 3.1  Ball  Geometry 

The  geometry  of  a ball  is  fully  defined  by  its  diameter. 

3. 3. 3. 2 Tapered  Roller  Geometry 

Item  1 - Roller,  Measured.  Large  End  Diameter 

For  a tapered  roller,  the  input  diameter  required  is  the 
largest  measureable  diameter.  (See  Fig.  3.4.)  Typically 
this  measurement  should  be  taken  where  the  large  end  corner 
radius  becomes  tangent  to  the  roller  surface  profile.  Within 
the  program  a "working"  large  end  diameter  is  calculated. 

This  diameter  is  shown  as  D„  in  Fig.  3.4.  All  bearing  geometrical 
relationships  are  calculated  based  on  Dw. 


(3.8) 


where:  D is  the  measured  large  end  diameter  fe.o  is 

the  distance  from  the  roller  end  to  the 
beginning  of  the  roller  effective  length  at 
the  outer  raceway  surface,  measured  parallel 
to  the  roller  surface.  If  Vo  is  the  corner 
radius  at  the  roller  large  end 

r-Q  Q V- 

m do5V  (3.9) 

IS  is  one  half  the  roller  included  angle. 


Item  2 through  6 are  shown  in  Fig.  3.4  and  are  self 
explanatory.  Note  that  the  program  can  handle  a nonzero  roller 
flat  length.  Item  6,  Most  tapered  rollers  are,  however,  fully 
crowned. 
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3. 3. 3. 3 Cylindrical  Roller  Geometry 

Typically,  cylindrical  rollers  are  partially  crowned  as 
shown  in  Fig.  3.5.  The  center  of  the  roller  is  flat.  Toward 
the  ends  the  roller  profile  is  formed  by  a crown  radius.  There 
are  usually  rounded  corners  at  the  roller  ends.  These  corners 
reduce  the  load  carrying  surface  of  the  roller  such  that  if 
there  are  no  raceway  undercuts  the  roller  raceway  effective 
length  equals  the  roller  total  length  less  the  two  corner  radii 
(see  Section  3. 3. 4. 2).  Note  that  a partially  crowned  roller 
is  specified  through  input  of  a non- zero  flat  length.  If  the 
flat  length  is  zero,  the  roller  is  fully  crowned  with  its  pro- 
file defined  by  the  roller  crown  radius. 

The  roller  end  sphere  radius  is  defined  in  Fig.  3.5.  For 
a cylindrical  roller,  the  roller  included  angle  is  zero. 


3.3.4  Card  Type  B4  - Rolling  Element-Ring  Geometry 

3. 3. 4.1  Ball  Bearing 

Items  1 and  2 refer  to  the  outer  and  inner  raceway  curva- 
ture s-respectTveTy  where  curvature  is  defined  as  the  cross 
groove  radius  divided  by  the  ball  diameter.  Typical  values 
range  from  0.515  to  0.57. 

3. 3.4. 2 Tapered  Roller  Bearing  Contact  Geometry 

Items  1 and  2 - Roller-Raceway  Effective  Length 

The  roller-raceway  load  bearing  surface  is  measured  para- 
llel to  the  roller  surface  such  that  if  there  were  no  relief  at 
the  roller  ends  the  effective  contact  length  would  be: 

iTe.  ■=  L (3.10) 

However,  since  the  roller  has  corners  at  the  large  and  small 
ends,  the  actual  effective  length  is  less  than  Ue  . The  con- 
sideration of  raceway  undercuts  at  the  inner  raceway  flanges 
may  result  in  a raceway  effective  length  less  than  the  roller 
effective  length  in  which  case  the  shorter  of  the  two  should  be 
input. 
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«(£  Flange  Angles 
Rq;  Crovm  Radii. 

Le  Roller-Raceway  Effective  Lengths 
Lf  Roller  Flat  Length 
RS  Roller  End  Sphere  Radius 
Ep  Raceway  End  Play 

FIG.  3.5  CYLINDRICAL  ROLLER  BEARING  GEOMETRY 
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Items  3 and  4 - Raceway  Crown  Radius 


The  present  analysis  permits  both  the  roller  and  race- 
way to  be  crowned.  If  the  raceway  is  crowned,  it  must  be  fully 
crowned  with  no  flat  specified.  If  the  raceways  are  flat  the 
input  crown  radius  may  be  left  blank,  in  which  case  a default 
value  of  1.  x 10+  10  inches  is  used.' 

Note  that  the  unloaded  roller-raceway  separation  along 
the  roller  profile,  (6c)  Fig.  3.6,  calculated  at  the  center  of 
each  roller  raceway  slice  is  comprised  of  the  sum  of  the  roller 
and  raceway  crown  drops . 


Items  5 and  6 - Roller  Large  End  Corner  Relief 

These  data  specify  the  distance  from  the  roller  large 
end  to  the  point  on  the  roller  surface  where  the  roller  effective 
length  begins.  For  the  outer  raceway  contact  this  distance 
may  be  calculated  using  Eq . 3.9.  For  the  inner  raceway  use  Eq,3.9 
or  the  width  of  the  inner  raceway  undercut  at  the  large  end. 

Item  7 - The  number  of  slices  into  which  the  roller 
raceway  contacts  are  divided. 

A maximum  value  of  twenty  (20)  is  permitted  a default 
value  of  eleven  (11)  is  used  if  a blank  or  zero  is  read. 

3. 3. 4. 3.  Cylindrical  Roller  Contact  Geometry 

3. 3. 4. 3.1  Card  Type  B4A  Roler- Raceway  Geometry 

Items  1 through  4 and  7 have  the  same  definitions  as 
they  had  for  the  tapered  roller  bearing. 

Items  5 and  6 - Roller  End  Corner  Relief  is  not  required 
input  for  the  cylindrical  roller  bearing  since  the  roller 
raceway  effective  length  is  assumed  to  be  centered  along  the 
roller.  This  is  not  the  case  for  the  tapered  roller  bearing. 

3. 3. 4. 3. 2 Card  Type  B4B  Roller  Flange  Geometry  for 

Cylindrical  Roller  Bearings 
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Figure  3.6  Roller- Race  Lamination  Showing  Relative  Approach 
& * (5u)  and  Crown  Drop  (6c) 


Items  1 - 4 - The  flange  angle  and  end  play  definitions 
can  be  seen  in  hig.  3.5.  Note  that  when  a ring  has  a single 
flange  such  as  Inversion  No.  3.  Fig.  3.7,  the  end  play  is  the 
distance  between  the  roller  and  the  flange  when  the  roller  is 
centered  on  the  raceway. 

Item  5 - The  Flange  Inversion  Index  - See  Fig.  3.7 
The  number  which  corresponds  to  tne  particular  Flange  inversion 
being  examined  must  be  input.  Note  that  since  the  inversions 
greater  than  eight  (8)  cannot  carry  axial  loading,  the  bearing 
carries  load  only  on  the  raceways,  and  thus  the  Program  resets 
the  inversion  index  to  one  if  the  input  value  is  greater  than 
eight.  The  inversion  index  must  be  input  as  a real  number,  (with 
a decimal  point) . r 

3.3.5  Roller-Raceway  Non-Uniform  Profile  Definition 
3. 3. 5.1  Card  Type  B5  - Inner  Roller  Raceway  Contact 

These  cards  are  used  to  input  the  inner  and  outer  race 
roller-race  separation  along  the  roller  profile.  With  the 
high  points  of  the  roller  and  race  in  contact,  i.e.,  with  all 
clearance  between  roller  and  raceway  removed.  These  cards 
must  be  omitted  if  item  7 of  the  Bearing  Data  Title  card  is 
zero  or  blank.  These  data  are  used  only  when  the  roller-raceway 
profile  geometry  cannot  be  defined  by  card  types  B3  and  B4A. 

3.3.6  Card  Type  B6  - Outer  Roller  Raceway  Contact 

Same  as  Card  Type  B5. 


3.3.7  Ring-Rolling  Element  Surface  Data 

3. 3. 7.1  Card  Type  B7A  - Raceway  - Rolling  Element  Surface'  Data 

Items  1 through  6 define  the  statistical  surface  micro- 
geometry parameters  of  the  rollers  and  raceways.  Items  1 
through  3 require  the  input  of  center  line  average  CLA  surface 
roughness.  Within  the  program  CLA  values  are  converted  to  RMS 
by  multiplying  by  0.9. 

Items  4 through  6 are  RMS  values  of  the  slopes  measured 
in  degrees,  of  the  surface  asperities  as  measured  in  a traverse 
across  the  groove  for  rings,  longitudinally  for  rollers  and 
in  any  arbitrary  direction  for  balls.  Typical  values  for  race- 
way and  rolling  element  surfaces  are  1 to  2 degrees.  This  card 
is  omitted  if  the  solution  level  is  NPASS  = 0. 


FIG.  3.7  CYLINDRICAL 
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3. 3. 7. 2 Card  Type  B7B  - Flange-Roller  End  Surface  Data 

The  data  are  identical  to  the  data  in  Card  Type  B7A  but 
refer  to  the  flange  and  roller  end  surfaces  rather  than  the 
raceway  and  roller  rolling  surfaces. 

Note  that  both  card  types  B7A  and  B7B  are  omitted  if  the 
solution  level  is  NPASS  =0. 


3.3.8  Card  Type  B8  - Cage  Data 

This  card  is  omitted  if  the  solution  level  is  NPASS  =0. 
These  data  are  self  explanatory.  Note  that  the  cage  weight  is 
an  input  item.  The  weight  however,  is  not  used  in  any  calcu- 
lation. It  is  included  only  for  future  consideration  of  cage 
stability  predictions. 


3.3.9  Card  Type  B9  - Shaft  and  Housing  Fit  Dimensions 

These  cards  are  to  be  included  only  if  the  change  in 
bearing  diametral  clearance  with  operating  conditions  is  to  be 
calculated,  i.e.  if  item  • 4,'ITFIT,  on  the  Bearing  Title  Card  is 
non-zero.  On  Card  Type  B9,  tight  interference  fits  bear  a 
positive  sign  and  loose  fits,  a negative  sign. 

Item  3 and  6 on  Card  No.  9 are  termed  the  shaft  and 
housing  effective  widths  respectively.  The  value  specified 
for  these  effective  widths  may  be  as  great  as  twice  the  ring 
width . 

Use  of  an  effective  width  is  an  attempt  to  account  for  the 
greater  radial  rigidity  of  a shaft  being  longer  than  the  ring 
that  is  pressed  on  to  it,  owing  to  the  fact  that  the  shaft 
deflects  over  a distance  that  extends  beyond  the  ring  width. 

In  the  program  the  calculated  internal  pressure  on  the  ring 
due  to  its  interference  fit  with  the  shaft,  is  distributed 
over  the  shaft  effective  width  and  this  (lower)  pressure  is 
used  in  computing  the  shaft  deflection.  Using  double  the 
actual  width  as  the  effective  width  is  customary. 
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3.3.10  Card  Type  BIO  - Shaft  Housing  Fit  Dimensions 
These  items  are  self  explanatory. 

Note:  Bearing  System  Components  Material  Properties.  Card 

types  Bll  through  B14. 

These  card  types  define  the  material  properties  of  the 
shaft,  inner  ring,  rolling  element,  outer  ring  and 
housing,  data  items  1 through  5 respectively.  This 
set  of  cards  is  to  be  included  if  either  the  bearing 
clearance  change  analysis  is  used,  i.e.  item  4,  ITFIT, 
Bearing  Title  Card  2 is  non  zero,  or  if  the  bearing 
rings  or  rolling  elements  are  not  steel,  i.e.  item  9, 
IMT,  Bearing  Title  Card  2 is  equal  to  1.  If  any  item 
on  card  types  Bll  through  B14  is  left  blank,  the 
program  inserts  the  appropriate  value  of  the  steel 
property. 


3.3.11  Card  Type  Bll  - Elastic-Moduli 

This  card  defines  the  elastic  modulus  for  the  shaft,  inner 
ring,  rolling  element,  outer  ring,  and  housing  respectively. 


3.3.12  Card  Type  B12  - Poisson's  Ratio 

This  card  defines  the  Poisson's  ratio  for  the  shaft,  inner 
ring,  rolling  element,  outer  ring,  and  housing  respectively. 


3.3.13  Card  Type  B13  - Density 

This  card  defines  the  density  for  the  shaft,  inner  ring, 
rolling  element,  outer  ring,  and  housing  respectively. 


3.3.14  Card  Type  B14  - Coefficient  of  Thermal  Expansion 

This  card  defines  the  coefficient  of  thermal  expansion  for 
the  shaft,  inner  ring,  rolling  element,  outer  ring,  and 
housing  respectively. 


3.3.15  Card  Type  B15  - Lubrication  and  Friction  Data 

This  card  is  omitted  if  the  solution  levels  NPASS  = 0. 
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Items  1 and  2 


Items  1 and  2 are  the  amounts  by  which  the  combined  thick- 
ness of  the  lubricant  film  on  the  rolling  track  and  rolling 
element  is  increased  during  the  time  interval  between  the  pas- 
sage of  successive  rolling  elements,  from  whatever  replenish- 
ment mechanisms  are  operative.  Item  1 applies  to  the  outer  and 
Item  2 to  the  inner  race-rolling  contacts  respectively.  If 
Item  1 is  zero  or  blank  the  mode  of  friction  is  assumed  to  be 
dry. 

At  the  present  time  the  magnitude  of  the  inner  and  outer 
raceway  replenishment  layers  has  not  been  correlated  to  lubri- 
cant flow  rate,  lubricant  application  methods  and  bearing  size 
and  speed  factors.  At  this  point  then,  the  user  is  forced  to 
establish  proper  values  for  the  replenishment  layer  thickness. 

As  a rough  guide  the  following  suggestions  are  made. 

1)  To  avoid  starvation,  the  replenishment  layer  thick- 
nesses should  be  one  to  two  times  the  EHD  film  thickness 
which  develops-  in  the  rolling  element  raceway  contacts. 

2)  Because  of  centrifugal  force  effects,  intuition 
suggests  that  the  outer  raceway  replenishment  layer 
should  be  several  times  thicker  than  that  prescribed 
at  the  inner  raceway. 

Item  3,  XCAV,  describes  the  percentage  of  the  bearing 
cavity,  estimated  by  the  user  to  be  occupied  by  the  lubricant. 
0<XCAV<100 . 

As  with  the  replenishment  layer  thicknesses,  the  amount 
of  free  lubricant  should  be  able  to  be  correlated  with  lubri- 
cant flow  rate,  lubricant  application  methods  and  bearing  size 
and  speed  factors.  At  this  time  such  correlations  do  not  exist. 
XCAV  values  less  than  five  percent  are  recommended  at  this  point. 

Item  4 is  the  coefficient  of  coulomb  friction  applicable 
for  the  contact  of  asperities.  If  Item  1 and  2 are  zero,  then 
Item  4 serves  as  the  coulomb  friction  coefficient  which  prevails 
in  all  contacts. 

Items  5 and  6 are  the  lubricant  replenishment  layer  thick- 
nesses for  the  outer  and  inner  ring  flanges  respectively. 

Item  7 is  the  coefficient  of  coulomb  friction  applicable 
for  the  asperity  interactions  at  the  roller  end- flange  contacts. 

3.3.16  Card  Type  B16 

This  card  is  omitted  if  Item  9 title  card  2 is  zero  or 
blank  or  if  Item  1.  card  15  is  zero  or  blank  which  implies  dry 
friction. 
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This  card  specifies  the  lubricant  type.  If  Item  1,  NCODE 
is  1,  2,  3,  or  4 the  Program  uses  preprogrammed  lubricant 
properties  as  presented  in  Table  1,  and  no  further  information 
is  required. 

NCODE  Lubricant 

1 A specific  mineral  oil 

2 A MIL-L- 7808G 

3 Polyphenyl-Ether 

4 A MIL-L- 23699 

NCODE  may  also  be  specified  as  negative  (-1  to  -4),  in 
which  case  the  traction  characteristics  of  the  respective  lubri 
cant  NCODE  noted- above  are  used  but  the  actual  properties 
specified  by  Items  2 through  7 override  the  hard  coded  data. 
This  option  is  most  useful  in  specifying  various  mineral  oils 
i.e.  NCODE  - -1. 

3.4  DATA .SET  III  - THERMAL  MODEL  DATA 

Appendix  B has  been  included  to  aid  the  user  in  data 
preparation  and  calculation  of  heat  transfer  coefficients 
required  at  input. 

3.4.1  Card  Type  T1 

Card  type  1 is  a control  card.  If  no  temperature  map 
is  to  be  calculated,  this  card  is  to  be  included  as  a blank 
card  followed  by  a Type  2 card  for  each  bearing  on  the  shaft. 

Card  Type  1 containes  control  input  for  both  steady  state  and 
transient  thermal  analyses.  It  is  not  intended  however,  that 
both  analyses  be  executed  with  the  same  run. 

Item  1:  The  highest  node  number  (M) . The  temperature 
nodes  must  be  numbered  consecutively  from  one  (1)  to  the 
highest  node  number.  The  highest  node  number  must  not  exceed 
one.  hundred  (100)  . 

Item  2:  Node  Number  of  the  Highest  Unknown  Temperature 
Node  (N) . This  number  should  equal  the  total  number  of 
unknown  node  temperatures.  It  is  required  that  all  nodes 
with  unknown  temperatures  be  assigned  the  lowest  node  numbers. 
The  nodes  which  have  known  temperatures  are  assigned  the 
highest  numbers. 

Item  3:  Common  Initial  Temperature  (TEMP)°C:  The  tempera 
ture  solution  iteration  scheme  requires  a starting  point,  i.e., 
guesses  of  the  equilibrium  tempratures.  Card  Type  3 allows 
the  user  to  input  estimates  of  individual  node  temperatures. 
When  a node  is  not  given  a specific  initial  temperature,  the 
temperature  specified  as  Item  3 of  Card  Type  1 is  assigned. 

Item  4:  Punch  Flag  (IPUNCH):  If  the  Punch  Flag  is  not 
zero  (0)  or  Blank,  the  system  equilibrium  temperatures  along 
with  the  respective  node  numbers  will  be  punched  according 
to  the  format  of  Card  T3.  This  option  is  useful  if,  for 
instance,  the  user  makes  a steady  state  run  with  lubrication, 
and  then  wishes  to  use  the  resultant  temperature  as  the  initia- 
tion point  for  a transient  dry  friction  run  in  order  to  assess 
the  consequence  of  lubricant  flow  termination. 
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TABLE  1 -A 

LUBRICANT  PROPERTIES  OF  FOUR  OILS  USED 


Oil  Oil  Kinematic  Walther  Equation  Density  Thermal  Thermal  Temp. 

No.  Type  Viscosity,  (cs)  Constants  @ 60°F  Conductivity  Coeff.  Viscosity 

100°F  210°F  A B gm/cm3  Kf  of  Coeff. 

. P Btu/Hr/Ft/°R  Expansion  ft  (°R~1) 

G 


1 

Mineral 

Oil 

64 

00 

• 

o 

10.349 

3.673 

0.8800 

0.0671 

3.52 

X 

io-4 

0.0193 

2 

MIL-L-7808G 

12.8 

3.2 

10.215 

3.698 

0.9526 

0.0879 

3.94 

X 

10‘4 

0.0132 

3 

C-Ether 

25.4 

4.13 

11.452 

4.113 

1.201 

0.0690 

4.15 

X 

io"4 

0.0168 

4 

MIL-L- 23699 

28.0 

5.1  • 

10.207 

3.655 

1.010 

0.0879 

4.14 

X 

10“4 

0.0161 

TABLE  X-B 


TABULATION  OF  CONSTANT  FOR  FOUR  OILS 


Oil  No. 

Oil  Type 

Aj 

*2 

Pl  x 105 

2l_ 

c2_ 

(n  r)  B 

(Vk)B 

1 

Mineral  (1) 

0.3 

3.42 

2 . 14 

1.5 

1.8 

115.9 

0.65 

0.10 

2 

Oil 

MIL-L-7808 

(2) 

0.94 

4.08 

1.48 

1.7 

17.6 

39.0 

0.68 

0.25 

v-n  3 

C-Ether 

0.70 

3.29 

1.50 

1.7 

26.2 

6.6 

0.65 

0.15 

•p- 

4 

MIL-L-23699 

(2) 

0.93 

6.88 

3.44 

2.2 

10.4 

47.3 

0.68 

0.25 

NOTES : 

(1)  obtained  from  Traction  Data  reported  by  Johnson  and  Cameron  (28) 

(2)  S K F tests 

(3)  obtained  from  Traction  Data  reported  by  Smith  et  al  (29) 


Item  5:  "Output  Flag"  (IUB).  If  the  "Output  Flag"  is 

not  zero,  the  bearing  program  output  and  a temperature  map  will 
be  printed  after  each  call  to  the  shaft  bearing  solution  scheme. 
This  printout  will  allow  the  user  to  observe  the  flow  of  the 
solution  and  to  note  the  interactive  effects  of  system  temper- 
atures and  bearing  heat  generation  rates.  Since  the  temperature 
solution  is  not  mathematically  coupled  to  the  bearing  solution 
the  possibility  exists  that  the  solution  may  diverge  or  oscil- 
late. In  such  a case,  study  of  the  intermediate  output  produced 
by  the  "Output  Flag"  option  may  provide  the  user  with  better 
initial  temperature  estimates  that  will  effect  a steady  state 
solution.  Two  levels  of  bearing  output  are  permitted.  If  IUB 
is  1,  the  rolling  element  output  is  not  printed.  If  IUB  is  2, 
full  bearing  output  is  obtained. 

Item  6:  "Maximum  Number  of  Calls  to  the  Shaft  Bearing 
Program"  (IT1) . IT1  is  the  limit  on  the  number  of  Thermal- 
Shaft-Bearing  iterations,  i.e.,  the  external  temperature  equi- 
librium calculation.  The  user  must  input  a non-zero  integer 
such  as  5 or  10  in  order  for  the  program  to  iterate  to  an 
equilibrium  condition.  If  IT1  is  left  blank  or  set  to  zero  (0) 
or  1,  shaft  bearing  performance  will  be  based  on  the  initially 
estimated  temperatures  of  the  system.  The  temperatures  printed 
out  will  be  based  on  the  bearing  generated  heats.  It  is 
unlikely  that  an  acceptable  equilibrium  condition  will  be 
achieved.  However,  the  temperatures  which  result  may  provide 
better  initial  estimates  for  a subsequent  run  than  those 
specified  by  the  user. 

IT1  also  serves  as  a limit  on  the  transient  temperature 
solution  scheme,  by  limiting  the  number  of  times  the  shaft- 
bearing solution  scheme  is  called.  Each  call  to  the  shaft- 
bearing scheme  will  input  a new  set  of  bearing  heats  to  the 
transient • temperature  scheme  until  a steady  state  condition 
is  approached  or  until  the  transient  solution  time  limit 
is  reached. 

Item  7:  "Absolute  Accuracy  of  Temperatures  for  the 
External"  TKermal  Solution"  (EPI).  In  the  steady  state  thermal 
solution  scheme,  each  calculation  of  system  temperatures  occurs 
after  a call  to  the  shaft-bearing  scheme  which  produces  bearing 
generated  heats.  After  the  system  temperatures  have  been 
calculated  for  each  iteration,  using  the  internal  temperature 
solution  scheme,  each  node  temperature  is  checked  against  the 
nodal  temperature  at  the  previous  iteration. 

If  £t(N)i”  t(N-l)i^<  EP1  £or  a11  nodes  1 then  equilibrium 
has  been  achieved  and  the  iteration  process  stops. 
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I tern  8 : "Iteration  Limit  for  the  Internal  Thermal  Solution" 

( I T 2 ) “ ATter  each  call  to  the  shaft  bearing  program,  the 
internal  temperature  calculation  scheme  is  used  to  determine 
the  steady  state  equilibrium  temperatures  based  on  the  calculated 
set  of  bearing  heat  generation  rates.  If  the  program  is  used 
to  calculate  the  temperature  distribution  of  a non  bearing 
system  it  is  the  internal  temperature  scheme  which  is  employed. 

If  IT2  is  left  blank  or  set  to  zero,  the  number  of  internal 
iterations  is  limited  to  twenty  (20). 

I tern  9 : "Accuracy  for  Internal  Thermal  Solution"  (EP2) . 

The  use  of  EP2  is  explained  in  Section  3.6.  If  EP2  is  left 
blank  dr  set  to  zero  (0),  a default  value  of  0.001  is  used. 

Item  10 : "Starting  Time"  (START)  is  a time  in  seconds  T at  which 
the  transient  solution  begins;  usually  set  to  zero  (0).  s 

Item  11;  "Stopping  Time"  (STOP)  is  the  time  in  seconds  at 
which  the  transient  solution  terminates , Tf  . The  transient 
solution  will  generate  a history  of  the  system  performance 
which  will  encompass  a total  elapsed  time  of 

(Tf  - Tg)  seconds 

Item  12:  "Calculation  Time  Step"  (STEPIN) . The  transient 
internal  solution  scheme  solves  the  system  of  equations 

t,  ♦ AT  (3.11) 

PC  pV 

AT  = STEPIN 

The  user  may  specify  STEPIN.  If  left  blank  or  set  to  zero  (0), 
the  Program  calculates  an  appropriate  value  for  STEPIN  using 
the  procedure  described  in  Section  3.6.  . 

Item  13:  "Time  Interval  Between  Printed  Temperature  Maps" 

(TTIMlfJ  seconds.  The  user  must  specify  the  length  of  time  in  seconds 
which  will  elapse  between  each  printing  of  the  temperature 
map.  The  interval  will  always  be  at  least  as  large  as  the 
"calculation  timestep"  (STEPIN) . 

Item  14 : "Time  Interval  Between  Calls  of  the  Shaft  Bearing 
Portion  of  the  Program"  (BTIME) . BTIME  will  always  have  a 
value  larger  than  or  equal  to  (STEPIN)  even  if  the  user 
inadvertently  inputs  a shorter  interval.  Computational  time 
savings  result  if  BTIME  is  greater  than  STEPIN,  however, 
accuracy  might  be  lost. 


'k+1 
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3.4.2  Card  Type  T2 


Card  Type  T2  is  required,  one  card  for  each  bearing  if 
no  thermal  analysis  is  being  performed.  The  temperature  data 
is  used  within  the  shaft-bearing  analysis  portion  of  the  pro- 
gram to  fix  temperature  dependent  properties  of  the  lubricant 
in  which  case  the  inner  race,  outer  race,  lubricant  bulk 
cavity  and  flange  temperatures  are  used.  Note  the  flange 
numbering  scheme  depicted  in  Fig.  3.5.  The 

tapered  roller  bearing,  inner  ring  flange  is  considered  flange  No.l. 
The  assembly  component  temperatures  at  each  bearing  location 
are  used  in  the  analysis  which  calculates  the  change  in  bearing 
diametral  clearance  from  "off  the  shelf"  to  operating  conditions. 


3.4.3  Card  Type -T  3 

In  the  steady  state  analysis  this  card  is  used  to  input 
initial  guesses  of  individual  nodal  temperatures  for  unknown 
nodes  as  well  as  the  constant  temperatures  for  known  nodes, 
such  as  ambient  air  and/or  an  oil  sump. 

In  the  transient  analysis.  Card  Type  T3  is  used  to  input 
the  nodal  temperatures  of  all  nodes  at  (START)  = Ts  i.e. 
at  the  initiation  of  the  transient  solution. 

3.4.4  Card  Type_T4 

With  this  card,  node  numbers  are  assigned  to  the  components 
of  each  bearing,  one  card  per  bearing.  With  this  information 
the  proper  system  temperatures  are  carried  into  each  respective 
bearing  analysis.  The  inner  race  and  inner  ring  node  numbers 
may  or  may  not  be  the  same  at  the  user's  discretion.  Similarly 
the  outer  race  and  outer  ring  node  numbers  may  or  may  not  be 
the  same. 

3.4.5  Card  Type  _T5 

The  shaft  bearing  system  analysis  accounts  for  frictional 
heat  generated  at  four  locations  in  the  bearing,  i.e.  at  the 
inner  race,  the  outer  race,  between  the  cage  rail  and  ring 
land,  and  in  the  bulk  lubricant  due  to  drag.  The  heat  genera- 
ted at  the  hydrodynamic  cage-rolling  element  contact  is  added 
to  the  bulk  lubricant.  Heat  generated  at  the  flange  is  not 
presently  considered.  This  card  allows  the  heat  generated  to 
be  distributed  equally  to  two  nodes.  For  instance  the  heat 
generated  at  the  inner  race-rolling  element  contact  should  be 
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distributed  half  to  the  rolling  element  and  half  to  the  inner 
race.  The  heat  developed  between  the  cage  and  inner  ring 
land  may  be  distributed  half  to  the  inner  ring  and  half  to 
the  cage  if  a cage  node  has  been  defined,  otherwise,  half  to 
the  bulk  lubricant. 

3.4.6  Card  Type  T6 

This  card  specifies  the  node  numbers  and  the  heat  genera- 
tion rate  for  those  nodes  where  heat  is  generated  at  a constant 
rate  such  as  at  rubbing  seals  or  gear  contacts. 

3.4.7  Card  Type  T7 

This  card  type  is  used  to  input  the  numerical  values 
of  the  various  heat  transfer  coefficients  which  appear  in  the 
equations  for  heat  transfer  by  conductivity,  free  convection, 
forced  convection,  radiation  and  fluid  flow.  Up  to  ten 
coefficients  of  each  type  may  be  used.  Separate  values  of 
each  type  of  coefficient  are  assigned  an  index  number  via 
card  T7  and  in  describing  heat  flow  paths  (Card  Type  T8  below) 
it  is  necessary  only  to  list  the  index  number  by  which  heat 
transfers  between  node  pairs. 

Indices  1-10  are  reserved  for  the  conduction  coefficient 
X,  11-20  for,the  free  convection  parameters,  21-30  for  forced 
convection,  31-40  for  emissivity  and  41-50  for  fluid  flow 
(product  of  specific  heat,  density  and  volume  flow  rate). 

As  an  example,  for  heat  transfer  by  conduction  with 
coefficient  X of  53.7  watts/M°C  one  could  prepare  a card  type 
T7  with  the  digit  1 punched  in  column  10  and  the  value  53.7 
punched  in  the  field  corresponding  to  card  columns  11-20. 

If  a conduction  coefficient  of  46.7  were  applicable  for 
certain  other  nodes  in  the  system  one  could  punch  an  additional 
card  assigning  index  No.  2 to  the  value  X = 46.7  by  punching 
a "2"  in  card  column  10  and  46.7  anywhere  within  card  columns 
11-20. 

Rather  than  inputting  constant  forced  convection  coefficients, 
optionally, these  coefficients  can  be  calculated  by  the  program 
in  one  of  three  ways.  If  the  calculation  option  is  exercised 
a pair  of  cards  is  used  in  place  of  a single  card  containing 
a fixed  value  of  a The  contents  of  the  pair  of  cards  depends 
upon  which  of  the  three  optional  methods  are  used. 

Option  1)  a is  independent  of  temperature  but  is  calculated 
as  a function  of  the  Nusselt  number  which  in 
turn  is  a function  of  the  Reynolds  number  Re, 
the  Prandtl  number  P as  follows,  (cf.  (17)  ) 
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where 


cl  ~ Nu  Ao:q/L 

Nu  = aRbPc 
e r 

A is  the  lubricant  conductivity,  L is  a 
characteristic  length  (with  a unit  of  meters) 
and  c,  a and  b are  constants. 

Option  2)  a is  a function  only  of  fluid  dynamic  viscosity 
and  viscosity  is  temperature  dependent. 

d 

a = cn 

Option  3)  a is  a function  of  the  Nusselt,  Reynolds  and 
Prandtl  numbers  and  viscosity  is  temperature 
dependent. 

3 .4  . 8 Card  Type  T8 

This  card  defines  the  heat  flow  paths  between  pairs  of 
nodes.  Every  node  must  be  connected  to  at  least  one  other 
node,  i.e.,  two  or  more  independent  node  systems  may  not  be 
solved  with  a single  Program  execution. 

The  calculation  of  heat  transfer  areas  is  based  on  lengths, 

Lj  and  L2  input  using  Card  Type  T8.  Additionally,  the 
type  of  surface  for  which  the  area  is  being  calculated  is  indicated 
by  the  sign  assigned  to  the  heat  transfer  coefficient  index. 

If  the  surface  is  cylindrical  or  circular  the  index  should 
be  positive,  if  the  surface  is  rectangular  the  index  should 
be  input  as  a negative  integer. 

In  the  case  of  radiation  between  concentric  axially 
symmetric  bodies,  L3  is  the  radius  of  the  larger  body.  For 
radiation  between  two  parallel  flat  surfaces  or  for  conduction 
between  nodes,  L3  is  the  distance  between  them. 

Fluid  flow  heat  transfer  accounts  for  the  energy  which 
the  fluid  transports  across  a node  boundary.  Along  a fluid 
node  at  which  convection  is  taking  place,  the  temperature 
varies.  The  nodal  temperature  which  is  output  is  the  average 
of  the  fluid  temperature  at  the  output  and  input  boundaries. 

If  the  emerging  temperature  of  the  fluid  is  of  interest,  it 
is  necessary  to  have  a fluid  node  at  the  fluid  outlet.  At 
this  auxiliary  node  only  fluid  flow  heat  transfer  occurs  and 
the  fluid  temperature  would  be  constant  throughout  the  node. 

Thus  the  true  fluid  outlet  temperature  will  be  obtained. 
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Conduction  of  heat  through  a bearing  is  controlled  by  index 
51.  The  actual  heat  transfer  coefficient  which  contains  a conduc- 
tivity, area  and  a path  length  term  is  calculated  in  the  bearing 
portion  of  the  program.  The  term  is  based  upon  conduction  through 
an  average  outer  race  and  inner  race  rolling  element  contact. 


3.4.9  Card  Type  T9 

This  card  inputs  data  required  to  calculate  the  heat 
capacity  of  each  node  in  the  system.  This  card  type  is  required 
only  for  a transient  analysis. 

3.5  DATA  SET  IV  - SHAFT  INPUT  DATA 

The  shaft-bearing  analysis  requires  all  loading  to  be 
applied  to  the  shaft.  The  loads  applied  to  each  bearing  are 
a product  of  the  shaft-bearing  solution.  There  is  no  need 
for  the  user  to  solve  the  statically  determinate  or  indeter- 
minate system  for  bearing  loads.  Even  if  a single  bearing  is 
being  analyzed, with  the  applied  load  acting  through  the  center 
of  the  bearing,  data  for  a dummy  shaft  must  be  supplied. 

In  the  analysis  the  housing  is  assumed  to  be  rigid.  . 
Provision  has  been  allowed  to  input  data  for  housing  radial 
and  angular  spring  characteristics.  However,  this  has 
been  done  for  future  consideration  of  an  elastic  housing  and 
is  therefore  currently  unavailable. 

The  shaft  input  data  consists  of  three  card  types: 

1)  Shaft  Geometry  and  Elastic  Modulus  Data 

2)  Bearing  Position  and  Mounting  Error  Data 

3)  Shaft  Load  Data 

3.S.1  Card  Type  SI 

This  card  type  is  used  to  describe  shaft  geometry  at  up 
to  twenty  locations  along  the  shaft.  The  user  must  place  his 
shaft  in  a cartesian  coordinate  system  with  the  end  of  the 
shaft  at  the  origin  and  with  the  shaft  lying  along  the  X-axis. 

The  shaft  may  have  stepwise  and  linear  diameter  variations. 
The  stepwise  variations  require  a single  card  which  specifies 
different  diameters  immediately  to  the  left  and  right  of  the 
relevant  X shaft  coordinate.  The  shaft  analysis  assumes  a linear 
diametral  variation  if  on  two  successive  cards,  i.e. 

two  successive  X coordinates,  the  diameter  to  the  right  of  the 
location  differs  from  the  diameter  to  the  left  of  the  location 


60 


of  the  following  card.  Complex  shaft  geometries  may  be  approxi- 
mated with  a set  of  linear  diameter  variations  spaced  at  close 
intervals . 

If  an  Elastic  Modulus  is  not  specified  at  the  designated 
input  location,  the  modulus  of  steel  is  assumed,  204083N/mm^. 

3.5.2  Card  Type  S2 

This  card  type  locates  the  bearing  inner  ring  on  the 
shaft  in  the  X-Y  and  X-Z  planes.  For  a ball  bearing,  the  X 
coordinate  specified  locates  the  inner  ring  center  of  curvature. 
For  cylindrical  roller  bearings  the  X coordinate  locates  the 
center  of  the  inner  race  roller  path. 

For  tapered  roller  bearings  in  the  strictest  sense  the 
X coordinate  locates  the  point  where  a line  from  the  roller 
center  of  gravity,  intersects  and  is  perpendicular  to  the  inner 
raceway,  with  all  bearing  end  play  removed.  It  is  sufficiently 
accurate  however  to  allow  the  X coordinate  to  locate  the  center 
of  the  inner  raceway. 

In  addition  to  specifying  bearing  location,  the  Type  2 
card  is  also  used  to  specify  housing  radial  and  angular 
mounting  errors.  As  mentioned  previously,  space  has  been  re- 
served for  inputing  housing  radial  and  angular  spring  charac- 
teristics, however,  these  characteristics  are  not  used  in  the 
system  analysis. 

Two  sets  of  Type  2 cards  may  be  required.  The  first 
set  is  always  required  and  defines  housing  alignment  errors  in 
the  shaft  X-Y  plane.  The  second  set  defines  the  housing  align- 
ment errors  in  the  shaft  X-Z  plane  and  is  required  only  if  non 
zero  errors  exist  for  the  particular  bearing  in  question. 

The  first  set  of  Type  2 cards  must  contain  a card  for  each 
bearing. 

3.5.3  Card  Type  S3 

Type  3 cards  are  used  to  specify  shaft  loadings  at  a 
given  X coordinate.  Loading  may  be  applied  in  the  X-Y  and 
X-Z  planes,  thus  requiring  two  distinct  sets  of  Type  3 cards. 
Applied  loads  may  have  the  form  of  concentrated  radial  forces, 
concentrated  moments,  linearly  distributed  radial  forces  and 
concentrated  axial  loads  which  may  be  eccentrically  applied. 

If  an  axial  load  is  eccentrically  applied,  the  moment  which 
results  must  not  be  separately  calculated  and  input  as  a con- 
centrated moment. 
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Variations  in  distributed  radial  loads  are  handled  at 
input  just  as  shaft  linear  diameter-  variations  are  handled. 

Note  that  each  set  of  Type  3 cards  must  be  followed  by  a 
blank  card. 

Also  note  that  in  order  for  symmetry  conditions  to  be 
considered  the  second  Type  3 card  must  be  void  of  any  loading 
data. 

4 . 0 COMPUTER  PROGRAM  OUTPUT 

4. 1 Introduction 

The  Program  Output  is  intended  to  provide  the  engineer  or 
designer  with  a complete  picture  of  the  shaft-bearing  system 
performance. 

In  addition  to  the  calculated  output  data,  the  input  data 
is  listed,  thus  producing  a complete  record  of  the  computer  run. 

A sample  output  is  included  in  Appendix  D.  The  multi- 
bearing supported  shaft  problem  demonstrates  the  full  capa- 
bilities of  SHABERTH  including  the  steady  state  thermal  and 
bearing  clearance  change  solutions. 

Key  output  items  are  discussed  briefly  below. 

4.2  Bearing  Output 

4. 2. 1.1  Linear  and  Angular  Deflections 

These  deflections  refer  to  the  bearing  inner  ring  relative 
to  the  outer  ring  and  are  defined  in  the  inertial  coordinate 
gys'tgni  of  Figure  2.4.  The  bearing  deflections  are  not  neces- 
sarily equal  to  the  shaft  displacements  since  the  bearing 
outer  ring  radial  or  angular  mounting  errors  may  be  specified 
as  non- zero  input. 
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4. 2. 1.2  Reaction  Forces  and  Moments 


These  values  reflect  bearing  reactions  to  shaft  applied 
loading  and  outer  ring  mounting  errors. 

When  the  bearing  inner  ring  has  achieved  an  equilibrium 
position,  the  summation  of  all  bearing  reaction  loads  should 
numerically  equal  the  shaft  applied  loading.  When  the  level 
of  solution  indicated  by  "NPASS"  - 2 is  employed,  as  discussed 
in  Section  5,  differences  between  shaft  applied  and  bearing 
reaction  loads  will  exist  but  will  typically  be  less  than  10%. 
This  difference  is  a consequence  of  friction  forces  contributing 
to  the  reaction  loads  whereas  the  inner  ring  equilibrium  posi- 
tion has  been  determined  considering  elastic  contact  forces  only. 

4.2.2  Fatigue  Life  Data 

The  L10  fatigue  life  of  the  outer  and  inner  raceways  as 

well  as  the  bearing  are  presented.  The  bearing  life  represents 
the  statistical  combination  of  the  two  raceway  lives.  These 
lives  reflect  the  combined  effects  of  the  lubricant  film  thick- 
ness and  material  life  factors.  The  lubricant  film  thickness 
life  factor  is  described  in  detail  in  Section  3. 


4 . 2 . 2 . 1 h/ sigma 

The  ratio  h/c,  also  referred  to  as^,  is  printed  for  the 
most  heavily  loaded  rolling  element.  The  variable  h,  represents 
the  EHD  plateau  film  thickness  with  thermal  and  starvation 
effects  considered.  The  variable  a,  represents  the  composite 
root  mean  square  surface  roughness  of  the  rolling  element  and 
the  relevant  raceway. 
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4. 2. 2, 2 Life  Multipliers 


4. 2. 2. 2.1  Lubrication  - This  life  multiplier  is  a function  of 

h/^  at  each  concentrated  contact  and  is  in  the  form  of  a derating 
factor.  Its  value  ranges  from  0.479  for  h/tj*  = 0 to  1.0  at  h/^ 

^.4.  Since  the  lubricant  life  multiplier  is  decremental  the 
normal  multiple  of  3 used  for  thick  film  lubrication  must  be 
multiplied  by  the  material  life  factor  normally  used  and  this 
product  should  be  specified  at  input.  This  subject  is  covered 
in  more  detail  in  Section  3.3. 1. 

4. 2. 2. 2. 2 Material  - This  output  simply  reflects  the  input 
value,  Again,  it  is  covered  in  Section  3. 

4.2.3  Temperatures  Relevant  to  Bearing  Performance 

These  temperatures  fully  describe  the  temperature  con- 
ditions which  affect  the  performance  of  a given  bearing.  If 
one  of  the  temperature  mapping  options  is  used,  the  tempera- 
tures printed  reflect  the  results  of  the  particular  option. 

If,  neither  temperature  option  was  used,  the  list  is  simply  a 
repeat  of  the  input  data.  Note  that  there  are  separate  tem- 
peratures for  outer  and  inner  raceways  and  flanges  and  ring 
temperatures.  The  raceway  and  flange  temperatures  are  used  to 
determine  lubricant  properties.  The  ring  temperatures  are 
used  in  the  bearing  dimension  change  analysis.  The  raceway, 
flange  and  ring  temperatures  may  be  the  same  value. 

4.2.4  Frictional  Heat  Generation  Rate  and  Bearing  Friction 
Torque 

4. 2.4.1  Frictional  Heat  Generation  Rate 


The  various  sources  of  frictional  heat  generated  within 
the  bearings  are  listed.  The  values  printed  for  "OUTER  RACE, 
OUTER  RING  FLANGES,  INNER  RACE,  INNER  RING  FLANGES,  R.E.DRAG 
AND  R.E.  CAGE"  represent  the  sum  of  the  generated  heats  for  all 
rolling  elements.  Additionally,  the  heats  printed  for  the 
outer  and  inner  raceways  and  flanges,  plus  the  rolling  element 
cage,  reflect  the  friction  developed  outside  the  concentrated 
contacts,  i.e.,.the  HD  friction  as  well  as  the  EHD  friction 
developed  within  the  concentrated  contacts.  The  raceway  and 
flange  data  also  includes  any  heat  generated  as  a consequence 
of  asperity  contacts.  "R.  E.  DRAG"  should  be  interpreted  as 
the  heat  resulting  from  lubricant  churning  as  the  rolling 
elements  plow  through  the  air-oil  mixture. 
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4. 2. 4. 2  Torque 


The  torque  value  is  calculated  as  a function  of  the  total 
generated  heat  and  the  sum  of  the  inner  and  outer  ring  rotational 
speeds.  The  intent  is  to  present  a realistic  value  of  the  torque 
required  to  drive  the  bearing.  Under  conditions  of  inner  ring 
rotation  the  torque  value  reflects  the  torque  required  to  drive 

the  inner  ring  Singe  the  inner  ring  torque  also  includes  that  ..friction 
torque  required  to  impart  an  angular  velocity  to  the  lubricant 
in  the  bearing.  A considerable  portion  of  the  lubricant  will 
come  to  rest  within  the  housing  and  not  at  the  outer  ring. 

Thus  the  measured  outer  ring  torque  may  not  equal  the  torque 
at  the  inner  ring. 

4.2.5  EHD  Film  and  Heat  Trans  ter  Data 

4. 2.5.1  EHD  Film  Thickness 

These  values  refer  to  the  calculated  EHD  plateau  film 
thickness  at  both  contacts  of  the  most  heavily  loaded  rolling 
element  and  include  the  effects  of  the  thermal  and  starvation 
reduction  factors. 

4. 2. 5. 2 Starvation  Reduction  Factor 

These  factors  give  for  the  inner  and  outer  ring  contacts, 
the  reduction  in  EHD  film  thickness  due  to  lubricant 
film  starvation  according  to  the  methods  of  Chiu,  (11). 

These  factors  pertain  to  the  EHD  film  thickness  for  both 
the  inner  and  outer  race  contacts  of  the  most  heavily  loaded 
rolling  elements,  but  are  applied  to  the  respective  inner  and 
outer  race  film  thickness  for  each  rolling  element  in  the  bear- 
ing. 

4. 2. 5. 3 Thermal  Reduction  Factor 

These  factors  are  calculated  according  to  the  methods  of 
Cheng,  (10)  and  pertain  to  the  EHD  film  thickness  for  both 
the  inner  and  outer  race  contacts  of  the  most  heavily  loaded 
rolling  elements,  but  are  applied  to  the  respective  inner  and 
outer  race  film  thickness  for  each  rolling  element  in  the  bear- 
ing. 

4. 2. 5.4  Meniscus  Distance 

These  factors  are  calculated  according  to  the  methods  of 
Chiu,  (11)  and  pertain  to  the  EHD  film  thickness  for  both 
the  inner  and  outer  race  contacts  of  the  most  heavily  loaded 
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rolling  elements,  but  are  applied  to  the  respective  inner  and 
outer  race  film  thickness  for  each  rolling  element  in  the  bear- 
ing. 

4. 2. 5. 5 Raceway- Rolling  Element  Conductivity 

These  data  reflect  the  amount  of  heat  transfer  between 
rolling  element  and  raceway  for  each  degree  centigrade  differ- 
ence between  the  two  components . These  data  reflect  the  average 
of  all  outer  and  inner  contacts  respectively. 

4.2.6  Fit  and  Dimensional  Change  Data 


4. 2.6.1  Fit  Pressures 

These  data  refer  to  the  pressures  built  up  as  a conse- 
quence of  interference  fits  between  shaft  and  inner  ring  and 
housing  and  outer. ring.  Pressures  are  presented  both  for  the 
standard  cold-static  condition  (16°C)  and  at  operating 
conditions . 

4. 2. 6. 2 Speed  Giving  Zero  Fit  Pressure  (Between  the  shaft  and 
inner  ring! 

This  is  a calculated  value  based  upon  operating  conditions 
and  provides  a measure  of  the  adequacy  of  the  initial  shaft 
fit. 

4. 2.6.3  Clearances 

"Original"  refers  to  cold  unmounted  clearance  which  is 
specified  at  input  if  the  diametral  clearance  change  analysis 
is  executed.  "Change"  refers  to  the  change  in  diametral 
clearance  at  operating  conditions  relative  to  .the  cold  unmounted 
condition.  A minus  sign  indicates  a decrease  in  clearance. 
"Operating"  refers  to  the  clearance  at  operating  conditions. 

For  all  types  of  ball  bearings  the  decrease  in  clearance  can 
be  combined  with  the  initial  diametral  clearance,  and  the  free 
operating  contact  angle  at  operating  conditions  may  be 
calculated.  Note  that  the  change  in  clearance  should  be  com- 
pared against  the  diametral  play  of  the  split  inner  ring  ball 
bearing  in  order  to  determine  if  the  possibility  for  three 
point  contact  exists.  The  Program  does  not  account  for  three 
point  contacts  even  though  the  change  in  clearance  might  suggest 
that  three  point  contact  is  obtained. 
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4.2.7  Lubricant  Temperatures  and  Physical  Properties 


The  lubricant  properties,  particularly  the  dynamic  vis- 
cosity and  to  a lesser  degree  the  pressure  viscosity  coefficient, 
are  heavily  temperature  dependent.  These  factors  enter  the  EHD 
film  thickness  calculation  and  the  HD  and  EHD  friction  models. 

The  lubricant  is  assumed  to  be  at  the  same  temperature  as  the 
relevant  raceway.  As  noted  elsewhere, these  temperatures  may 
be  either  input  directly  or  calculated  by  the  Program. 

The  physical  properties  printed  are  self  explanatory. 

The  units  are  enumerated. 

4.2.8  Cage  Data 

4. 2. 8.1  Cage-Land  Interface 

The  cage  data  indicates  the  performance  parameters  at  the 
interface  between  the  cage  rail  and  the  ring  land  on  which  the 
cage  is  guided.  The  torque,  heat  rate  and  separating  force 
require  no  explanation.  The  eccentricity  ratio  defines  the 
degree  to  which  the  cage  approaches  the  ring  on  which  it  is 
guided  at  the  point  of  nearest  approach.  The  radial  displace- 
ment of  the  cage  relative  to  the  bearing  axis  is  divided  by  one 
half  the  cage-land  diametral  clearance.  An  eccentricity  ratio 
of  one  indicates  cage- land  contact.  A ratio  of  zero  indicates 
that  the  cage  rotation  is  concentric  with  the  bearing  axis. 

Only  the  cage- land  and  rolling  element  pocket  forces  are 
considered  in  determining  the  cage  eccentricity.  The  cage 
weight  and  centrifugal  force  which  result  from  the  eccentricity 
although  available, are  not  considered  in  the  analysis.  The 
omission  of  these  considerations  helps  reduce  convergence 
problems . 


4. 2. 8. 2 Cage  Speed  Data 

Cage  speed  data  presents  the  comparison  between  the  cage 
speed  calculated  based  upon  the  quasidynamic  equilibrium  con- 
siderations and  the  speeds  calculated  with  raceway  control 
theory  for  ball  bearings  and  the  epicyclic  speeds  of  the 
roller  bearing  components. 

4. 3 Rolling  Element  Output 

4.3.1  Rolling  Element  Kinematics 
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4. 3. 1.1  Rolling  Element  Speeds 

All  of  the  rolling  element  speeds  tend  to  vary  from  posi- 
tion to  position  when  the  bearing  is  subjected  to  combined 
loading. 

The  total  rolling  element  speed  is  with  reference  to  the 
cage  and  represents  the  vector  sum  of  the  three  orthogonal 
components . 

4. 3.1. 2 Speed  Vector  Angles 

The  rolling  element  speed  vector  angles,  Arctan  (JUy/iQx) 
and  Arctan  are  presented  in  order  to  show  a clearer 

picture  of  the  predicted  ball  kinematics.  The  ball  speed 
vector  tends  to  become  parallel  with  the  bearing  X axis  with 
increasing  shaft  speed  and  decreasing  contact  friction. 

4.3.2  Rolling  Element  Raceway  Loading 

4. 3. 2. 1 Normal  Forces 

The  normal  forces  acting  on  each  rolling  element  are 
printed.  The  rolling  element  race  normal  forces  are  self 
explanatory.  The  cage  force  is  calculated  only  when  the 
friction  solution  is  employed  and  is  always  directed  along  the 
rolling  element  Z axis.  If  the  rolling  element  orbital  speed 
is  positive,  a positive  cage  force  indicates  that  the  cage  is 
pushing  the  rolling  element,  tending  to  accelerate  it.  Cage 
force  is  a function  of  rolling  element  position  within  the 
cage  pocket.  Its  magnitude  is  derived  using  hydrodynamic  lubri- 
cation assumptions,  when  the  distance  between  the  rolling  element 
and  cage  web  is  large,  and  EHD  assumptions  when  the  separation 
is  of  the  order  of  the  EHD  film  thickness  or  when  rolling 
element- cage  web  interference  exists. 

4. 3.2.2  Hertz  Stress 

The  stress  printed  represents  the  maximum  normal  stress 
at  the  center  of  each  ball  race  contact  or  at  the  most  heavily 
loaded  slice  of  the  roller  raceway  contact. 

4. 3. 2.3  Load  Ratio  Qa3p/Qt0t  - 

If  the  EHD  film  thickness  is  small  compared  to  the  RMS 
composite  rolling  element-race  surface  roughness,  the  rolling 
element-race  normal  load  will  be  shared  by  the  EHD  film  and 
asperity  contacts.  The  load  ratio  reflects  the  portion  of  the 
total  load  carried  by  the  asperities. 
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4. 3. 2. 4 Contact  Angles 

A ball  bearing,  subject  to  axial  loading,  misalignment 
or  mounted  such  that  the  inner  ring  is  always  displaced  axially 
relative  to  the  outer  rings,  (i.e.  a duplex  set  of  angular 
contact  ball  bearings)  will  have  non- zero  contact  angles.  At 
low  ball  orbital  speeds  the  inner  and  outer  race  angles  will  be 
substantially  the  same.  At  high  speeds,  ball  centrifugal  force 
will  cause  the  outer  race  contact  angle  to  be  less  than  the 
inner  race  angle. 


4.3,3  Roller  End-Flange  Contact  Data 

For  a tapered  roller  bearing  a single  set  of  roller 
end-flange  data  is  printed.  For  the  cylindrical  roller  bearing, 
which  may  have  up  to  four  flanges,  the  program  examines  the 
data  and  prints  the  results  for  the  most  heavily  loaded  outer 
and  inner  ring  flanges.  In  addition  to  the  data  listed  below 
for  a cylindrical  roller,  the  semimajor  contact  axis  as  well  as 
the  concentrated  contact  and  hydrodynamic  heat  generation  rates 
are  printed  for  the  tapered  roller. 

4. 3. 3.1  Normal  Force 


The  interference  between  the  roller  end  and  flange  is 
determined  from  the  solution  for  the  relative,  ring  and  rolling 
element  positions.  Hertz  theory  is  used  to  calculate  the  load 
which  results  from  this  interference. 

4. 3. 3. 2 Hertz  Stress 

The  Hertz  stress  printed  is  the  maximum  normal  stress 
which  occurs  in  the  contact. 

4. 3. 3. 3 EHD  Film  Thickness 


The  plateau  film  thickness  is  calculated  using  the 
Archard  Cowking  equation  for  point  contacts.  This  result  is 
then  modified  to  account  for  starvation  (11)  and  thermal  (10) 
reduction  effects.  This  modified  film  thickness  is  printed. 

4. 3. 3.4  Sliding  Velocity 

The  sliding  velocity  is  defined  as  the  difference  be- 
tween the  flange  and  roller  end  linear  velocities  at  the  center 
of  the  contact. 
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4. 3. 3. 5  Rolling  Velocity 


Rolling  velocity  is  defined  as  one  half  the  sum  of  the 
roller  end  and  flange  linear  velocities  at  the  center  of  the 
contact. 


4. 3. 3. 6  Contact  Ellipse  Semiminor  Axis 

To  help  assess  the  severity  of  the  roller  end  flange 
contact  and  the  possibility  for  edge  loading,  the  semiminor 
axis  of  the  contort  is  nrinted. 


4.4  Thermal  Data 


As  in  the  case  for  bearing  output,  all  of  the  input  data  are 
printed.  The  calculated  output  data  are  presented  in  the  form 
of  a temperature  map  in  which  a node  number  and  the  respective 
node  temperature  appear.  The  appearance  of  the  steady  state 
and  transient  temperature  maps  are  identical.  The  transient 
temperature  map  also  includes  the  time  (T)  at  which  the 
temperature  calculations  were  made. 

4.5  Shaft  Data 

These  data  simply  reflect  the  input  information. 

4 . 6 Program  Error  Messages 

4.6.1  From  Subroutine  ALLT 

''Steady  State  Solution  with  CEP1)  degrees  accuracy  was  not 
obtained  after  CITl)  Iterations". 

This  message  pertains  to  the  external  temperature  itera- 
tion scheme  in  which  system  temperatures  and  bearing  generated 
heats  are  being  solved  for  an  equilibrium  condition. 

4.6.2  From  Subroutine  SHABE 

"It  was  not  possible  to  obtain  the  change  of  clearance 
with  an  accuracy  of  CERFIT)  times  the  rolling  element  diameter 
in  CITFIT)  iteration's^ 
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This  message  pertains  to  the  bearing  diametral  clearance 
change  iteration  scheme.  The  solution  may  be  converging  in 
which  case  the  number  of  iterations  (ITFIT)  should  be  increased. 
This  can  be  checked  with  an  NPRINT  - 1 intermediate  printout. 

The  intermediate  print  may  indicate  that  the  solution  is 
oscillating.  The  most  likely  cause  of  oscillation  is  the 
alternate  prediction  of  bearing  preload  with  all  rolling  elements 
loaded,  and  then  in  the  next  iteration,  only  a subset  of  the 
rolling  elements  loaded.  This  problem  can  usually  be  over- 
come by  either  of  two  methods. 

1)  In  subroutine  FIT  remove  the  GO  TO  20  statement. 

This  will  cause  the  inner  ring  load  distribution 

to  have  no  effect  on  the  change  in  diametral  clear- 
ance . 

2)  The  solution  can  be  damped  by  redefining  the  solution 
damping  factor  FA,  such  that  it  would  take  on  a value 
O^FA^l.  FA  is  presently  set  to  1 in  subroutine  SHABE. 

If  this  damping  technique  is  used,  the  number  of  FIT 
iterations  should  be  increased  as  the  value  of  FA  is 
decreased.  An  upper  limit  of  10  iterations  is  recom- 
mended. 

4.6.3  From  Subroutine  SOLVXX 

1)  "SINGULAR  SET  OF  EQUATIONS" 

This  message  might  occur  when  the  thermal  input  data 
is  not  input  properly. 

2)  "THE  LIMIT  FOR  NUMBER  OF  ITERATIONS  IS  REACHED" 

This  message  might  occur  either  during  a steady  state 
temperature  solution  or  bearing  solution.  Before 
increasing  the  number  of  iterations  check  the  equa- 
tion residue  values.  If  they  are  low,  the  solution 
may  be  good  enough. 

3)  "THIS  IS  THE  BEST  WE  CAN  DO.  IT  MAY  BE  USEABLE" 

This  message  reflects  the  fact  that  the  next  iteration 
will  result  in  divergence.  The  iteration  procedure  is 
thus  terminated.  The  equation  residue  values  should 
be  checked.  If  they  are  low,  the  solution  may  be  useable. 


As  noted  above  the  occurrence  of  messages  2)  and  3),  do 
not  necessarily  mean  that  the  solution  is  not  good.  Generally 
the  messages  indicate  that  the  solution  has  not  converged  as 
tightly'  as  the  user  has  requested. 

Note:  The  XX  suffix  on  SOLV  specifies  the  version  of  SOLV 
contained  in  the  user's  program.  As  of  this  writing  the  current 
version  is  S0LV1 3.  The  suffix  is  changed  each  time  improve- 
ments are  made  which  require  a change  in  the  SOLVXX  calling 
sequence. 

4.6.4  From  Subroutine  INTFIT 

"Singular  matrix  on  tight  shaft  fit" 

"Singular  matrix  on  loose  shaft  fit" 

These  messages  reflect  an  error  in  the  input  data  usually 
as  a consequence  of  inconsistent  component  diameters,  such  as 
the  shaft  inside  diameter  being  greater  than  the  outside 
diameter.  ; 
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5.0  GUIDES  to  program  use 


The  Computer  Program  is  a tool.  As  with  any  tool  the 
results  obtained  are  at  least  partially  dependent  upon  the 
skill  of  the  user.  Certainly  the  economics  of  the  Program 
usage  are  highly  dependent  upon  the  user's  technical  need 
and  discriminate  use  of  Program  options. 

Some  general  guides  for  efficient  use  of  the  Program  are 
listed  below: 

1.  Attempt  to  use  the  lowest  level  of  solution  possible. 
For  instance  if  the  prime  object  of  a given  run  is  to 
obtain  bearing  fatigue  lives,  execute  only  the  elastic 
solution  (NPASS  = 0) . If  an  estimate  of  bearing  . 
frictional  heat  is  required,  execute  the  low  level 
friction  evaluation  (NPASS  = 1) . Execute  the  friction 
solution  (NPASS  = 2)  only  if  rolling  element  and  cage 

• kinematics  are  of  interest.  Execute  the  highest  level 
of  solution  (NPASS  = 3)  if  kinematics  are  of  interest 
and  the  bearing  reaction  loads  deviate  substantially 
from  the  shaft  applied  loading,  i.e.  a deviation 
greater  than  ten  percent. 

2.  Attempt  to  input  bearing  operating  diametral  clearance 
rather  than  calculate  it.  Or,  execute  the  diametral 
clearance  change  analysis  once  for  a group  of  similar 
runs  and  use  the  output  from  the  first  run  as  input 

to  the  subsequent  runs,  omitting  the  clearance  change 
analysis . 

3.  Attempt  to  input  accurate  operating  temperatures 
rather  than  calculate  them. 

4 . The  more  non-linear  the  problem  the  more  computer  time 
required  to  solve  it.  In  the  bearing  friction  solution 
large  coefficients  of  friction  seem  to  increase  the 
degree,  of  nonlinearity.  In  the  thermal  solutions,  if 
possible,  eliminate  nonlinearities  by  omitting 
radiation  terms  and  by  using  constant  rather  than 
temperature  dependent  free  and  forced  convection 
coefficients . 

5.  In  the  transient  thermal  solution,  space  the  calls  to 
the  shaft-bearing  solution  (BTIME-)  to  as  large  an 
interval  as  prudently  possible.  Be  careful  however, 
too  long  an  interval  will  produce  large  errors  in  heat 
rate  predictions. 
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6.  In  the  steady  state  thermal  analysis  attempt  to  esti- 
mate nodal  temperatures  on  a node  by  node  basis. 

Nodes  which  are  heat  sources  should  have  highel’ 
temperatures  than  the  surrounding  nodes. 

The  above  suggestions  are  intended  to  encourage  the  use 
of  the  Program  on  a cost  effective  basis.  The  intent  is  not 
to  discourage  the  use  of  important  program  capabilities  but 
to  emphasize  how  the  program  should  be  most  effectively  used. 

It  is  suggested  that  the  user  take  a simple,  axially 
loaded  ball  bearing  problem  and  execute  the  program  through 
the  full  range  of  options  beginning  with  a frictionless  solu- 
tion proceeding  to  the  three  levels  of  friction  solution  with 
a low  (0.013  and  high  (0.1)  friction  coefficient.  The  dia- 
metral clearance  change  analysis  and  the  thermal  solutions 
should  also  be  executed  on  an  experimental  basis.  This  exercise 
will  provide  the  user  with  some  insight  into' economics  of  the 
Program  usage  on  his  computer  as  well  as  the  results  obtained 
from  various  levels  of  solution  of  the  same  problem. 

It  is  also  suggested  that  a constant  user  of  the  program 
should  study  the  hierarchical  Program  flow  chart  Appendix  A, 
along  with  the  Program  listing  to  gain  an  appreciation  of  the 
program  complexity  and  the  flow  of  the  problem  solution.  The 
Program  is  comprised  of  many  small  functional  subroutines. 
Knowledge  of  these  small  elements  may  allow  the  user  to  more 
easily  piece  together  the  philosophy  of  the  total  problem 
solution. 
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APPENDIX  - A 


S K F COMPUTER  PROGRAM  AT77Y001  "SHABERTH" 


HIERARCHIAL  FLOW  CHART 
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APPENDIX  A 

S K F COMPUTER  PROGRAM  AT77Y001  FLOW  CHART 


Flow  Chart 


The  hierarchical  flow  chart  presents  the  program  structure, 
listing  the  program  elements  in  the  order  in  which  they  would 
be  called  to  solve  the  shaft-bearing  dynamic,  as  well  as  steady 
state  and  transient  temperature  distribution  problems.  The 
various  solution  loops  are  indicated,  as  well  as  notes  which 
indicate  the  functions  of  various  subroutine  groupings. 

Each  line  in  the  flow  chart  represents-  a program  element, 
subroutine,  function  or  the  main  program  ALWAYS.  The  call  of 
one  subroutine  by  another  is  denoted  by  indenting  the  called 
subroutine  relative  to  the  routine  doing  the  calling.  As  an 
example,  subroutine  SKF  calls  subroutines  FLAGS,  TYPE,  PROPST, 
LUPROP,  LUBCON,  DATOT,  CNVRT,  CONS  and  SPRING.  Subroutine 
CONS  calls  CONBRI,  BCON,  TCON,  CRCON  and  CONBRZ.  BCON  calls 
ABDEL.  Both  TCON  and  CRCON  call  ABDEL  and  SLICES. 

The  first  mention  of  a subroutine  within  the  flow  chart 
includes  the  entire  list  of  subordinate  program  elements.  At 
subsequent  calls  to  that  subroutine  the  list  of  subordinate 
elements  is  omitted-.  As  an  example  the  first  call  to  subroutine 
AXLBOJ  is  followed  by  the  subordinate  elements,  JMVIKT,  SNITMT, 
NUMLOS,  DUBSIM,  MEIE,  MEIL  and  SIMQ.  After  the  call  of  AXLBOJ 
from  INDEL,  the  subordinate  elements  are  not  listed  but  are,, 
nevertheless,  employed.  The  list  of  subordinate  program  ele- 
ments are  omitted  in  repeated  calls  of  subordinate  GUESS,  BEAR, 
SOLVXX  and  DELIV3  as  well  as  AXLBOJ. 

As  noted  earlier,  rolling  equilibrium  is  calculated,  first 
without,  then  if  required,  with  friction  forces  included. 

Whether  or  not  friction  is  considered  is  highlighted  with  the 
words  Frictionless  or  Friction. 

If  the  Program  is  too  large  to  fit  in  its  entirety  on  the 
user's  computer,  segments  of  the  program  may  be  "overlaid". 

For  this  purpose  the  Program  is  subdivided  into  ten  (10) 
moduleswhich  can  be  sequentially  "overlaid". 

The  Program  segments  SKF,  TEMPIN,  SHAFT  and  GUESS  all  per- 
form initiation  functions  and  with  the  exception  of  GUESS,  are 
called  only  once  per  program  execution. 

The  real  problem  solving  portion  of  the  program  is  embodied 
in*segment  ALLT.  Within  this  segment  the  shaft  bearing  solu- 
tion is  obtained  through  the  call  to  SHABE,  then  the  steady 
state  or  transient  temperature  distributions  are  obtained. 
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This  scheme  is  repeated  until  the  end  objective,  steady  state 
thermal  equilibrium  or  time  up  for  the  transient  scheme,  is 
realized. 

The  nonlinear  equation  solver  S0LV13  is  central  to  the 
program  and  deserves  special  discussion  as  related  to  the  flow 
chart.  The  first  call  to  SOLV'13  is  from  BEAR.  Only  for  this 
first  call  are  all  of  the  S0LV13  subordinate  subroutines 
listed  as  noted  earlier.  These  include  INSOLV,  EQS.  PARDER, 
SIMQ,  EQCHEK,  ERWRIT  and  ERCHEK.  In  the  subsequent  call  to 
S0LV13  in  which  the  steady  state  temperatures  are  being  calcu- 
lated, the  above  listed  subroutines  are  again  called  but 
these  calls  with  the  exception  of  EQS  are  not  listed  on  the 
flow  chart. 

EQS  is  the  name  given  by  S0LV13  to  a subroutine  which  sets 
up  the  system  of  equations  to  be  solved.  EQS  is  brought  into 
S0LV13  through  the  argument  list.  When  the  bearing  equations 
are  being  solved,  subroutine  BRGGEQ  is  brought  into  S0LV13 
and  within  S0LV13  is  referenced  by  the  name  EQS.  When  the  heat 
transfer  equations  are  being  solved  as  a consequence  of  the 
call  of  SOLV'13  from  ALLT,  NET  is  brought  into  SOLV'13  and  is 
referenced  as  EQS. 
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Evaluate  cage  speed,  ball-cage  normal  and  friction 
forces  and  neat  rates  plus  cage  rail-ring  land 
normal  and  friction  forces  plus  heat  generation  rates. 


Evaluate  ball,  position  dependent,  accelerations. 


Calculate  bearing  fatigue  life 


Calculate  ball-raceway  heat  transfer  coefficients 


Evaluate  roller -raceway  and  roller  - flange  contact 
loads 


Evaluate  roller  raceway  EHD  film  thickness  inlet  and 
concentrated  contact  friction  forces  and  heat  generation 
rate 


Evaluate  rol ler- flange  EHD  film  thickness  inlet  and 
concentrated  contact  friction  forces  and  heat  genera- 
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Evaluate  tapered  roller  bearing  cage  speed,  cage 
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Solve  shaft,  rolling  element  and  cage  equilibrium 
equations  with  friction  included 


Write  shaft-bearing  output 


Solve  for  the  steady  state  temperature  distribution. 
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desired . 
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transient  temperature  distribution.  Write  shaft- 
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APPENDIX  - B 


HEAT  TRANSFER  COMPUTATION  NOTES 


APPENDIX  B 


HEAT  TRANSFER  INFORMATION 


B.l  BACKGROUND 

The  temperature  portion  of  Program  AT77Y001  is  designed 
to  produce  temperature  maps  for  an  axisymmetric  mechanical  system 
of  any  geometrical  shape.  The  mechanical  system  is  first  approxi- 
mated by  an  equivalent  system  comprising  a number  of  elements  of 
simple  geometries.  Each  element  is  then  represented  by  a node 
point  having  either  a known  or  an  unknown  temperature.  The  environ- 
ment surrounding  the  system  is  also  represented  by  one  or  more 
nodes.  With  the  node  points  properly  selected,  the  heat  balance 
equations  can  be  set  up  accordingly  for  the  nodes  of  unknown 
temperature.  These  equations  become  non-linear  when  there  is 
convection  and/or  radiation  between  two  or  more  of  the  node  points 
considered.  The  problem  is  therefore  reduced  to  solving  a set 
of  linear  and/or  non-linear  equations  for  the  same  number  of 
unknown  nodal  temperatures.  It  is  obvious  that  the  success  of 
the  approach  depends  largely  on  the  physical  subdivision  of 
the  system.  If  the  subdivision  is  too  fine,  there  will  be  a large 
number  of  equations  to  be  solved;  on  the  other  hand,  if  the 
subdivision  is  too  crude,  the  results  may  not  be  reliable. 

In  a system  consisting  of  rolling  bearings,  for  the  sake 
of  simplicity,  the  elements  considered  are  usually  axially 
symmetrical,  e.g. , each  of  the  bearing  rings  can  be  taken  as  an 
element  of  uniform  temperature.  For  an  element  which  is  not 
axially  symmetrical,  its  temperature  is  also  assumed  to  be 
uniform  and  its  presence  is  assumed  not  to  distort  the  uniformity 
in  temperature  of  a neighboring  element  which  is  axially  symmetri- 
cal. That  is,  the  non- symmetrical  element  is  represented  by 
an  equivalent  axially  symmetrical  element  with  approximately 
the  same  surface  area  and  material  volume.  This  kind  of  approxi- 
mation may  seem  to  be  somewhat  unrealistic,  but  with  properly 
devised  equivalent  systems,  it  can  be  used  to  solve  complicated 
problems  with  results  satisfying  some  of  the  important  engineer- 
ing requirements. 

The  computer  program  can  solve  the  heat-balance  equations 
for  either  the  steady  state  or  the  transient  state  conditions 
and  produce  temperature  maps  for  the  mechanical  system  when  the 
input  data  are  properly  prepared. 

B.2  BASIC  EQUATIONS 

B.  2 . 1 Heat  Conduction 

The  rate  of  heat  flow  qci } j (W}  that  is  conducted  from  node 
i to  node  j may  he  expressed  hy, 
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and  tj  are  the  temperatures  at  i and  j,  respectively,  Aj^j 
the  areanormal  to  the  heat  flow,  (m^)  Lij  the  distance 
(m)  and  Ajj  the  thermal  conductivity  between  i and  j , (W/m°C)  . 

Assuming  that  the  structure  between  point  i and  j is 
composed  of  different  materials,  an  equivalent  heat  conductivity 
may  be  calculated  as  follows: 


The  calculation  of  the  areas  will  be  discussed  in  Section  B.2.5. 
B.2.2  CONVECTION 

The  rate  of  heat  flow  that  is  transferred  between  a solid 
structure  and  air  by  free  convection  may  be  expressed. by 
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where 


- a 


vi, J i,J  i,J  ' i J 


t,  - tj  I 1,25  • SIGN  (t.  - t.) 


SIGN  = 


> ° 

ti)<  0 


in  which 


• 2 2 12^ 

♦10  W/m  - (degC)  for  hot  surfaces  facing  upward 

and  cold  surfaces  facing  downward 

• 10~2  W/m2-  (degC)1,25  for  hot  surfaces  facing  downward 

and  cold  surfaces  facing  upward 

• 10~2  W/m2-  (degC)**2**  for  vertical  surfaces 


For  other  special  conditions,  a-jj  must  be  estimated  by  referring 
to  heat  transfer  literature. 


The  rate  of  heat  flow  that  is  transferred  between  a solid 
structure  and  a fluid  by  forced  convection  may  be  expressed  by 


q..-  a.  .A.  .ft.  - t . ) 

ni,3  i,j  i,l  v i r 


in  which  is  the  heat  transfer  coefficient. 

Now,  with  a=  Ojj  , introduce  the  Nusselt  number 


the  Reynolds  number 
and  the  Prandtl  number 


where 
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L is  a characteristic  length  which  is  equal  to  the  diameter 
in  the  case  of  a cylindrical  surface  and  is  equal  to  the 
plate  length  in  case  of  a flat  surface  (m) 
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U is  a characteristic  velocity  which  is  equal  to  the  difference 
between  the  fluid  velocity  at  some  distance  from  the  surface 
and  the  surface  velocity  (fo/sec  ) 

X is  the  fluid  thermal  conductivity (W/M°  C) 

v is  the  fluid  kinematic  viscosity (M^/sec) 

p is  the  fluid  density  (kg/m  ) 

c is  the  fluid • specif ic  heat(J/kg°c) 

sr  .i  * ' 

For  given  values  of  Rg  and  Pp  the  Nusselt  number  and 
thus  the  heat  transfer  coefficient  may  be  estimated  from  one 
of  the  following  expressions: 

Laminar  flow  along  a flat  plate:  Re  < 2300 

Nu  = 0.323  fTe  • 

Laminar  flow  of  a liquid  in  a pipe: 

Nu  - 1.36  3/Re  • Pr  (f) 

where  D is  the  pipe  diameter'  and  L the  pipe  length 
Turbulent  flow  of  a liquid  in  a pipe:  * 

Nu  = 0.027  • Re°’8  . 

Gas  flow  inside  and  outside  a tube: 

Nu  = 0.3  R°’57 


" Liquid  flow  outside  a tube: 


N 


u 


0.6  R, 


0.5 


p0  • 31 
r 


Forced  free  convection  from  the  outer  surface  of  a 
rotating  shaft  : • » 

Ny  = 0.11  [ 0.5  R0  2 -•  pj0*35 

where  the  Reynolds  number  P ■ is  developed  by  the  shaft 
rotation. 

C07T 

Re  = "v 
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in  which  w is  the  angular  velocity  (rad/sec) 

D is  the  roller  diameter  (m) 

The  average  coefficient  of  forced  convection  to  the  lubri- 
cating oil  within  a rolling  contact  bearing  may  be  approximated 
by, 


a = 0.0986 


*(PT) 


1/3 


using  + for  outer  ring  rotation 
- for  inner  ring  rotation 


in  which  N is  the  bearing  operating  speed  (rpm) 

D is  the  diameter  of  the  rolling  elements  (mm) 
dm  is  the  bearing  pitch  diameter  (mm) 

3 is  the  bearing  contact  angle  (degrees) 

5.2.3  FLUID  FLOW 


The  rate  of  heat  flow  that  is  transferred  from  fluid 
node  i to  fluid  node  j by  fluid  flow  is 


• 

Vj[ j is  the  volume  rate  of  flow  from  i to  j.  It  must  be  observed 
that  the  continuity  of  mass  requires  the  following  equation  to 
be  satisfied 


provided  the  fluid  density  is  constant.  The  summation  should 
be  extended  over  all  nodes  i within  the  fluid  which  have  heat 
exchange  with  node  j by  fluid  flow. 

B.2.4  HEAT  RADIATibN 

The  rate  of  heat  flow  that  is  radiated  to  node  j from 
node  i is  expressed  by 

«Ri,j  - si,i 

where 

Tj  - t.  + 273.16 
Ti  = ti  + 273.16 
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and  the  value  of  the  coefficient  j depends  on  the  geometry 
and  the  emissivity  or  the  absorptivity  of  the  bodies  involved. 

For  radiation  between  large,  parallel  and  adjacent  surfaces 
of  equal  area,  Aj^j  and  emissivity,  ei,j>  ^i,j  is  obtained  from 
the  equation 


where 


a 

9 


6.  . 
1,3 


e . . 

i,3 


the  Stefan-Boltzmann  constant, 


is 


a = 5.76  • 10"8  W/m2  /(degK)4 


For  radiation  between  concentric  spheres  and  coaxial  cylinders 
of  equal  emissivity,  ei,j»  ^i,j  is  given  by  the  equation 


1J 


l + C1- 


ei , j ^ 


& 


i,3 


where  o is  as  above  j is  the  area  of  the  enclosed  body 
and  A*.j^j  is  the  area  of  the  surrounding  body,  i.e.  Aitj  <A*i>j, 

Expressions  for  6j.  j that  are  valid  for  more  complicated 
geometries  or  for  different  emissivities  may  be  found  in  the 
heat  transfer  literature. 


B . 2.5  CALCULATION  OF  AREAS 

In  the  case  of  heat  transfer  in  the  axial  direction  Aj^j  is 
given  by  the  equation  (Fig.  B-3) 

A • • = 2rr  . Ar 
i,3  m 

Referring  to  the  input  instructions,  Section  V,  but  recalling 
L must  be  input  in  mm  not  m. 

t r,  + r-, 

L.,  = r_  = 1 2 

1 m 2 


L2  = Ar  = r2  - r1 
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In  the  case  of  heat  transfer  in  the  radial  direction,  A-j^j 
is  obtained  from  the  expression 


2tt  r 


m 


Hi  4 " 


rm » 
m 


H 


and  similarly  for  the  radiation  term  above 


A*.  . a 2i  r*  H 
m 


h = 2H 

in  which  H is  the  length  of  the  cylindrical  surface;  where  heat 
is  conducted  between  i and  j , rm  is  given  by  the  same  equation 
as  above  (Fig.  B-4'' (a)  );  where  heat  is  convected  between  i and  j, 
rm  is  the  radius  of  the  cylindrical  surface  (Fig.  B-4(b)  );  where 
heat  is  radiated'  between  i and  j , r^  is  the  radius  of  the  enclosed 
cylindrical  surface  and  rm*  the  radius  of  the  surrounding  cylindri- 


9^ 


** ► 

oJ 

. r 

^ ► 

/////$//// 

1 

Fig. B -4(b)  Fig.  B -4(c) 


B.3.1  TRANSIENT  ANALYSIS 

For  the  transient  analysis  all  of  the  data  pertaining  to 
the  node  to  node  heat  transfer  coefficients  must  be  provided  by 
the  input.  Additionally,  .the  volume  and  the  specific  heat  at 
each  node  is  required.  For  metal  nodes  this  input  is  straight- 
forward. However,  when  fluid  flow  is  being  considered  there  is 
no  easy  way  to  approximate  the  fluid  nodal  volume  in  a free  space 
such  as  the  bearing  cavity.  However,  through  use  of  the  Program, 
the  user's  ability  to  make  appropriate  estimates  will  improve. 
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APPENDIX  - C 


S K F COMPUTER  PROGRAM  AT77Y001  "SHABERTH" 


INPUT  FORMAT  FORMS 
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* » 


• * 


T 1 1 1 c C a rd 
Card  Type  TTI 


Tl  1 1c  ( 20 ) 29  30  31  32  33  34  33  30  37  30  39  40  41  12  43  44  4 3 Ah  47 


D 

a 

a 

a 

a 

a 

a 

a 

a 

B 

B 

B 

B 

fl 

a 

a 

B 

S 

1 

3 

3 

a 

i 

a 

a 

a 

a 

a 

9 

a 

a 

1 

a 

a 

a 

a 

a 

a 

a 

a 

g 

1 

1 

9 

3 

a 

3 

a 

m 

3 

B 

E 

2 

B 

1 

B 

a 

B 

B 

a 

a 

E 

S 

g 

a 

a 

s 

a 

s 

0 

a 

E 

a 

a 

a 

5 

D 

B 

s 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

8 

1 

1 

8 

i 

1 

i 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

I 

1 

1 

8 

1 

1 

1 

1 

1 

1 

i 

i 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

■ 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

I 

8 

I 

1 

3 

D 

B 

B 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

■ 

1 

1 

i 

i 

1 

i 

1 

1 

8 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

8 

8 

i 

■ 

1 

1 

1 

1 

1 

1 

■ 

■ 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

■ 

1 

1 

1 

1 

8 

1 

8 

8 

1 

1 

8 

8 

8 

B 

1 

Title  io  be  printed  og  each  page. 
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TITLE  CARD 

CARD  TTPE  TT2  - One  Card 
1FLAG  (1 


)IRJLG<2>  I FLAG (3)  I FLAG (4) 

MnrmT  hfit  ithair 


naaaaraQaDnaBaaaaaaaaaEasaEjaasaaaaEaasasnEEisssECE 

MilWBIilllBMMIMMiHBMi 


Eh**T  XffLft  REVOLUTION* 


IF  M9t  K>M  U*C,  Ki«  ITWITIM 

OEvim  ouxrut  win  w)  flic.  . 

mu  ie  fit  cneu-  Positive  it  holler  clement  **iEr  the  utm  ™»rci| 

a>T  t IMiB.  UTWk  M»  M,  MEHII  DtirfTCR  w BCFIECHOH 

|r  = 1 MO  HERAT  tOH.  FRE-UT  ►'ll  I*  EEtl  »«»■  Ef*l  OR 

fUiR  loop  Positive  mo.  I*  BUS  = If  left  tL*«* , ►«•»«»  EP52  • ielrinc  kfle»- 

J*IUC  OHlT.  Hein'S  FIT  l SICCATIVE  ACCUR4CT  OF  0.0®1  l«  tlO*.  OR  THE  (ClRlMI 

If  = 2 ITEHITIDN  WITH  *0.  «*M  USED.  AWWUR  BtFUXTlOH  It 

FULL  OEKM  HI  :5  - ITEMT I OR  W*S  Than  ? 1 HT*  |M 

OUTWIT.  HEOiTIVE  HO.  LIMIT  = LIRElR  OfFLECT  IOVII 

1 ME**:  FIT  OlVCM  W.  LESS  Twin  5 X 1Q~®  THEM 

ITEHITIOH  HIM  a UEWTIM  H 

Ml  = LIVEN  1 5 ITERATION 

mo.  If  left  him.  f*t*m 

lOCmUCT  OF  D.OOl  IS 

H(4)  E(5)  E{6)  E(7)  E(l>  W“'2) 

! Shaft 

Fit  Call 
Flag 

0 » None 
♦ ■ Smax 
, -Us  e r 

•NPASS,  IFLAGi 5) , Beating  Solution  Level  Flag  * 

0-  Elastic  Contact  Forces  are  calculated,  no  lubrication  or  friction  effects  are  considered 

1 Elastic  Contact  Forces  sre  calculated.  Lubricstiou  and  friction  effects  are  considered  using  rncewaj  control  (ball  bearing) 
or  eplcjrclic  (roller  bearing)  assumptions  to  estimate  rolling  element  and  cage  speeds 

A 

2 Inner  Equilibrium  is  satisfied  considering  only  the  Elastic  Contact  Forces,  Using  the  inter  ring  positions  thus  obtained, 
rolling  element  and  cage  equilibrium  are  determined  considering  friction 

3 Complete  Solution.  The  inner  ring,  rolling  clement  and  cage  equilibrium  Is  determined  considering  all  elastic  and 
friction  forces 

‘IKT  If  one  or  more  of  the  bearing  ring  or  rolling  element  materials  is  other  than  steel,  set  IKT  to  1 and  include  card 
types  Bll  through  B14  for  all  bearings. 


Fit  IMF  ICCVUOT  IT  It 

OIF  I MED  it  cuiNur  •«  01*- 
HCTHtV  CU*R»HCf  BtVlDED 
■T  HOLLER  CLEMENT 
• lire  UR 


Material  Property  Flag 


TOT 


Jti'ar  Iifif  Unfa,  One  Sri  pt?r  hearing.  All  applicable  r arils  hi  lit  rough  lit  l*  for  hearing  1<  f nil  nvril  by  U|  tbrangli  III**  fur  hear  I 
CARfl  TYPE  B1 

hi>M)  mims)  htn-m 


hearing  Type 
BALL 

CYLINDRICAL 

TAPERED 


Hut  I Deg  In  In 
Coluonc  1 


ft  i n -fl  . W_a_L_e_T. 

Designation 111  i n n Life Factors 1 

Inner  Ring 

Ottiar  Ring  Outer 

1 

Inne**  I 

Uricnlai.mii 
Angle  in  De- 
gree! of 
Rolling  Ele- 
ments 0.0  or 
Oland  aeons 
first  *»|i;nru 
is  nil  1 hi; 

Y-ax i ■. 


NC  Nnraclljzero  or  blank,  If  act  to  1,  non  onlfora  rol ler-r acena j profile  gooaetry  ia  Input  on  card  types  DS  and  B6. 
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hearing  data 

CARD  TYPE  U2,  Cue  Card 


DD  (51 M 

T7 

i) HI 

BJ 

HI 

01 mil 

191 BD 

12 

LI 

D 

a 

a 

a 

0 

D 

B 

a 

B 

a 

a 

El 

a 

9 

a 

a 

a 

a 

a 

a 

a 

m 

2 

a 

2 

9 

m 

s 

a 

a 

m 

a 

m 

a 

w 

m 

a 

aa 

m 

m 

9 

m 

m 

a 

s 

9 

9 

m 

m 

a 

a 

a 

a 

m 

iatf 

a 

s 

a 

!* 

m 

a 

a 

a 

a 

a 

■ 

a 

a 

a 

m 

m 

a 

a 

a 

a 

a 

m 

u 

a 

a 

2 

a 

a 

a 

a 

a 

a 

a 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

■ 

1 

1 

1 

1 

1 

1 

1 

1 

i 

1 

1 

1 

1 

1 

1 

i 

1 

i 

i 

i 

II 

I 

1 

i 

1 

1 

1 

1 

i 

1 

1 

1 

l 

1 

1 

1 

i 

1 

1 

1 

1 

1 

1 

1 

1 

i 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

■ 

1 

1 

1 

1 

1 

1 

1 

1 

1 

i 

1 

1 

1 

1 

1 

1 

1 

i 

1 

i 

■ 

i 

1 

n 

i 

1 

i 

1 

1 

1 

1 

i 

1 

1 

1 

1 

1 

1 

1 

i 

1 

1 

1 

1 

1 

1 

1 

1 

1 

i 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

P 

0 

_a 

1 t 
Of 

jlL 

eb 

t 

is. 

Dio.  (a*) 
for  Tapered 
i. Bearing  . . 

No.  Of  Rolling 
Elements 

D i one  trol 
Clearonce  <oa) 

Contact  Angle 
(Deg) 

Flange  Angle 
(Tapered  Only) 

See  Flg.3.9  (Deg) 

Outer  Ring  Speed 
<8P*) 

Axial  Play  for 
Tepered.  Slagle 
Bearing  Only  (on) 

Yt  Cap  Iacl. 

Angle  for 

Tapered . 

' 

BEARING  DATA  CARO  TJPE  One  Card 


jflnonBB555H5nnBHn5Hn33B33B3B33S^SB?P^SBS^^E35BE53B3333S3333Bi!533SES33SB35EB^^^^^^^^M 

liiiiiiiiiiiiiiiiiiiiiiiiiliiiiiiiiiiiiimiiNHinniiiiMiniMiiiiiiiiiiiiiiil 


BALL  BEARING 


&%11  DIa.  (be) 


Roller  DJaneter 


D1 looter  of  Roller 
litrgr  Kml  For 
Tapered 


Roller  Length  Roller  End  Roller  Inclndod 

End  to  End  tee)  Sphere  Rodim  (an)  Angle  (Deg) 


Roller  Crown  | Roller  Flat 


Radial  (on) 


Leogth  (nn) 


BEARING  DATA 
; ARD  TYPE  DdA 


[flQ9b9D9DQHBSQQBQQ3Q53£29S3533!MBHHHlHHHHHHHHHIHHHBHHHHHHIH| 

ImMummmmiimmmiiiiiiHBiiiimimiiimimimiimiiiiiin 

limimmiiMiimmHiimmmmmigimHmmMmiimmiimminl 


Outer  Baccway 
Curvature 
»o/n 


Roller  Btceway 
Effective  Length 
Outer  (on) 


Roller  Raceway 
Effective  Leng 
Inner  (ea) 


Raceway  Crowe 

Raceway  Crown 

Roller  Corner  or 

Roller  Corner  or 

No.  of  Slices 

1 

Sadias  Outer 

Sadias  Inner 

Raceway  Undercut 

Raceway  Undercut 

(do  ) 

(ea) 

Belief  Outer  (no) 

Relief  Inner  (no) 

♦ > 


* * 


Bearing  Data 

Card  Type  B7A,  one  card.  Omit  if  NPASS,  Title  Card  2,  is  aero  or  blank. 


Bearing  Data 

Car*!  tjfjitf  rani*  <1*11  if  NpASS , Tlljr  Curil  2 , |a  i,i?ro  or  li  | « <»  k 

D0( 36 ) riD(37l 


Cifje  Type 

-1 . for  oatcr-rlng 

Rail-Land 

Single  Hall 

Rail-Land 

Rolling  Eleaenl- 

Cage  Weight 

M 

land  riding 

Di aoeter 

Width 

Diaoetral 

Cage  Pocket 

*8. 

O 

+1.  for  loner-ring 
lend  riding 

( sa) 

(no) 

Clearance 

Diaoetral 

VO 

inn) 

Clearance 

(aa) 

0.  for  rolling 
eleaeot  riding 
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Rearing  On  In 

Uni*  farili  D«  omlll  10  IT  operating  clearances  ore  not  to  bn  calculated 
Card  Type  B'J,  one  eard 

rtutni)  mnr.z) ruxr.n) nntMj. 


itncvr.) 


\z  U M*  |fc|7  6 It  |to|Mjte|^|>4tl»|»<>|o|cs|t9|w|£t[^|n|»t|n[HlCT|zg|gf|Ml%l |uj*s| 


w 3^36  5ibch7UokkzkjH^K?k«  f7l44  >9 he  si  bz  rs  m 55bt  y»  w b? ko|6t  w wW  t5 b 47  « h |w  mjwb*  7*  7/  ip l»>p 


Shaft  Fit  Potitire 
ilf  interference 

(aa) 


Hooting  Fitf 
Positive  If 
InlerfercncoCstsi) 


Shaft  Ef fee  life 
[Leng  th 
(am) 


Bearing  Inner 
Ring  Width 
<ma> 


Bearing  On ter 
Bing  Width 
( mm) 


Routing  Effective 
Width 

( mm) 


Card  Type' BIO, 

one  card 

BD:r.7) 

RI>  < r>fl ) 

ni)  ( 60 ) 

111)  (61  ) 

DD (62 ) 

81)  (64 ) 

♦ / 


1,  ■ 


[tearing  Dan 

0**1  corc|  Type l HU  Through  Bid  if  Operating  Clearance*  are  aot  to  be  Calculated 
Card  Type  Nil,  one  card 


Rolling  Element  [Outer  Ring 


Cord  Type  012,  one  cord 
H 1)  ( 7 1 1 Bit  172  1 


I*  «i  i s >on  * s Ratio 


Inner  Ring 


Rolling  Element  | Outer  Bing 


Card  Type  B13,  one  cord 

J BD (77 ) 


HHMgngBgEiaBgBaBflaBB3BagPPsagB3gsig3g3B!SEaEggsE3BE3gB3g3S3SiaBg5BB333B3BasffiaaBaaag 


Density  ‘JnA" 


Inner  Bing 


Ca.'d  Typo  bid,  one  card 
r.OUU)  BO  (02) 


Bolling  Element  { Outer  Ring 


I'icien!  Thermal  Expansion  I/' C 


Inner  Ring  | Roliiug  Element  Outor  Ring 
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BEARING  DATA 

CARD  TTPE  B15,  one  card,  Onlt  if  NPASS,  Title  Card  2,  ia  Zero  or  Blank 


one  ctrd.  Oalt  •*  HP***,  title  Cord  2,  i*  2ero  or  Blm» 
Lubricant  Properties 
DD  (93  )‘ 


Lohri cant 

Designation 

VISi 

Dynastic 

Viscosity  Us) 

0 100aF 

V1S2 

Dyoanie 

Vtscosiiy  (c») 

0 210°F 

RII060 

Density 

gp/co'* 

O is.s'c 

c 

Coefficient  of 
Thermal 

Expaos 1 on 

1/-C 

If  NCODE  It  aero  retd  thit  dtlt 
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Thetftt I Data  - Bearing  Te&pe ra to  rot 


* See  Fig.  3.5  fii  tbo  FI 


Nm  tbit  « uporad  roller  bearing  Flngt  It  fltnge  *1. 


Soil  Oeatlagi  art  nil  cinildand  to  lift  Fliagti 


Thermal  Data  - Individual  Initial  Temperature t 
73,  As  many  cards  as  ooeded«  followed  by  a blank  card 


D 

a 

a 

a 

a 

a 

a 

a 

a 

B 

B 

B 

a 

D 

B 

B 

B 

B 

B 

B 

B 

B 

S 

g 
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B 

B 

g 

g 

g 

a 

B 
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g 

B 

g 

B 

g 

9 

3 

B 

9 

g 

B 

B 

0 

SB 
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3 

a 
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a 

B 

g 

g 
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IS 

a 

a 

a 

3 
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B 

S 

s 

B 

B 

9 

E 

□ 

B 

9 

B 

B 

B 

2 

1 

1 

1 

1 

B 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

Mode 

Number 

Initial 
Tempera  t 

Etc . 
ure 



1 

1 

1 

r i 

1 

1 

i 

1 

1 

; 

Node 

Number 

Initial 
Teape  ra  t u rc 

1 1 

1 



1 

L 1 

1 
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The  rail  Dati  - fiodft  there  Belting  frletloaalviieit  if  Generated 
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Thermal  Data  - Node  heat  capacities,  only  for  transient  calculation! 


T9,  One  Cord/Node 
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Shaft  Data  ■*  Hearing  position  and  support*  X-Z  plane 

S2,  at  nan jr  cards  as  needed.  Card  needed  only  if  any  data  la  different  fron  X-T  plane 
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APPENDIX  - D 


COMPUTER  PROGRAM  SHABERTH 


SAMPLE  OUTPUT 
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The  sample  SHABERTH  output  represents  a steady  state 
thermal  analysis  of  a helicopter  input  pinion  shaft  which 
delivers  turbine  power  to  the  transmission. 

The  system  contains  one  216  size  cylindrical  roller  and  a 
stack  of  four  6216  size  split  inner  ring  ball  bearings.  The 
bearing  heat  rates  were  calculated  executing  SHABERTH  at  the 
NPASS=1  level,  with  two  clearance  change  iterations  permitted 
for  each  each  thermal  iteration.  The  bearings  act  as  heat 
sources  which  supply  heat  to  the  temperature  nodes  representing 
inner  and  outer  rings,  rolling  elements  and  lubricant. 

Additional  heat  sources  include  a gear  mesh  and  a seal  which 
are  specified  by  constant  heat  rate  values  at  input.  The  nodal 
arrangement  and  the  system  equilibrium  temperatures  are  shown 
in  Fig.  D-l. 

All  system  temperatures  were  set  at  an  estimated  value  of 
75°C.  Bearing  performance  is  determined  based  upon  these 
temperatures.  The  bearing  heat  rates  serve  as  input  to  the 
temperature  calculation  scheme  which  in  turn  produces  a new  set 
of  system  temperatures.  These  new  temperatures  affect  bearing 
clearance  and  most  notably  lubricant  viscosity  in  the  next 
calculation  of  shaft  bearing  system  performance  parameters. 

Six  iterations,  calculations  of  bearing  system  performance,  are 
required  before  the  equilibrium  condition  is  achieved.  Equi- 
librium is  satisfied  when  a check  of  all  system  temperatures, 
comparing  two  sequential  iterations  reflect  less  than  a 1°C 
change . 

Although  it  is  not  included,  shaft-bearing  system  output 
data  may  optionally  be  printed  after  each  calculation  of  bearing 
frictional  heat  rates.  This  data  allows  the  user  to  observe 
system  performance  as  a function  of  temperature. 

Since  the  program  input  is  printed,  the  thermal  data  input 
serves  as  a good  example,  to  be  reviewed  when  preparing  thermal 
data  for  another  system. 
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Fig.  D-l  Sample  Problem  Temperature  Node 
Arrangement  and  Steady  State 
Temperatures . 
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00E 

TEMPERATURE 

MOOE 

TEMPERATURE 

NODE 

TEMPERATURE 

NODE 

TEMPERATUR E 

NODE 

TEMPERATURE 

l 

93*103 

2 

94.710 

3 

00.748 

4 

113.494 

5 

117.314 

6 

113.694 

7 

106.149 

a 

102.353 

9 

107.245 

10 

106.375 

11 

95.107 

1 2 

101 .027 

13 

70.735 

14 

82.319 

15 

06.455 

16 

84,951 

17 

66.153 

18 

88.457 

19 

79.535 

20 

84.628 

21 

94.332 

22 

60.004 

23 

67.308 

24 

80.050 

23 

85.233 

26 

94  .609 

27 

80.393 

20 

07.313 

29 

86.495 

30 

79.6i3 

31 

81.989 

32 

87.267 

33 

95.320 

34 

80.343 

35 

86*504 

36 

85.144 

37 

103.561 

38 

04.396 

39 

83.695 

40 

92.442 

41 

101.703 

42 

72.266 

43 

78.765 

44 

90.902 

45 

72.494 

46 

84.750 

47 

94*705 

48 

72.632 

49 

85.078 

50 

95.285 

51 

73.103 

52 

06.022 

53 

95.677 

54 

93.415 

55 

95.130 

56 

93.922 

57 

101.770 

58 

87.757 

KNOWN  BOUNDARY  TEMPCRATURES 

NODE 

TEMPERATURE 

NOOE 

TEMPERATURE 

NODE 

TEMPERATURE 

NODE 

TEMPERATURE 

NODE 

TEMPERATURE 

59 

75.000 

60 

24  .000 

61 

70.000 

131 


SHUEUH/mi  **  TECHNOLOGY  DIVISION  S K F INDUSTRIES  INC,  ♦•SHA9ERTH/ARR* 
POWER  INPUT  MODULE  ‘STEADY  STATE  LUBRICATED  ANALYSIS*  1*50  SHP 


THIS  OATA  SET  CONTAINS  ^ HEARINGS 
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***  SHA8ERTH/A3R  **  TECHNOLOGY  OIVISION  $ K F INDUSTRIES  INC.  **  3HABERTH/A9R  **. 

POWER  IN»UT  MODULE  *STCAOY  STATE  LUBRICATED  ANALYSIS*  1450  SHD 

UNLESS  OTHERWISE  STATED*  LINEAR  DIMENSIONS  ARE  SDEC I F 1 ED  IN  MILLIMETERS*  TEMPERATURES  IN  DEGREES  CENTIGRADE*  FORCES  IN  NEWTONS* 
WEIGHTS  IN  KILOGRAMS*  PRESSURES  AND  ELASTIC  MODULI  IN  NEWTONS  PER  SQUARE  MILLIMETER*  ANGLES  AND  SLOPES  IN  DEGREES  * SURFACE 
ROUGHNESS  IN  MICRONS*  SPEEOS  IN  REVOLUTIONS  PER  MINUTE*  DENSITY  IN  GRAMS  PER  CUBIC  CENTIMETER*  KINEMATIC  VISCOSITY  IN  CENTISTOKES 
AND  THERMAL  CONDUCTIVITY  IN  WATTS  PER  METER  -DEGREE  CENTIGRADE.' 
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• 9 52  50  0 
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RIDING 
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.6  55 

•laoooo 
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RIDING 
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**♦'$  HABERTH/ABR  **  TECHNOLOGY  DIVISION  S K F INDUSTRIES  INC.  **S  HABER  X H / A 9 R *■  * * 
POWER  INPUT  MODULE  * STE AOY  STATE  LU^RICATEO  ANALYSIS*  1*50  SHP 
ROLLING  ELEMENT  0 A T A 

BEARING  NUMBER  ( 1 ) TYPE  - CYLINDRICAL  ROLLER  BEARING 


ROLLER 
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ROLLER  END 

ROLLER  INCL. 

AXIAL  PLAY 

FLANGE  ANGLE 

DIAMETER 
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NO.  OF  AXIAL 

EFF .LENGTH 

FLAT  LENGTH 

CROWN  RAO. 

EFF. LENGTH  FLAT 
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CROWN  RAO. 

laminae 

17.0U0Q 

8.7500 

1524.000 
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8 . 

7500 

1524.000 

10 

TYPE  - BALL  REARING 

OUTER  RACEWAY  CURVATURE 
.515 

TYPE  - BALL  BEARING 

OUTER  RACEWAY  CURVATURE 
.515 

TYPE  - BALL  BEARING 

OUTER  RACEWAY  CURVATURE  INNER  RACEWAY  CURVATURE 


INNER  RACEWAY  CURVATURE 

.520 

INNER  RACEWAY  CURVATURE 
• 520 


BEARING  NUMBER  (2) 

BALL  DIAMETER 
19.0500 


LO 


BEARING  NUMBER  (3) 

BALL  DIAMETER 
19.0500 


SEARING  NUMBER  (4) 

BALL  diameter 


19*0500 


*515 


.520 


4 * 


***  S'  H A 9 E R T H / A 9 R **  TECHNOLOGY  DIVISION  S K F INDUSTRIES  INC*  **SHA9ERTH/ABR  *** 
POWER  INPUT  MODULE  *STEAOY  STATE  LUBRICATED  ANALYSIS*  1450  SHP 
ROLLING  ELEMENT  DATA 

BEARING  NUMBER  (5)  TY»E  - BALL  BEARING 

9 ALL  DIAMETER  OUTER  RACEWAY  CURVATURE 

19.0500  .515 


INNER  RACEWAY  CURVATURE 
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H 

UJ 
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*•*  SHABERTH/A8R  **  TECHNOLOGY  DIVISION  S K F INDUSTRIES  INC.  * ♦ S ’h  A 8 E R T H / A 3 ^ . 

POWER  INPUT  MODULE  *$TEADY  STATE  LUBRICATED  ANALYSIS*  1450  SHP 


SURFACE  DATA 


BEARING 

CLA  ROUGHNESS 

RMS  ASPERITY 

SLOPE 

NUMBER 

OUTER 

INNER 

ROLL.  ELM. 

OUTER 

INNER 

ROLL.  ELM. 

1 

.20 

.20 

.20 

2.000 

2.000 

2..U00 

2 

.15 

.15 

.15 

2.000 

2.000 

2.000 

3 

.15 

.15 

.15 

2.000 

2.000 

2.000 
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. 15 

.15 

.15 
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2.000 

2.000 

5 

.15 

.15 

.15 
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2.000 

2.000 

U B R I 

C 

A N 

T DATA 

SEARING 

DESIGNATION 

KINEMATIC 

VISCOSITY 

DENSITY  AT 

THERMAL  EXPAN. 

thermal 

NUM3ER 

(37.78  Cl 

(9R.89  C> 

(15.56  C) 

COEFFICIENT 

CONDUCT  I VI TY 

1 

M1L-L-7908G 

12.76 

3.20 

.9526 

7. 09-04 

.152 

2 

MIL-L-70OBG 

12.76 

3.20 

.9526 

7.09-04 

.152 

3 

MU.-L-7308G 

12.  76 

3.20 

.9526 

7. 09-04 

.152 

4 

MIL-L-7SD8G 

12.76 

3.2  0 

• 9526 

7.09-04 

.152 

5 

M IL-L- 7R08G 

12.  76 

3.20 

• 9526 

7.09-04 

.152 

U B R I 

C 

A T 

ION  AND  c R 

I C T I 0 N 0 

A T A 

BEARING 

PERCENT  LUBE 

c ILM  REPLENISHMENT 

ASDER I T Y 

NUMBER 

IN  CAVITY 

LAYER 

THICKNESS 

FRICTION 

(ROLL. ELM. 

♦ RACEWAY  ) 

COEFFICIENT 

OUTER 

INNER 

1 

.75 

.3000-03 

.1000-03 

• 10 

2 

• 75 

.1500-02 

.5000-03 

• 10 

3 

1.50 

, 50:10-0? 

• iogn-02 

• 10 

k 

1.50 

•iflOO-02 

• 10 

5 

1*50 

.3000-02 

•1000-02 

• 10 

***  SHABERTH/ABR  **  TECHNOLOGY  DIVISION  S K F INDUSTRIES  INC*  •*  S H A BERTH/A38  ♦** 
POWER  INPUI  MODULE  *STE  ADY  STATE  LUBRICATED  ANALYSIS*  1 «»S0  SHP 
FIT  OATA  AND  M A T F R I AL  PROPERTIES 


BEARING  COLD  FITS  I HR  TIGHT)  EFFECTIVE  WIDTHS 

NUM3ER  SHAFT  HOUSING  SHAFT  INNER  RING  OUTER  RING  HOUSING 

1 .0380  - * 0 0 7 50.0000  26.0000  26.0000  50*0000 

2 .0203  -.0076  40*0030  26*0000  26.0000  40.0000 

3 .0203  -.0076  26.0000  26.0000  26.0000  26.0000 

4 . .0203  -.0076  26.0000  26.0000  26.0000  26.0000 

5 .0203  -.0076  40.0000  26.0000  26.0000  40.0000 


EFFECTIVE  DIAMETERS 


BEARING 

SHAFT 

BEARING 

INNER  RING 

OUTER  RING 

BEARING 

HOUSING 

NUMBER 

I . D. 

BORE 

AVE.  O.D. 

AVE.  I'D. 

O.D. 

a.O. 

1 

64.000 

BO. 000 

93.000 

127.000 

140.000 

184.000 

2 

64.000 

BO  .oon 

99.070 

124. 1 70 

1 40.000 

iR4.nnn 

3 

64.004. 

BO. 000 

99.070 

124.170 

140.000 

184. 000 

4 

64.000 

BO. 000 

99.070 

124.170 

140.000 

184.000 

5 

56.000 

30.000 

99.070 

124.170 

140.000 

184,000 

137 

BEARING  NUMBER  (1) 

SHAFT 

INNER  RING 

ROLL.  ELEM. 

OUTER  RING 

HOUSING 

MODULUS  OF  ELASTICITY 
POISSONS  RATIO 

WEIGHT  density 

COEFF  . 0*  THERMAL .EKP . 

234083.0 
• 3000 

7.306 

.00001224 

204033*0 

.3000 

7.306 

.00001224 

204033.0 
• 3000 

7.806 

.00001224 

2040R3.0 
• 3000 

7.306 

.00001224 

366.9 

.3500 

1.770 

•30002520 

BEARING  NUMBER  (21 

SHAFT 

INNER  RING 

ROLL.  CLEM. 

OUTER  RING 

HOUSING 

modulus  of  elasticity 
pgissons  MTh) 

WEIGHT  DENSITY 

COEFF.  OF  THERMAL  EKP. 

♦ 3000 

7.806 

.00001224 

• 3000 

7.806 

•00001224 

? 04ft 8 SO 
• 30  00 

7*806 

.00001224 

204  Out.  fi 
• 5000 

7.806 

.00001224 

4 1 % 0 * 9 

. 3500 
1.770 
.003025 20 

***  S H H i U H / U R 


***  SHA9ERTH/A9R 

POWER  INPUT  NOPULE 

••  TECHNOLO 

♦STEADY  state 

GY  DIVISION  S K F 

LU3RICATCD  ANALYSIS* 

INDUSTRIES  INC. 

1450  SH° 

• * S H*  A B £ R 

T H / A 9 R *** 

SCARING  NUMBER  (3) 

shact 

inner  RING 

ROLL.  elem. 

OUTER  RING 

MOUSING 

MODULUS  OF  ELASTICITY 

204093.  D 

204083.0 

204083.0 

204083.0 

41369.9 

POISSONS  RATIO 

• 3000 

. 3000 

.30  00 

. 3000 

.3500 

WEIGHT  DENSITY 

7.906 

7.806 

7.806 

7.806 

1 . 770 

COEFF.  OF  THERMAL  EXP* 

*00001224 

.00001224 

•00001224 

• 00001224 

•00002520 

BEARING  NUMBER  <4> 

S 4 A c T 

INNER  RING 

ROLL.  ELEM. 

OUTER  RING 

HOUSING 

MODULUS  OF  ELASTICITY 

204093*0 

204083.0 

204083.0 

204083.0 

41369.9 

POISSONS  RATIO 

.3000 

.3000 

♦ 3000 

. 3000 

.3500 

WEIGHT  DENSITY 

7 .90S 

7.806 

7.806 

7.806 

1,770 

COEFF.  OF  THERMAL  EXP. 

.00001224 

.00001224 

.00001224 

•00001224 

« DO  0 D2  520. 

REARING 

NUM3ER  C 5) 

shaft 

INNER  RING 

ROLL.  ELEM. 

OUTER  RING 

HOUSING 

MOOULUS 

OF  ELASTICITY 

204093.0 

204083.0 

204093.0 

204083.0 

41368-9 

POISSONS  RATIO 

.3000 

• 3000 

• 3000 

.30on- 

. 3530 

WEIGHT 

OENSITY 

7.836 

7.806 

7.806 

7.806 

1.770 

COEFF. 

of  thermal. exp. 

.00001224 

.00001224 

•00001224 

.00001224 

.00002520 

***  s H A 9 E ft  T H-/  ARP  ♦♦  TE  CH  NQLOGY  DIVISION  S K F INDUSTRIES  INC*  **SHABCRTrt/A0R** 

POWER  INPUT  MODULE  *STE  AU  Y 3TATC  LUBRICATED  ANALYSIS*  1450  SH? 

STEADY  STATE  TEMPERATURE  CALCULATION*  ITERATION  LIMIT  10*  ABSOLUTE  ACCURACY  1.00  DEGREES 
INTERMEDIATE  OUTPUT  WILL  RE  OBTAINED 

UNLESS  OTHERWISE  STATEO*  INTERNATIONAL  UNITS  ARE  USED 


NODE  POINTERS 


3RG 

SHAFT 

I.  RING 

I*  RACE 

ROLL  EL. 

0.  RACE 

o*  RING 

HOUSING 

SUL  ^ 

PL  ANGE 

1 

A 

5 

5 

S 

7 

7 

3 

41 

0 

2 

13 

14 

14 

15 

IS 

16 

17 

43 

0 

3 

19 

2 0 

20 

21 

22 

22 

23 

46 

0 

<1 

24 

25 

25 

25 

27 

27 

2B 

49 

0 

5 

31 

32 

32 

33 

34 

34 

35 

52 

0 

NODES  WHERE 

SCARING 

HEAT 

IS  GENERATED 

3RS 

INNER  RACE 

OUTER 

RACE 

CAGE 

ORAG 

flange 

1 

5 

6 

6 

7 

7 

7 

6 

41 

0 

0 

2 

14 

15 

15 

16 

14 

14 

15 

43 

n 

0 

3 

20 

21 

21 

22 

20 

20 

21 

46 

0 

0 

4 

25 

26 

26 

27 

25 

25 

26 

49 

0 

0 

5 

32 

33 

33 

34 

32 

3? 

33 

52 

0 

0 

CONSTANT  GENERATEO  HEATS 

NODE 

GEN*  HEAT 

NODE 

GEN.  HEAT 

node 

GEN.  HEAT 

NODE  GEN.  HEAT 

NODE  GEN.  HEAT 

10 

2500*00 

37 

100.00 

58 

2500.00 

• **  $ H A B.  F.  R T H / A R R **  TECHNOLOGY  DIVISION  S K F INDUSTRIES  INC-  ♦ * $ H 4 3 E ft  T H / A B R • *• 
POWER  INPUT  MODULE  ‘STEADY  STATE  LUBRICATED  ANALYSIS*  1450  Sh’P 
HEAT  TRANSFER  COEFFICIENTS 


H 

■F- 

O 


TY  PE 

INDEX 

COEFFICIENTS 

CONDUCTION 

1 

5.5. 65  W 

CONDUCTION 

2 

46.  72R9 

CONDUCTION 

* 

■in.  32 11 

KUKCEO  CONVECTION 

21 

• 100030 

.570000 

.000000 

• 1 RQOOO-Ol 

34.0000 

• 170000 

.000000 

900.000 

2000.00 

25. 3000 

.183000-01 

219.330 

. 733000-03 

103.434 

.421463*11 

-4 .0.9921 

•ROD 000-03 

-.570000 

FORCED*  CONVECTION 

22 

-300000 

.570000 

.000000 

.250000-01 

4 ft  . 0 0 0 0 

.294000-01 

.220000-04 

1.21000 

1000.00 

.000300 

.000030 

.000000 

.OUUOOO 

18  7.707 

.003003 

.0000:10 

.000000 

.000000 

FORCEO  CONVECTION 

23 

37. 00  00 

FORCED  CONVECTION 

24. 

5000.00 

FORCEO  CONVECTION 

25 

.270000-01 

.800000 

.333000 

.500000-33 

40.0030 

.170000 

.000000 

900.000 

2000.00 

25.3030 

.1*0000-01 

218.330 

.700000-03 

2101.67 

.421460*11 

-4 .01921 

.900000-03 

-.467000 

FORCED  CONVECTION 

26 

.300000 

.570000 

.000000 

.250000-01 

4R.  033  0 

.1  70000 

•000000 

900.000 

2000.00 

25.0000 

.190000-01 

219.330 

.700000-03 

109.315 

.421460*11 

-4 .00921 

.900000-03 

-.570000 

FLUID  FLOW 

41 

28.0000 

FLU  10  FLOW 

42 

42.0000 

FLUID  FLOW 

43 

21.0000 

FLUID  FLOW 

44 

63.0000 

« h 


¥ k. 


• * * jHA  3EKTH  / A R R ♦*  TECHNOLOGY  DIVISION  S * F INDUSTRIES  INC*  **  S HA  SERTH/ARR*** 
POWER  INPUT  MODULE  *S  TE  ADY  STATE  LUBRICATED  ANALYSIS*  l*5G  SH» 

HEAT  TRANSFER  COEFFICIENTS 


T^E 

INDEX 

COEFFICIENTS 

fluid  flow 

*5 

iorj*ono 

fluid  flow 

A 6 

1 A 7 • 00  0 

fluid  flow 

V**»  .000 

-luid  flow 

*8 

136.000 

H* 

-F~ 


***  SHA3ERTH/ABR  **  TECHNOLOGY  DIVISION  S K F INDUSTRIES  INC « **SHA3ERTh/A3R 
POWER  INPUT  HOOULE  ‘STEADY  STATE  LUBRICATED  ANALYSIS*  1*50  $HP 
DESCRIPTION  OP  THE  GEOMETRY  AND  INDICATION  OF  the  TYPES  AND  PATHS  OF  HEAT  TRANSFER 


ALL  LENGTHS 

ARE  IN 

MILLIMETERS.  A 

NES 

ATI  VE 

SIGN  OF 

THE  INOEX  MEANS 

NO  rotational 

SYMMETRY 

TYPE  OF  HEAT 

TR. 

INDEX 

NODE 

NODE 

1ST  LENGTH 

2ND  LENGTH 

3RD  LENGTH 

CONDUCTION 

1 

■BETWEEN 

1 

ANO 

2 

60.0003 

24.C090 

A * . 00  00 

CONDUCTION 

-1 

BETWEEN 

? 

AND 

3 

2 0.0  COO 

149  .9000 

153.3900 

CONDUCTION 

3 

BETWEEN 

p 

AND 

8 

77.0000 

26.0000 

8*. 03  DO 

CONDUCTION 

1 

BETWEEN 

g 

ANO 

12 

96.0000 

ft.OPQO 

70. 0900 

CONDUCTION 

1 

BETWEEN 

12 

AND 

18 

100.0000 

ft. 0000 

74. 0000 

CONDUCTION 

l 

BETWEEN 

18 

AND 

29 

i o * . o o n o 

23.0000 

53.3303 

CONDUCTION 

-1 

3 E TWEEN 

29 

ANO 

30 

20.0000 

149.3000 

150.0030 

CONDUCTION 

1 

BETWEEN 

29 

AND 

3= 

l CA.GOOO 

22 • 0 u 0 0 

2b. 30  GO 

CONDUCTION 

1 

BETWEEN 

2* 

AND 

28 

92.0000 

y 

2ft • 0 0 DO 

22.3300 

CONDUCTION 

1 

BETWEEN 

3ft 

AND 

35 

92.0000 

2ft .0000 

22.0030 

CONDUCTION 

1 

BETWEEN 

36 

ANO 

39 

109. 0900 

22.0003 

26.0090 

CONDUCTION 

1 

between 

*9 

ANO 

3ft 

92.0000 

2ft. 0000 

22.0900 

CONDUCT  ION 

1 

BETWEEN 

3 

AND 

30 

1 70.00CO 

fa.  0 00  0 

276.0030 

CONDUCTION 

l 

Be  tween 

17 

AND 

23 

81.0000 

22.0000 

26. 9030 

CONDUCTION 

1 

BETWEEN 

23 

ANO 

26 

91.0000 

22.0000 

26.0000 

CONDUCTION 

1 

BETWEEN 

28 

ANO 

35 

ai.ooco 

22.0200 

26.3930 

CONDUCTION 

1 

BETWEEN 

35 

AND 

36 

81.0000 

22.0000 

26.0000 

CONDUCTION 

2 

BETWEEN 

8 

and 

7 

70.0000 

26.0000 

10. 0090 

CONDUCTION 

2 

BETWEEN 

17 

AND 

1ft 

7*. 0000 

26.0000 

13.0300 

CONDUCTION 

2 

BETWEEN 

23 

AND 

22 

7*. 0000 

26.0000 

13. 9000 

CONDUCTION 

2 

BETWEEN 

2ft 

ANO 

27 

7*. 0000 

26.0000 

13.0000 

CONDUCTION 

2 

BETWEEN 

35 

AND 

3* 

7 A • 0 00  0 

26.0000 

13.0000 

CONDUCTION 

3 

BETWEEN 

A 

ANO 

5 

90.0000 

26.0000 

7.0000 

CONDUCTION 

3 

BETWEEN 

4 

AMO 

9 

36.0000 

8.0000 

36.9000 

CONDUCTION 

3 

BETWEEN 

9 

AND 

10 

56.0000 

12.0000 

40 * 00 3 0 

4 » 


« « 


***  S H A 8 E R T H / A 8 R **  TECHNOLOGY  DIVISION  S K F INDUSTRIES  INC.  **SH48ERTH/AaR*** 
POWER  INPUT  MODULE  * STEADY  STATE  LUBRICATED  ANALYSIS*  1450  Sh'p 
DESCRIPTION  OP  THE  GEOMETRY  AND  INDICATION  Of  THE  TYPES  AND  PATHS  OF  HEAT  TRANSFER 


ALL  LENGTHS  ARE  IN 

MILLIMCTERS,  A 

NEG 

ATI  VE 

SIGN  OF 

THE  INDEX  MEANS 

NO  rotational 

SYMMETRY 

TY°C  Of  HEAT  TR. 

INDEX 

NODE 

NOOE 

1ST  LENGTH 

2ND  LENGTH 

3RD  LENGTH 

CONDUCT  ION 

3 

BETWEEN 

in 

AND 

11 

5 h * o n o o 

1 4 . 3 u 0 0 

44.0900 

CONDUCTION 

3 

BETWEEN 

5 

AND 

9 

43.0000 

3.0000 

36.9000 

CONDUCTION 

3 

BETWEEN 

11 

AND 

14 

42.0000 

5.3000 

36.9000 

CONDUCTION 

3 

BETWEEN 

11 

and 

13 

38.0000 

7.0000 

36.0300 

CONDUCTION 

3 

BETWEEN 

13 

and 

14 

40.0000 

26.0000 

6.0390 

CONDUCTION 

T 

RETWrEN 

13 

AND 

19 

36.0000 

8.3000 

26.0030 

CONDUCT  ION 

3 

BETWEEN 

14 

AMO 

20 

43.0000 

fc.0000 

26.0000 

CONDUCTION 

3 

BETWEEN 

19 

AND 

20 

40.0000 

26.0000 

3. D 00  0 

CONDUCT: DU 

3 

BETWEEN 

IS 

AND 

22 

56.0000 

8.0000 

26.0*330 

CONDUCTION 

3 

BETWEEN 

2? 

AND 

27 

56.0000 

3.0000 

26.3300 

CONDUCTION 

3 

BETWEEN 

20 

AND 

25 

43.0000 

8.0000 

26. GOOD 

CONDUCTION 

5 

BETWEEN 

1* 

and 

24 

36.0000 

8.0000 

26.0030 

CONDUCTION 

3 

BETWEEN 

?4 

AND 

25 

49.0000 

26.3000 

8.0000 

CONDUCTION 

? 

BETWEEN 

27 

AND 

34 

66.0000 

8.0003 

26.0000 

CONDUCTION 

5 

BETWEEN 

25 

AND 

. .32 

43.0000 

3.0000 

26.0000 

CONDUCTION 

3 

BETWEEN 

2* 

AND 

31 

36.0000 

8.0000 

26.0000 

CONDUCTION 

5 

BETWEEN 

31 

ANO 

32 

43  .0000 

26.0000 

a. 3000 

CONDUCTION 

T 

BETWEEN 

31 

AND 

57 

43.0000 

B.0000 

40.3000 

CONDUCTION 

3 

BETWEEN 

32 

AND 

37 

4b. 0000 

5.3000 

32.0000 

FORCED  CONVECTION 

21 

BETWEEN 

5 

AND 

41 

46.0000 

26.3000 

FORCED  CONVECTION 

21 

BETWEEN 

6 

and 

41 

17.0000 

220.0000 

FOR CEO  CONVECTIun 

21 

BETWEEN 

7 

and 

41 

63.0000 

26.0000 

FOHCZQ  CONVECTION 

21 

BETWEEN 

14 

AND 

43 

46.0000 

26.0000 

FORCED  CONVECTION 

21 

BETWEEN 

15 

ANO 

43 

142.5000 

19.0000 

FORCED  CONVECTION 

21 

BETWEEN 

16 

AND 

43 

53.0000 

26.0000 

• **  SHA9ERTH/ABR  **  TECHNOLOGY  DIVISION  S < F INDUSTRIES  INC 


*SHA  BERTH/A  3.  R** 


0OUER  INPUT  MODULE  *STEaDY  STATE  LUBRICATED  ANALYSIS*  1450  SHP 
DESCRIPTION  Oe  THE  GEOMETRY  AND  INDICATION  Oc  THE  TY°E$  AND  PATHS  OF  HEAT  TRANSFER 

ALL  LENGTHS  ARC  IN  MILLIMETERS*  A NEGATIVE  SIGN  OF  T HE  INDEX  MEANS  NO  ROTATIONAL  SYMMETRY 


I— 1 
•P- 
4=- 


TYPE  OF  HEAT  TR.  INOEX  NODE  NODE  1ST  LENGTH  2ND  LENGTH  3R 3 LENGTH 


FORCED 

CONVECTION 

21 

■AT  TWEEN 

20 

AND 

46 

46.0000 

26.0000 

FORCED 

CONVECTION 

21 

BETWEEN 

21 

AND 

46 

142.5000 

19.0000 

FORCED 

CONVECTION 

21 

BETWEEN 

22 

AND 

45 

63.0000 

26.0000 

FORCED 

CONVECTION 

21 

BETWEEN 

25 

AND 

49 

46.0000 

26.0000 

FORCED 

CONVECTION 

21 

BETWEEN 

26 

AND 

49 

142. 5000 

19.0000 

FORCED 

CONVECTION 

21 

BETWEEN 

2? 

AND 

49 

63.0000 

26.0000 

FORCED 

CONVECTION 

21 

BETWEEN 

32 

AND 

52 

46.0000 

26.0000 

FORCED 

CONVECTION 

21 

BETWEEN 

33 

AND 

52 

142.5000 

19.0000 

FORCED 

CONVECTION 

21 

BETWEEN 

34 

AND 

52 

63.0000 

26.0000 

FORCED 

CONVECTION 

22 

BETWEEN 

1 

AND 

40 

46.0000 

23.0000 

FORCED 

CONVECTION 

22 

BETWEEN 

2 

AND 

40 

60.0000 

20.0000 

FORCED 

CONVECTION 

22 

BETWEEN 

3 

AND 

40 

70.0000 

140.0000 

FORCED 

CONVECTION 

22 

BETWEEN 

ft 

and 

40 

40.0000 

60.0000 

FORCED 

CONVECTION 

22 

BETWEEN 

R 

AND 

40 

60.0000 

40.0000 

FORCED 

CONVECTION 

22 

BETWEEN 

o 

and 

40 

40.0000 

40.0000 

FORCED 

CONVECTION 

22 

BETWEEN 

10 

AND 

40 

66.0000 

100.0000 

FORCED- 

CONVECTION 

22 

BETWEEN 

11 

AND 

40 

42.0000 

30.0000 

FORCED 

CONVECTION 

22 

BETWEEN 

12 

AND 

40 

96.0000 

100.0000 

FORCED 

CONVECTION 

22 

BETWEEN 

17 

AND 

40 

80.0000 

4 Q . 0 0 0 0 

FORCED 

CONVECTION 

22 

BETWEEN 

IB 

AND 

40 

110.0000 

20.0000 

FORCED 

CONVECTION 

22 

BETWEEN 

29 

AND 

40 

120.0000 

10.0000 

FORCED 

CONVECTION 

22 

BETWEEN 

30 

and 

40 

30.0000 

130.0000 

FORCED 

CONVECTION 

23 

BETWEEN 

1 

AND 

60 

46.0000 

23.0000 

FORCEO 

CONVECTION 

25 

between 

42 

AND 

13 

32.0000 

5.0000 

FORCED 

CONVECTION 

25 

BETWEEN 

45 

AND 

19 

32.0000 

5.0000 

4 


***$HABERTH/A8R  **  TECHNOLOGY  OlVISION  S K F INDUSTRIES  INC-  **SHA8ERTH/A3R* 
POWER  INPUT  MODULE  * STE  AD  Y STATE  LUBRICATED  ANALYSIS*  1450  SHP 
OESCR IP  T I ON  OF  THE  GEOMETRY  A NO  INDICATION  OF  THE  TYPES  AND  3 ATH  S OF  HEAT  TRANSFER 

ALL  LENGTHS  ARE  IN  MILLIMETERS*  A NEGATIVE  SIGN  OF  THE  INDEX  MEANS  NO  ROTATIONAL  SYMMETRY 


H* 

\J1 


TYPE  OF  HEAT  TR. 

INDEX 

NODE 

FORCED 

C IN Vf C T T ON 

75 

BETWEEN 

4" 

FORCED 

CONVECTION 

25 

BF TWEEN 

51 

FORCED 

CONVECTION 

26 

BETWEEN 

10 

FORCED 

CONVECTION 

2 3 

BETWEEN 

<•> 

FORCED 

CONVECTION 

'2  3 

BETWEEN 

FORCED 

CONVECTION 

23 

BETWEEN 

2D 

FORCED 

CONVECTION 

23 

BETWEEN 

30 

FORCED 

CONVECTION 

23 

BETWEEN 

36 

FORCED 

CONVECTION 

23 

BETWEEN 

38 

FORCED 

CONVECTION 

23 

BETWEEN 

39 

FORCED 

CONVECTION 

23 

between 

12 

FORCED 

CONVECTION 

?3 

BETWEEN 

R 

FORCED 

CONVEC  T ION 

23 

BETWEEN 

ia 

FORCED 

CONVECTION 

24 

BETWEEN 

l? 

FLUID 

FLOW 

41 

.FROM  . 

61 

FLUID 

FLOW 

41 

FROM 

41 

FLUIO 

FLOW 

42 

FROM 

61 

FLUID 

flow 

42 

-ROM 

61 

FLUID 

flow 

. 42 

FROT 

61 

fluid 

flow 

42 

FROM 

61 

FLUID 

FLOW 

42 

cROM 

42 

FlUI  D 

FLOW 

42 

FROM 

45 

FLUID 

FLOW 

4? 

FROM 

48 

FLUIO 

FLOW 

42 

FROM 

51 

FLUIO 

FLOW 

42 

FROM 

43 

NODE 

1ST  LENGTH 

2NO  LENGTH 

AND 

24 

32.0000 

5.0000 

AND 

31 

52.0000 

5.0000 

AND 

58 

66.0000  * 

100.0000 

AND 

60 

80.0000 

20.0000 

AND 

60 

70.0000 

160.0000 

AND 

60 

120.0000 

i o . n o o o 

AND 

60 

80.0000 

160.0300 

AND 

60 

U C ♦ 0 0 0 0 

30.0000 

AND 

60 

82.0000 

20.0000 

and 

60 

110.0000 

45.0000 

AND 

60 

113.0000 

50.0000 

and 

60 

100.0000 

60.0003 

AND 

60 

50.0000 

25.0000 

AND 

57 

PS. 0000 

100.0000 

T D 

4 1 

(index  41) 

_ 

TO 

57 

(INDEX  47) 

TO 

42 

(INDEX  42) 

TO 

45 

(INDEX  42) 

TO 

48 

(IN^EX  4*> 

TO 

51 

(INDEX  42) 

TO 

43 

( INDEX  42) 

TO 

46 

(INDEX  42) 

TO 

49 

(INDEX  42) 

TO 

52 

( INDEX  42) 

TO 

44 

(INDEX  43) 

3RD  LENGTH 


**•  SH40CRTH/A3R 


TECHNOLOGY  DIVISION  S K F INDUSTRIES  INC 


SHA9ERTH/ABR 


POWER  INPUT  MODULE  - * STE  4 DY  STATE  LUBRICATED  ANALYSIS*  1450  SHP 
DESCRIPTION  Oc  THE  GEOMETRY  ANO  INDICATION  OF  THE  TYPES  AND  n AT  MS  OF  HEAT  TRANSFER 

ALL  LENGTHS  ARE' IN  MILLIMETERS*  A NEGATIVE  SIGN  OF  THE  INDEX  MEANS  NO  ROTATIONAL  SYMMETRY 


TYPE. OF  HEAT  TR. 

INDEX 

NODE 

NODE 

1ST  LENGTH 

2ND  LENGTH 

3RD  LENGTH 

FLUID  FLOW 

4 3 

FROM 

43 

TO 

57 

(INDEX  47) 

fluid  FLOW 

42 

CR9M 

46 

TO 

44 

(INDEX  43) 

fluid  flow 

A ? 

FROM 

46 

TO 

47 

( INDEX  43) 

FLUID  FLOW 

42 

FROM 

49 

TO 

47 

(INDEX  43) 

, 

FLUID  FLOW 

42 

FROM 

49 

TO 

50 

(INDEX  43) 

fluid  flow. 

40 

FROM 

52 

TO 

50 

( INDEX  43) 

FLUID  FLOW 

42 

from 

52 

TO 

53 

(INDEX  43) 

FLUID  FLOW 

42 

FROM 

47 

TO 

55 

( INDEX  45) 

fluid  flow 

43 

FROM 

53 

TO 

54 

(INDEX  44) 

FLUID  FLDW 

42 

FROM 

50 

TO 

54 

(INDEX  *4) 

FLUID  FLOW  • 

44 

FROM 

54 

TO 

55 

(INDEX  45) 

FLUID  FLOW 

4? 

FROM 

44 

TO 

56 

(INDEX  46) 

FLUID  FLOW 

45 

FROM 

55 

TO 

56 

(INDEX  46) 

FLUID  FLOW 

46 

FROM 

56 

TO 

57 

(INDEX  47) 

fluid  flow 

47 

e ROM 

57 

TO 

61 

(INDEX  47) 

fluid  FLOW 

43 

FROM 

61 

TO 

53 

(INDEX  43) 

FLUID  FLOW 

40 

FROM 

58 

TO 

57 

(INDEX  47) 

BEARING  CONDUCTION 

51 

BETWEEN 

6 

AND 

5 

1*0000 

1*0000 

1 .0009 

BEARING  CONDUCTION 

51 

BETWEEN 

6 

AND 

7 

1*0000 

-1.0000 

1.0009 

BEARING  CONDUCTION 

51 

BETWEEN 

15 

AND 

14 

2*0000 

1.0000 

1.39Q0 

bearing  CONDUCTION 

51 

BETWEEN 

15 

AND 

16 

2.0000 

-1.9000 

1 .3000 

BEARING  CONDUCTION 

51 

BETWEEN 

21 

AND 

20 

3*0000 

1.0000 

1.0000 

BEARING  CONDUCTION 

51 

BETWEEN 

21 

AND 

22 

3.0000 

-1.0000 

1.0000 

BEARING  CONDUCTION 

si 

BETWEEN 

26 

AND 

25 

4.0000 

1.0090 

1.9000 

REARING  CONDUCTION 

51 

BETWEEN 

26 

ANO 

27 

4.0000 

-1.0000 

1.0900 

BEARING  CONDUCTION 

51 

BETWEEN 

33 

AND 

32 

5.0000 

1.0000 

1.0009 

BEARING  CONDUCTION 

51 

BETWEEN 

33 

ANO 

54 

5.0000 

-1.0000 

1.0000 

t 4 


•**  SHA8ERTH/A9R  **  TECHNOLOGY  DIVISION  S K F INDUSTRIES  INC.  “SHA3ERTH/ASR**. 

POWER  INPUT  MODULE  ‘STEADY  STATE  LU3RICATED  ANALYSIS*  1450  SHP 
temperature  WAP 

TEMPERATURES  ARE  IN  DEGREES  CELSIUS.  THE  FIRST  SA  TEMPERATURES  ARE  CALCULATED,  THE  OTHERS  ARE  KNOWN 

STEADY  STATE  TEMPCRATURE  CALCULATION,  INITIAL  TEMPERATURES 


CALCULATED  TEMPERATURES 


NODE 

TEMPERATURE 

NODE 

TEMPERA  TURE 

NODE 

temperature 

NODE 

TEMPERATURE 

NODE 

TEMPERATURE 

1 

75.000 

*3 

75.000 

3 

75.000 

4 

75.000 

5 

75.030 

6 

75,000 

7 

75.000 

8 

75.000 

9 

75.000 

10 

75.000 

U 

75.000 

12 

75.000 

13 

75.000 

14 

75.000 

15 

76.030 

15 

75,000 

17 

75.333 

18 

75.000 

19 

75.000 

20 

75.030 

21 

75.000 

72 

75.310 

23 

75.000 

24 

75.000 

25 

75.030 

26 

75.000 

27 

75.030 

2P 

75.000 

29 

75.003 

30 

75.000 

31 

75.000 

32 

75.030 

53 

75.000 

34 

75.000 

35 

/5. 000 

56 

75.000 

37 

75.300 

38 

75.000 

39 

75.000 

40 

75.030 

41 

75.003 

4? 

75.030 

43 

75.000 

44 

75.000 

45 

75.003 

45 

75.UUU 

47 

75.030 

4 a 

75.0  DO 

49 

75.000 

50 

75.030 

51* 

75.000 

52 

75.000 

53 

■75.000 

54 

75.000 

55 

75.000 

56 

75.000 

57 

75.000 

58 

75.000 

KNOWN  SOUNDARY  TEMPERATURES 

NODE  TEMPERATURE  NODE  TEMPERATURE  NODE  TEMPERATURE  NODE  TEMPERATURE  NODE  TEMPERATURE 

51  ' 7R,n(io  En  <to  f>\  7 n • n o ti 


•**  SHABERTH/ABR  **  TECHNOLOGY  DIVISION  S K F INDUSTRIES  INC.  **$HA9ERTH/ABR  *** 
POWER  INPUT  MODULE  ‘STEADY  STATE  LU3RICATEO  ANALYSIS*  1450  SHP 


SHAFT  GEOMETRY.  BEARING  LOCATIONS  AND  SHAFT  l DAD*  DL  ANE  X * Y * 

12  GEOMETRIC  SECTIONS  1 LOAD  SECTIONCSW  5 HEARINGS.  MODULUS  OF  ELASTICITY  = 2.041*05 


POSI- 

INNER 

01  AN. 

OUTER 

91AM. 

3OI  NT 

POINT 

LOAD  INTENSITY 

BEARING 

SEAT 

TION 

LEFT 

RIGHT  ‘ 

LEFT 

RIGHT 

FORCE 

MOMENT 

LEFT  RIGHT 

POS.ERR 

DEFL/FOR 

AN6.ERR 

DEFL/MON 

1 

.0 

.0 

65.4 

• 0 

86.0 

1 

2 

13-0 

55.4 

65.4 

36.0 

96.0 

2 

.000 

0.00 

. 0 Q 0 3 

0.00 

3 

33-0 

85.4 

65.4 

86.0 

R4.0 

3 

A 

40.0 

85.4 

65.4 

84.0 

34  .0 

4 

5 

49.0 

76.0 

76.0 

100.0 

122.0 

5 

6 

58.0 

*2.0 

92.0 

124.0 

124.0 

6 

7 

73.0 

92.0 

92.0 

130.5 

130.5 

3300.0 

318500.0 

7 

8 

104.0 

92.0 

92.0 

139.0 

139.0 

3 

9 

109.0 

96.0 

86.0 

14  1.0 

141  .0 

9 

10 

112.0 

78.0 

78.0 

96.0 

96.0 

id 

11 

119.0 

66.0 

66.0 

84.0 

84.0 

n 

12 

122.0 

64.0 

64.0 

84.0 

84.0 

12 

13 

128.0 

84.0 

64.0 

84.0 

90.0 

13 

14 

149.0 

64.0 

64.0 

90.0 

90.0 

14 

• 000 

0.00 

.0000 

0.00 

15 

174.0 

64.0 

64.0 

90.0 

90.0 

15 

• 000 

0.00 

.0000 

o.oo 

lb 

200.0 

64.0 

64.0 

90.3 

90.0 

16 

. .000 

0.00 

.0000 

0.30 

17 

228.0 

64.0 

64.0 

90.0 

90.0 

17 

.000 

0.00 

.0000 

0.00 

18 

239.0 

64.0 

• 0 

90.0 

.0 

18 

-P- 

OD 


f * 


**  + $HABERfH/A6R  *♦  TECHNOLOGY  DIVISION  S K F INDUSTRIES  INC.  + *SHA3ERTH/a3R  *** 
^OUER  1 N3UT  MODULE  *STC  AflY  STATE  LUBRICATED  ANALYSIS*  1450  SHP 
SHAFT  GEOMETRY,  BEARING  LOCATIONS  AND  SHAFT  LOAD.  PLANE  X - Z. 

12  GEOMETRIC  SECTIONS  1 LOAD  SECTIONS),  5 BEARINGS.  MODULUS  OF  ELASTICITY  = 2.041+05 


THRUST  LOAD  - A. 900+03 


PQSI- 

INNER 

DI4M. 

OUTER 

D I A 4 . 

POINT 

POINT 

LOAD  INTENSITY 

BEARING 

SEAT 

T I ON 

LEFT 

RIGHT 

LEFT 

RIGHT 

FORCE 

MOMENT 

LEFT  RIGHT 

POS .ERR 

DEFL /C0R 

ANG.ERR- OEFL/MO 

1 

.0 

.0 

65.4 

.0 

86.0 

1 

2 

1 3.0 

65.4 

65.  4 

86.0 

86-0 

2 

.OCO 

0.00 

.0000  0.00 

3 

33.0 

65.4 

65.4 

86.0 

3*.0 

S 

4 

40-0 

65.4 

55.4 

ft  4 . .1 

94.0 

4 

5 

<*9.0 

76.0 

76.0 

100.0 

122.0 

5 

6 

5ft. 0 

92.0 

92.0 

124.0 

124.0 

6 

7 

79.0 

92.0 

9 ? • 0 

13  0.5 

no. 5 

11500.0 

7 

9 

104.0 

42. C 

9 ? ■ 0 

159.0 

1 39  . o 

8 

9 

IDS. 0 

96.0 

86. 0 

141.0 

141.0 

4 

10 

112.0 

7ft.  0 

78.0 

96,0 

95.0 

10 

11 

11  7.1 

66.0 

66.0 

*4 .0 

84 . 0 

11 

12 

122.0 

64.0 

64.0 

fl4.  o 

84.0 

12 

13 

129.0 

64.0 

64.0 

84.0 

90.0 

13 

14 

143.0 

64.0 

64.0 

90.0 

90.0 

14 

.000 

0.00 

.0000  0.00 

15 

174.0 

64.0 

64.0 

90.0 

90.0 

15 

• 000 

0.00 

.0000  0.00 

16 

200.0 

64.0 

64.0 

90.0 

90.0 

16 

.000 

0.00 

♦0000  0.00 

17 

226.0 

64.0 

64.0 

90.  o 

90. 0 

17 

.000 

0.00 

•0000  0.00 

ie 

239. D 

64. D 

• 0 

90.0 

.0 

19 

+**  SHABERTH/ABR  + * TECHNOLOGY  OIVISION  S K F INDUSTRIES  INC 


“SHABERTh/AB**** 


POWER  INPUT  MODULE  ‘STEADY  STATE  LUBRICATED  ANALYSIS*  1450  SH*> 

BEARING  SYSTEM  OUTPUT  METRIC  UNITS 


LINEAR  < MM  > AND  ANGULAR  (RADIANS)  DEFLECTIONS  REACTION  FORCES  (N)  AND  MOMENTS  (MM-N) 


3RD  • 

OX 

D Y 

D 2 

GY 

G? 

FX 

FY 

F2 

M < 

M2 

1 

2.581-03 

5*993-04 

9.599-03 

-3.579-05 

5.621-05 

0.000 

459. 

6.902*03 

-319. 

355. 

2 

2*591-0  3 

7.524-03 

1.205-02 

3.314-05 

3.545-05 

2.389*03 

670. 

1.252*03 

3.  5574Q* 

-1.934404 

5 

2. 591-03 

8.570-05 

U 188-0  2 

ft. 389-06 

3.025-05 

2.380*05 

769. 

1.243*05 

3.54  3*0  4 

-2.128404 

A 

2*581-03 

9.109-03 

1.152-02 

1.129-05 

2*709-05 

2.4  10*0  5 

869. 

1.245*03 

3.54  7*04 

-2.405*04 

5 

2. 591-03 

9.780-03 

1.129-02 

1.370-05 

2.472-05 

-2.279*03 

1.079*03 

1.11  7*03 

-3.10**04 

3.071*04 

FATIGUE 

LIFE  (HOURS) 

H/SIGMA 

lube-life 

FACTOR 

MATERI AL 

FACTOR 

BRG. 

0.  RACE 

I.  RACE 

REARING 

0.  RACE 

I.  RACE 

0.  RACE 

1.  RACE 

0.  RACE 

I.  RACE 

1 

2.594405 

2.009*05 

1-221+05. 

1.05 

.914 

• 4 55 

.432 

5.00 

5.00 

2 

7.958*04 

1.575*05 

5.632*04 

1.95 

1.74 

• 639 

♦ 602 

5.00 

5.00 

3 

7.569*04 

1.492+05 

5.350*04 

1 . 76 

1.68 

• 611 

♦ 596 

5.00 

5. 00 

4 

7.3.4  0*0  4 

1.376+05 

5.122*04 

1.75 

1.66 

.607 

.579 

5.00 

5.02 

5 

7.605*04 

1.419*05 

5.279*04 

1.75 

1.61 

.607 

.565 

5.00 

5.0  0 

TEMPERATURES  RELEVANT  TO  BEARING  PERFORMANCE  (DEGREES  CENTIGRADE) 


S • 

SHAFT 

I.  RING 

l.  RACE 

T.  PING. 

ROLL.  EL. 

0.  FLNG. 

0.  RACE 

0.w  RING 

HSG. 

BULK  L 1 10  E 

1 

1 16. 

118. 

118. 

118. 

115. 

107. 

107. 

107. 

1C3. 

102. 

2 

78.7 

92.3 

92.3 

32.3 

86.5 

85.0 

35.0 

85.0 

86.2 

78.8 

3 

79.6 

84.6 

84.6 

84.6 

94.3 

86.1 

88*1 

88.1 

87.4 

84.8 

4 

80.1 

85.3 

85.3 

35.3 

94.6 

88.4 

8 B*  4 

88.4 

87,3 

95.1 

i 

82.0 

•17. 3 

9 7.3 

81,3 

95.3 

88. 4 

98.4 

HP  , 4 

0 ♦_  t . l 

*'••0 

f 


4 i 


* * 


♦**  S H A B E R T H / * 0 R **  TECHNOLOGY  DIVISION  S K F INDUSTRIES  INC.  **SHA9ERTH/A9R  **• 
POWER  INPUT  MODULE  *$ TE  AOY  ST ATE  LUBRICATED  ANALYSIS*  1450  Shp 

BEARING  SYSTEM  OUTPUT  METRIC  UNITS 


FR ICTIONAL 

HE  4 T 5 

ENERATION 

RATE  (WATTS) 

AND  FRICTION  TORQUE 

(N “MM ) 

3RG  • 

0.  RACE 

0.  FLNS 

S.  I.  RACE 

I.  flngs. 

R.E.DRAG 

R .E  .-CAGE 

CAGE-LAND 

TOTAL 

TORQUE 

1 

1.033*03 

o .ooo 

405. 

0.000 

62  6. 

20.5 

10.4 

2.096*03 

1.472*03 

' 2 

266. 

0.0  00 

ISO. 

D.0OC 

293. 

36.6 

39.6 

794. 

557. 

* 

303. 

0.000 

194. 

o .non 

555. 

32.1 

34.3 

1.113*03 

73  5. 

A 

SOI. 

0.000 

115. 

0.030 

553. 

31.5 

34  .1 

1 .115*03 

733. 

5 

301  . 

O.ODO 

190. 

0.000 

549. 

31.1 

33.5 

1.104*03 

776. 

EHD 

rilN  THICKNESS 

* MLM 

REDUCTION  FACTORS  ANO 

HEAT  CONDUCTIVITY 

DATA  FOR 

THE  OUTER 

AND  INNER  ra 

CEWAYS  RESPECTIVELY 

9 RG  . 

FILM  (MICRONS) 

STARV4  TION 

F ACTOR 

THERMAL 

FACTOR 

MENISCUS 

DIST.  (MM) 

CONDUCT! VI 

*TY  (J/DEG.C) 

i 

.268 

.207 

• 998 

.969 

• 9 15 

.932 

.767 

♦ 225 

13.5 

9.18 

2 

.354 

.33? 

.999 

.993 

.967 

• d 69 

1 .16 

.527 

7.1  2 

2.98 

3 

.537 

.321 

1.000 

• 998 

.974 

.375 

1.69 

.813  . 

7.41 

3.05 

4 

.335 

.317 

1 .000 

.993 

.975 

.876 

1.69 

.314  * 

7.46 

3.10 

5 

.335 

.30-7 

1.000 

• 993 

.9  75 

.880 

1.69 

.815 

7.40 

3.08 

c I T PRESSURES  ( N/MM  2 ) 

REARING  clearances 

(MM)  $P 

EEO  GIVING 

ZERO  FIT  PRESSURE 

BRG  • 

SHAFT-COLD* 

3° ER  . 

HSG .-COLO* 

OPER. 

original  change 

OPERATING  SHAFT-INNER 

RING  (RP*> 

1 

10.5 

8.41 

0.000 

0.000 

3. BOO-02  -4  .339-02 

-5.390-03 

3.390*34 

2 

6.63 

3.37 

0.000 

0.000 

2. 053-02  -1.169-02 

1.061-02 

2.113*04 

3 

5.57 

2.42 

0.000 

0.000 

2.053-02  -1.126-02 

1 .104-02 

2.021*04 

4 

5.57 

2.33 

0.000 

0 .000 

2.053-02  -1,154-02 

1.076-02 

2.012*04 

5 

7,67 

2.75 

o.non 

ft  .non 

2.16  3 -0,1  -1,500-02 

1 , A “« 1 * 0 4 

***  S H A 9 E R T H / A 9 R + + T ECHNOLO  S Y OIVISION  S K F INDUSTRIES  INC.  + *"  SHA8ERTH/A3R  * + + 
POWER  INPUT  MODULE  * s r E AD  t state;  LUBRICATED  ANALYST S*' 1450  SHft 

PEARING  SYSTEM  output  metric  units 


LUBRICANT  TEMPERATURES  and  PHYSICAL  PROPERTIES 


LOCATION 

TEMPERATURES 

DENSITY 

VISCOSIT* 

PRESSURE  VISCOSITY 

(DEGREES  C.) 

CGM/CM3) 

KINEMATIC  <CS>  DYNAMIC  (C^) 

COEFFICIENT  (MM2/N) 

8RG. 

I 

OUTER 

107.085 

.8877 

2.82’ 

2.509 

.9817-02 

INNER 

117.977 

.8800 

2.433 

2.141 

.9235-02 

BULK. 

102.361 

.8910 

3.032 

2.702 

.1003-01 

5RG. 

2 

OUTER 

84.962 

.9034 

4.052 

3.660 

.1114-01 

INNER 

82.322 

.9052 

4.255 

3.352 

.1131-01 

VJLK 

78.787 

.1078 

4.552 

4.132 

.1155-01 

BRG. 

3 

OUTER 

88.014 

• 9312 

3.831 

3.452 

.1094-01 

INNER 

84 .621 

• 9 036 

4.077 

3*684 

.1116-01 

BULK 

84.752 

.9035 

4.068 

3.675 

♦1115-01 

BRG. 

4 

OUTER 

8%. 408 

.9009 

3.810 

3.432 

.1092-01 

INNER 

85.257 

.9032 

4.030 

3.640 

.1112-01 

BULK 

85.083 

.9035 

4.043 

3.652 

.1113-01 

8 R G » 

5 

OUTER 

88.360 

.9010 

3.813 

3.435* 

.1092-01 

INNER 

87.273 

.9017 

3.887 

3.505 

.1099-01 

-ULK 

88.050 

.9026 

3.974 

3.587 

.1107-01 

CAGE  DATA  METRIC  UN  I TS 

CAGE  RAIL  - R I NS  LAND  DATA  CAtir.  SPCED  DATA 


TORDUE 

HEAT  RATE 

SEP. FORCE  ECCENTRICITY 

EPICYCLIC 

SPEED 

CALCULATED 

SPEED 

CALC/E^IC 

CAGE/SHAFT 

HRS  • 

(MM-N) 

(WATTS) 

(N.’W  TONS) 

RATIO 

< RAD/SEC) 

( R?M  ) 

(RAD/SEC) 

( RPM  ) 

RATIO 

RATIO 

1 

-17.2 

10.4 

0.000 

0.000 

802. 

5. 749.03 

602. 

5.749+03 

1 .000 

.423 

2 

49.1 

35.6 

1.618-03 

1.000-02 

638. 

6.093+03 

638  • 

6.093+33 

1.00 

.448 

X 

43.7 

34.3 

1.439-03 

1.000-02 

638. 

6.097+03 

633. 

6.097+03 

1.00 

.448 

4 

43.4 

34.1 

1.430-03 

1.000-02 

633  • 

6.095+03 

6 3 8. 

6.095+03 

1.00 

.443 

5 

42.6 

33.5 

1.405-03 

1.000-02 

638. 

6.095+03 

638. 

6.095+03 

1.00 

• 443 

. # * 


♦ •♦SHA8ERTH/ABR 


♦ * 


TECHNOLOGY  division  $ k f industries  INC 


♦SHABERTH/A9R 


>0 \4K*  IAIPUT  MODULE;  ‘STEADY  STATE  LU4R1CATC0  ANALYSTS*  1450  SH? 


rolling  element  output  for  scaring  number  i metric  units 


Azimuth  INSULAR  SPEEDS  (RAOI AN3/SEC0N1I  SPEED  VECTOR  ANGLES  (DEGREES) 


angle  OES.) 

7X 

WY 

wz 

TO  1 AL 

0991 TAL 

TAN-1  I UY/yX) 

TAN-1  ( VZ/UX) 

*00 

-4477.613 

-.115 

.000 

4497.618 

6 3 2.043 

-160.00 

160.00 

22.50 

-4497,613 

-.143 

• 000 

4497,618 

602.043 

-160.00 

160.00 

45.00 

-4447.619 

-.093 

.000 

4497.618 

602.043 

-180.00 

190.00 

67.50 

*449  7.618 

-.076 

.000 

4497.618 

602.043 

-190.00 

HO. 03 

40.00 

-449 7.613 

-.069 

.000 

4497.613 

6 0 2.043 

-180.00 

193.00 

112.50 

-44  *7.618 

-.017 

.0  00 

4497.61 8 

6 3 2.043 

-180  .00 

190.00 

135.00 

-4497.613 

.021 

.000 

44*97.61  8 

6 0 2.043 

190.00 

180.03 

157.50 

-4417.613 

• 08  3 

.000 

4497.618 

602.043 

180.00 

190.00 

180.00 

-4497.618 

• 000 

.000 

4497.518 

602.043 

180.00 

190.00 

202.50 

-4497.613 

.000 

• 0 00 

4497.613 

602.0*3 

160. 00 

180.00 

225.00 

-4497.618 

.000 

.000 

4497.618 

602.043 

180.00 

190.00 

247.50 

-4497.618 

.000 

• 000 

4497. 6ie 

602.043 

180. OU 

190.00 

270.00 

-4497.618 

• 000 

• 000 

4497.618 

602.043 

160.00 

180.00 

292.50 

-4497*618 

.000 

• 000 

4497.618 

602.043 

180.00 

193.00 

315.00 

-4497. 6l9 

.000 

• 000 

4497.618 

602.043 

190.00 

180.00 

337.50 

-4497.618 

• ODD 

• 000 

4497.618 

602.043 

180.00 

180.00 

SHA9ERTH/49R  “ TECHNOLOGY  DIVISION  S K F INDUSTRIES  INC 
POWER  INPUT  MODULE  ‘STEADY  STATE  LUBRICATED  ANALYSIS*  1450-  SHP 


• •SHA3ERTH/AVR 


ROLLING  ELEMENT  OUTPUT  FOR  SEARING  NUMBER  1 METRIC  UNITS 


AZI MUTH 

normal 

FORCES  (NEWTONS) 

• 

HZ  STRESS 

<N/MM“2> 

LOAD  RATIO 

OASP/OTOT 

CONTACT 

ANGLES  < DEG. ) 

ANGLE  (DCS.) 

CAGE 

■JUTER 

INNER 

OUTER 

INNER 

OUTER 

INNER 

OUTER 

INNER 

• 90 

.000 

91  0.450 

1 39.426 

549.366 

292.364 

.3104 

.6986 

.00 

.00 

22-50 

.000 

1344.194 

721.827 

716.860 

611.850 

. 2746 

.4611 

.00 

.00 

4 5*00 

• 000 

1 127.221 

1 255.971 

826.604 

802.054 

• 2 55  6 

. 3M5 

.00 

.20 

67,50 

. 000 

2219.878 

1629.  .549 

386 .688 

906.016 

.2465 

.3587 

.00 

.00 

90,00 

-.000 

2404 .8*9 

. 17  33.873 

902.574 

925.105 

.2442 

. 35  38 

.00 

• 00 

112,50 

-.000 

220«.612 

1537. 00R 

R67.997 

868.406 

.2492 

• 3689 

.00 

.00 

135.00 

-.000 

17C  3.491 

1092. 359 

792.25  8 

753.966 

• 361  3 

.4113 

.00 

.00 

157.50 

-.000 

1 20  142 

538.3*5 

663.845 

524.562 

.2849 

.5061 

• 00 

.00 

1R0.00 

-•000 

671.13? 

.000 

489.194 

.000 

. 3258 

• 0000 

.00 

.00 

202  .50 

-.000 

671.050 

• 000 

489.162 

.000 

. 325R 

.0000 

.00 

.00 

225.00 

-.000 

671.050 

.000 

489.162 

• 000 

.3258 

.0000 

.00 

.00 

247.50 

-.000 

671.050 

* 9 90 

489.162 

.000 

• 3258 

.0000 

.00 

.00 

270,00 

-.000 

671.050 

.000 

489. 162 

• 000 

. 3?53 

• 0 000 

.00 

.no 

292,50 

-.000 

671.050 

• 000 

489.162 

• 000 

.3258 

.0000 

.00 

.00 

315.00 

-.000 

671.050 

• 000 

489.162 

• 000 

.3258 

• 0 0 00 

• 00 

.00 

337.50 

-.000 

671.050 

• 000 

489.162 

.030 

.3258 

• .3030 

• 0 0 

.00 

* • 


a 


***  SHABERTH/A9*  **  TECH  N0LO3 Y DIVISION  S K F INDUSTRIES  INC*  **SHABERT4/A9R*** 
POWER  INPUT  MODULE  *STE  A DY  STATE  LUBRICATED  ANALYSIS*  1450  SHP 


ROLLING  ELEMENT  OUTPUT  FOR  BEARING  NUMBER  2 METRIC  UNITS 


AZ?MUTH  ANGULAR  SPEEDS  ( R AD  I ANS / SECO NO ) SPEEO  VECTOR  ANGLES  (DEGREES) 


1LE  (DEG*  > 

4< 

WY 

9Z 

TOTAL 

C95ITAL 

tan-i (uy/wx) 

TAN-1 (NZ/WK) 

• 00 

-4157.236 

721.758 

-4.559 

4 >20.323 

624.693 

170.04 

-1  71.94 

24.00 

-4  H 92 . 2 31 

915.425 

-5.046 

4172.380 

619.054 

163.72 

-171.9  5 

43.00 

-4057.495 

966.523 

-5.327 

4149.195 

614.807 

167.14 

-179  .92 

72.00 

-4055.935 

971. 385 

-5.361 

4148.672 

614.799 

107.87 

-179.12 

96.00 

'-4037.756 

832.495 

-5. 144 

4171.667 

617.959 

163.41 

-1  71.  9 3 

120.00 

-4150.739 

753.940 

-4.703 

4218.716 

624.494 

169.71 

-1 79.9  3 

1 A A * •'  0 

-4235.714 

652.725 

-4.139 

4295.714 

633.124 

171.24 

-1 75  .14 

168.90 

-4326.513 

541.307 

-3.545 

4361.309 

644.694 

172.76 

-179.95 

192.00 

-4406.380 

462.336 

-3.027 

4430.575 

654.648 

174.01 

-171.16 

216.00 

-4461.453 

40  2. 039 

-2.660 

4479.512 

661.715 

174.85 

-171.97 

240.00 

-4481.653 

376.364 

-2.504 

449  7.486 

664.511 

175.19 

-1 71.97 

264.00 

-4463.094 

399.766 

-2.580 

4480.979 

661.904 

175.01 

-171.97 

288.00 

-4411.415 

439.990 

-2.832 

4433.304 

654.192 

174.30 

-17?. 96 

312.00 

-4333.271 

521.645 

-3.365 

4364.557 

645.034 

173.14 

-1 79.96 

336.00 

-4242.902 

624.210 

-3.959 

4288.575 

634.236 

171.63 

-179.95 

**#$HA8ERTH/ABR  **  TECHNOLOGY -DIVISION  S K F INDUSTRIES  INC 


**SHA3ERTH/A8R  **♦ 


POWER  INPUT  NOOULE  *STEI9Y  STATE  LU2R 1C  A TED  ANALYSIS*  1450  SHP 


ROLLING  ELEMENT  OUTPUT 


FOR  BEARING  NUMBER  2 METRIC  UNITS 


VJ1 

G\ 


AZIMUTH 

NORMAL 

FORCES  (NEWTONS) 

HZ  STRESS 

<N/Mtt**2  ) 

LOAD  RATIO 

DAS°/QT3T 

CONTACT 

ANGLES  (DEG 

ANGLE  ( DEG*  > 

CAGE 

OUTER 

INNER 

OUTER 

INNER 

OUTER 

INNER 

OUTER 

INNER 

• 00 

-1.137 

939,437 

366.464 

1 1C8. 148 

994.243 

.0403 

.05*0 

11.67 

31.24 

21*  *00 

-.4*3 

1033.7*50 

475.546 

1144.056 

1034.464 

.0418 

• 05  73 

13.22 

29.85 

43*00 

-.141 

1106.363 

555.072 

1170.239 

1141.828 

.0428 

.0572 

14.13 

29.11 

72*00 

.133 

1111.429 

580.47ft 

1172*022 

1143.523 

.0426 

.0571 

14.22 

29.18 

’6.00 

.470 

1046.267 

489.586 

1148.654 

1094.272 

.0413 

.0577 

13.  “9 

29.97 

120  *00 

1.103 

964.615 

382.948 

1114.084 

10  0 B . 9 35 

.0404 

.0533 

12.07 

31.42 

144.00 

2.636 

87*. 690 

286.328 

1 083.728 

915.740 

• 0 385 

,0604 

1C. 27 

3 3.14 

16(5.00 

2.239 

930. 9 71 

213.933 

1063.744 

831.019 

.0366 

.0624 

8.49 

34.93 

192.00 

.920 

805.959 

165.379 

1052.963 

762.625 

. 0350 

.0644 

7.03 

36.33 

216.00 

• 36 1 

794.514 

133.197 

1047.954 

718.319 

. 0 335 

.0680 

8.04 

37.23 

240.00 

• 006 

790.222 

127.390 

1046.064 

699.033 

. 0330 

. 0 666 

5.64 

3 7.50 

264.00 

-.347 

790.853 

133.796 

1046.342 

710.610 

. 0335 

.0661 

5.95 

37.14 

288.00 

-.996 

798.670 

157.610 

1049.778 

750.490 

.0349 

.0647 

6.68 

36.13 

312*00 

-2.196 

820.347 

202*464 

1059.191 

815.831 

.0365 

.0  627 

8.05 

34.73 

336.00 

-2.792 

864.848 

270.781 

1078.006 

898.856 

• 0384 

.0807 

9.31 

32.99 

■* 


< U- 


* 


***  S H A 0 E ft  T H t A 3 R **  TECHNOLOGY  DIVISION  S K F INDUSTRIES  INC, 
POVER  INPUT  MODULE  *STE  AH  Y STATE  L U3RI C A TED  ANALYSIS*  1*50  SH° 


**$HA3ERT  H ./  A B R **< 


R '0  L L I N G 


ELEMENT 


OUTPUT 


FOR  SEARING  NUMBER  3 METRIC  UNITS 


AZIMUTH 


ANGULAR  SPEEDS  < R A D I ANS / SECOND ) 


SPEED  VECTOR  ANGLES  (DEGREES)  — 


bJ 

VJl 


LE  (DE3.) 

vx 

'•JY 

UZ 

TOTAL 

ORBITAL 

TAN-1 (UY/WX) 

TAN- 1 < MZ/WX1 

• 00 

-4145.807 

745.555 

-4.649 

4212.314 

623.5  1 2 

169.81 

-179.94 

24.00 

-4082.302 

330.625 

-5.126 

4166.442 

617.180 

168.50 

-1 79.93 

48.00 

-4051.950 

375.709 

-5.380 

4145.503 

614.344 

167.80 

-179.92 

72.00 

-4055.339 

374.166 

-5.374 

4143.490 

614.775 

167.54 

-174.92 

96.00 

-4092.940 

326. 370 

-5.111 

4175.532 

618. 496 

168.59 

-179.93 

120.00 

-4162.U2 

739.763 

-4.629 

4227.545 

625.699 

169.92 

-179.94 

144.00 

-4251.739 

633.001 

-4.024 

4298.603 

635.747 

171.53 

-179.95 

168.00 

-4345.053 

527.339 

-3.412 

4376.945 

646.928 

173.08 

-1  79.95 

192.00 

-4424.527 

440.488 

-2  .894 

4446.400 

655.924 

174.31 

-179.96 

21*. 00 

-4476.445 

333.461 

-2.545 

4492.840 

663.640 

175.10 

-179.97 

240.00 

-4491 .909 

363.266 

-2.413 

4506.575 

665.623 

175.38 

-179.97 

2*4.00 

-4488.208 

382.808 

-2.536 

4484.577 

662.419 

175.10 

-1  79.97 

2BR.00 

-4409.355 

440.545 

-2.384 

4431.318 

654. 70G 

174.29 

-179.96 

312.00 

-4325.977 

529.303 

-3.410 

4358.239 

644.176 

173.02 

-179.95 

336.00 

-4232. 191 

637.509 

-4.03b 

4279.929 

633.013 

171.43 

-179.95 

*•»  S MB  E R T H./  H R **  TECHNOLOGY  OIVISION  S K F INDUSTRIES  INC.  **$HA3ERTH/A8R  *** 
POWER  INPUT  NODULE  * S TE  AOY  STATE  LUBRICATED  ANALYSIS*  1450  $H® 


ROLLINS  ELEMENT  OUTPUT  .FOR  REARING  NUMBER  3 METRIC  UNITS 


AZIMUTH 

NORMAL 

FORCES  (NEWTONS  > 

HZ  STRESS 

(N/MM**2  > 

LO AO  RATIO 

Q ASP/1  TOT 

CONTACT 

angles  <oe; 

ANGLE  (DEG.) 

CAGE 

OUTER 

INNER 

OUTER 

INNER 

OUTER 

INNER 

OUTER 

INNER 

.00 

-.881 

955.357 

383.103 

1113.595 

1009.074 

. 0 497 

.0669 

11.95 

- 31 .02 

24.00 

-.374 

1051.870 

495.551 

1130.623 

1099.4  70 

.0514 

.0655 

13.48 

29  .67 

48.00 

-.089 

* 1119.976 

569.970 

1175.019 

1151.953 

.0522 

.0  <*4  7 

14.29 

29.02 

73.00 

.183 

1113.026 

562.317 

1172.583 

1146.775 

.0521 

. 0 6'*8 

14.R6 

2 .9 . 1 7 

96. QQ 

.469 

1036.771 

478.141 

1145.169 

1036.453 

.0511 

.0656 

13.38 

30.11 

120.00 

. 1.099 

941.507 

367.442 

1109.961 

995  .'1  32 

.049? 

. 0672 

11.32 

31.66 

144,00 

2.488 

867.491 

270.842 

1 0 79. 1 C4 

898.924 

.0467 

.0694 

9.93 

33.49 

168.00 

1.71-0 

823.458 

201.421 

1060.529 

814.428 

.0441 

.0719 

8. 12 

35.24 

n?.oo 

.720 

811.272 

155.231 

1050.917 

746.696 

.0419 

. 0 74  3 

6.67 

36.61 

216.00 

.270 

791 .71 4 

129.515 

1046.722 

702.950 

.0406 

.0760 

5.74 

37.43 

240.00 

-.036 

788.503 

122.031 

1045.305 

689.140 

.0402 

.0786 

5.42 

37.63 

\J1 

rn 

2&4.00 

-.360 

739.969 

131.479 

1045.952 

706.435 

.0403 

.0753 

5.74 

37.17 

288.00 

-.893 

799.006 

15R.-279 

1049.926 

751.552 

.0424 

.0740 

6.69 

* 36.12 

312.00 

-2.133 

823.302 

207.423 

1060.462 

822.439 

.0447 

.0715 

8.18 

34.59 

336.00 

-2.123 

972.427 

281.059 

1081.148 

910.088 

.0473 

.0690 

10.04 

32.78 

» t 


c..  i 


***  SHA3ERTH/A3R  *•  TECHNOLOGY  DIVISION  S K F INDUSTRIES  INC.  **SHA8ERTH/A3R  *** 
POWER  INPUT  MODULE  ‘STEADY  STATE  LUBRICATED  ANALYSIS*  1450  S H=> 


ROLLING  ELEMENT  OUTPUT  FOR  REARING  NUMBER  4 METRIC  UNITS 


AZIMUTH  ANGULAR  SPEEDS  ( R A D I ANS /SECOND ) SPEED  VECTOR  ANGLES  (DEGREES) 


LE  (QES.) 

VK 

JY 

WZ 

TOTAL 

CR9ITAL 

TAN-1 (WY/WX) 

T AN  - 1 ( VIZ  / WX  ) 

• 00 

“4131.371 

76  5.  *.31 

-4.761 

420i .682 

622.041 

169.50 

-1 7 7.  93 

24,00 

-4071 .246 

847.645 

-5.222 

4158.594 

616.1 16 

163.24 

-179.93 

4 8.00 

-4044.746 

837,042 

-5.444 

4140.375 

613.727 

167.63 

-179.92 

72.00 

-4052.526 

879.221 

-5.403 

414^,809 

614.556 

167.76 

-179.92 

96.00 

-4094.B82 

324. 709 

-5.103 

4177.108 

613.719 

163.61 

-17  1.93 

12  0.00 

-4168 .807 

731.  30  3 

-4  .585 

4232.570 

626.432 

170.04 

-1 79.94 

144.00 

-4251.743 

621.306 

-3.957 

4 3 0 6.796 

636.91  0 

171.71 

-1  7*3 .93 

163.00 

-43  :3.069 

514.741 

-3.357 

4336.377 

649.280 

173.26 

-179.96 

192.00 

-4434 .212 

479.304 

-2.825 

4454.946 

658.157 

174.47 

-175.96 

216.00 

-4432.891 

375.578 

-2.496 

4498.597 

664.472 

175.21 

-179.97 

240  .00 

-4493.250 

360.741 

-2.402 

4507.709 

665.785 

175.41 

-179.97 

264.00 

-4454.109 

386. 245 

-2.556 

4400.793 

661.868 

175.05 

-179.97 

2S3.00 

-4403.233 

450.012 

-2.941 

4423.235 

653.532 

174.16 

-179  .96 

312.00 

-4312.725 

544.691 

-3.500 

4346.989 

642.566 

172.90 

-179.95 

336.00 

-4217.164 

656.545 

-4.145 

4267.967 

631.325 

171.15 

-173.94 

♦**  SHA3ERTH/A9R  **  TECHNOLOGY  DIVISION  $ K F INDUSTRIES  INC.  * * S*  H A 3ER  TH/A  B R *.* 
POJER  INPUT  M30ULE  *STEADY  STATE  LUBRICATED  ANALYSIS*  1450  SH^ 


ROLLING  ELEMENT  OUTPUT  FOR  BEARING  NUMBER  A METRIC  UNITS 


AZIMUTH 

NO  RM  AL 

FORCES  (NEUTQNS) 

HZ  STRESS 

( N/ SM**  2 ) 

LOAD  RATIO 

QASP/OTOT 

CONTACT 

ANGLES  <DE( 

angle  <de$*> 

CAGE 

OUTER 

INNER 

OUTER 

INNER 

OUTER 

INNER 

OUTER 

INNER 

.00 

-.7  75 

472.367 

4 0 5.433 

1 120. 94  4 

1023.352 

.0511 

.053? 

12.30 

30.73 

24.00 

-.32b 

1075.177 

521.551 

1154.133 

1113.364 

.0527 

.0473 

13.73 

29.44 

43.00 

-.057 

1139.245 

591.409 

1181. 719 

1166.219 

. 0534 

.0671 

14.50 

28  .9  3 

7?. 30 

.143 

J 120.253 

570.595 

1175.116 

M^2.3  74 

-0532 

.067  ? 

14.35 

2 9,13 

94.30 

.514 

1033.236 

474*339 

1143.466 

1033.545 

• 0520 

.0633 

1 3.35 

30.17 

120.00 

1.135 

334  .246 

358.863 

1106.103 

937.330 

.0500 

.0  702 

11.47 

31.91 

144.00 

2.47fl 

361.300 

2 6 2.340 

1 076.531 

8.39.418 

.0473 

. 0 726 

9.73 

33.68 

isa.oo 

1.491 

819.498 

194.393 

1054.826 

804.343 

.0446 

.0753 

7.90 

35.42 

142.00 

.638 

799.179 

150.072 

1050.302 

738.330 

. 0425 

.0779 

6.49 

36.76 

216.00 

.225 

790.456 

127.007 

1 046.344 

693.332 

.041  2 

.0796 

5.62 

37.52 

240.00 

-.073 

788.294 

120.474 

1045.213 

696.197 

.0410 

.0900 

5.38 

37.63 

264.00 

-.407 

790.534 

132.496 

1046.224 

703.302 

.0417 

.0790 

5.80 

37.09 

283.00 

-.931 

301 .382 

163.19? 

1050.965 

759.233 

.0435 

.0769 

6.85 

35.94 

312.00 

-2.244 

32 H. 79 3 

215.660 

1 062,81  4 

833.184 

.0459 

.0742 

A. 44 

34.35 

333.00 

-1.827 

883.910 

296.330 

1035.870 

926.231 

.0437 

.0715 

10.33 

32 .50 

e*  - 


SH43ERT4/A3R  *♦ 


TECHNOLOGY  OIVISION  S K c I NOUST  R I £ $ ’ I NC 


**SHAbERTM/ABR‘* 


POWER  INPUT 


ROLLING 


4Z  I MiJT H 
ANGLE  OE3. 

.3 : 

?£» . DO 
4 8 • 0 0 

72.00 

96.00 
120.00 
1**A  . 0 0 
163.  Du 

192.00 

216.00 
2 A 0 . 0 0 
264  .00 
2^8.00 

312.00 

536.00 


MODULE  * STE  AD Y STATE  LUBRICATED  ANALYSIS*  1450  SHO 


ELEMENT  OUTPUT  FOR  HEARING  NUMBER  5 METRIC  UNITS 


ANGULAR  SPEEDS  (RAO  IANS/ SECOND) 


u 

VY 

« 1 T 

A 4 

TOTAL 

-4109.180 

- 79  6.36  7 

4.873 

4183*787 

-4  3 58 i 059 

-454.67  > 

5.252 

4147.096 

-4042.335 

-R7'-i  .?A? 

5.559 

4155.796 

-4061.716 

-p.a  5.  89  5 

5.1  87 

4143.461 

'4116.435 

- 765. 722 

A. 747 

4 187.051 

-APQ0.0A1 

-S54, 650 

4.1  16 

4250.756 

-4,297.467 

-535. 014 

3.425 

4330.643 

-4390.275 

-429.555 

2.799 

4411.238 

-4462.078 

-553.696 

2.538 

4476.075 

-4500.'40a 

-316.919 

2.112 

4511.553 

-4498.993 

-322.719 

2.150 

4510.463 

-4457.1 74 

-370.712 

2.449 

4472.565 

-4392. R06 

-45 A .36  2 

2.958 

4406.296 

-4289*106 

-56  3 . 7 A 7 

3.605 

4325.011 

-4190.714 

-681.961 

4.284 

4245.842 

SPEED  VECTOR  ANGLES  (DEGREES) 


ORBITAL 

TAN-1  (IJY/Wxj 

TAfJ-1  (WZ/WO 

619.509 

-169.16 
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